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A class of Interactive control systenis is derived in Tart I by 
ge. ./alizing interactive manipulator control systems. Tasks of inter- 
active control systems can be represented as a network of a finite set 
of actions which have specific operational characteristics and specific 
resource requirements, and which are of limited duration. This has en- 
abled the de;:umposition of the overall control algorithm into simul- 
taneously and asynchronously. 

The general objective for development of Part It is to evaluate 
the performance benefits of sensor-referenced and computer-aided control 
of Tuanipulators in a complex environment. This report represents the 
first phase of the CURV Arm Control System (CACS) software development, 
and gives the basic features of the control algorithms and their soft- 
ware Implementation. 

Part III investigates the problem of finding an optimal solution 
for a production scheduling problem that will be be easy to implement In 
practical situations. The results show that the optimal solution Is 
very easy to Implement In real life problems as the jobs have to be 
arranged according to mcnotinic increasing processing times. 

Part IV is an initial Investigation and Is the first In a series 
leading to a fully-developed model, ROBECON, which may be used for speci- 
fying the economic consequences of robot systems acquisitions. The model 
will be computer-based and user interactive 
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AISTXACT 

A claaa of IncaractiTO control ayoCMo ia darlTod 
by ganaralixing intaractiva nanipulator control oyatana* 
Tha ganaral atructural propartiaa of auch ayataan ara 
diacuaaad and an appropriata ganaral aofewara iaplo- 
nantation ia propoaad. Thia ia baaad on tha fact that 
taaka of iLcaractira control ayatana can ba rapr^ 
aantod aa a nacwork of a finita aat of actiona idiich 
hara apacific oparational charactariatica and apacific 
raaourca raquiraaanta, and which ara of linitad dura- 
tion. Thia haa anablad tha dacoapoaition of tha orar- 
all control algoritbn into a aat of auhalgoritloM, 
callad auhcontrcllera, which can oparata ainultanaoualy 
and aaynchronoualy. Coordinate tranafomationa of 
aanaor faadback data and actuator aat-pointa hava 
anablad tha further tinplification of tha aubcontrol- 
lara and hava raducad thair conflicting raaourca 
raquirananta. Tha nodulaa of tha dacoapoaad control 
ay^tan ara laplanan;.ad aa parallel procaaaaa with dia- 
joint aanory tpaca coanunlcatlng only by 1/0. Tha 
•ynchronlaation nachaniana for dynamic raaourca allo- 
cation among aubcontrollara and other aynchronlxatlon 
nachaniana ara alao diacuaaad In thia paper. Such a 
aoftwara organization ia auitabla for tha ganaral fom 
of nultiprocaaaing uaing coaputar natworka with dia- 
tributad atoraga. 

DITKODUCTIGM 

In recant yaara, tha anarging raquirananta for 
interactive conputar-aidad control of ayatana bacaaa 
prograaaivaly nor a evident. Tha nor a obvloua arena 
whara auch raquirananta appear ara apace, undaraaa 
exploration, nuclear energy-producing facilitiaa, and 
tha like, where tha work apace ia Inaccaaaibla or dan- 
garoua for nan to oparata in, yet where conplax taaka 
need to ba acconpliahad. Many of thaaa taaka ara gan- 
arally not suacaptibla to parfomanca by conplataly 
autonatad ayatana at tha praaant atata of autonatad 
daciaion-naking tachnology. It la therefore required 
to develop ayatana that project certain human capabil- 
Itiea, auch aa aanaing and handling, into tha work 
apace, thereby enabling rwota oparationa and procaaa 
control. 

Racota oparationa, fraquently alao called tele- 
oparatlooa, benefit tranendoualy from conputar inter- 
active control, where low-laval control and daciaion- 
naking functiona ara dona by tha conputar, while tha 
higher level deciaion fimctiona ara parfomad by tha 
human operator. The efficient allocation of functiona 
between nan and aachina ia remote oparationa haa bean 
tha central aubjact of intanaiva raaaarch activltiea at 
varioua inatitutiona during tha paat decade. For han- 
dling and aaaenbly in apace, tha development require- 
nenta and the atata ^f advanced technology for the 
control of remote nanipulatora hava been outlined in 
Uf 1. 


With thia backgromd ia nlad, we build on and 
axtand in thia paper acme of the devalopn en ta for ayn- 
chronoua control of ■aaipulatom ia lafa 2 and 3 and 
conatruct a nore ganaral framework for aaynchronoua 
control of oparationa and procaaaaa requiring the 
arheduling of nany aingle-actuator controllara by a 
** aup erviaory** i.ontrol conputar. However, it ahould 
be recognized here that a nanipulator ayatan can aerve 
aa a convenient, yet aufficiantly coaplez, framr of 
reference for the diacuaaion of nore general ayatana. 

Aaynchronoua control differa fron aynchronoua con- 
trol in that initiation of a aingle-actuator controller 
ioea not depend on cyclic interrupta of equal tine 
intarvala but, rather, on interrupta baaed on overall 
procaaa requiranenta. Thia anablaa tight acheduling 
with nlninal idle tine for each actuatsr. 

Tha inplanentation of aaynchronotia control ia baaed 
on the deconpoaition of the overall control activity 
into a number of control aubactivitiea which affect 
different actuatora or groupa of actuatorn. In tha 
interactive control procaaa, theae actuatora or groupa 
of actuatora ara active only for a limited duration. 

The control aubactivitiea ara generally performed aiaul- 
taneoualy but mutually tine-independently , i.a., 
aaynchronoualy • 

In thia paper, we firat diacuaa tha general phyai- 
cal atructura of tha interactive control ayatem. Then, 
we atudy the control functiona to derive their deconpo- 
aition into control aubfunctiona. And, finally, we 
propoaa a general aoftwara inplanentation of tha decon- 
poaad ayatan. It la ahown that all aubayatan conponenta 
can ba nada aa parallel, aaynchronoua procaaaaa which 
ara dia joint in the addreaa apace and which conmunic^ta 
between each other only by I/O. Thia approach builda on 
the worka of Hoara (laf 4) and Brinch-Hanaen (Raf 5), 
and ia oriented toward the uae of a general coaiputar 
anvironmenr, auch aa nultiproceaaor neeworka with dia- 
tributed atoraga. Accepting degradation of reliability 
and efficiency, the aana approach can be uaed in nulti- 
proceaaor natworka with cowon atoraga or In a aingle- 
proceaaor anvlronnant. 

GENEKAL SYSTEM DESCtIPTION 

Sinca our concam hare ia primarily with the con- 
trol conputar in Fig 1, we abatract tha ayatan aa ahown 
in Fig 2. Only unita coimunicating directly with the 
control computer ara ahown with tha correaponding 
incoming and outgoing infotmation flown. 

Coemand Device 

Tha conBMnd unit racaivaa coanmnda fron tha hixMn 
operator and tranalataa thaaa into a control computar- 
accaptable fom. Cnananda hava charactariatica of 
diacrata avanta and are generally iaaued through kay- 
boarda, voice, awitchaa, puahbuttona, and tha like. 

Hare, we conaidar ennnanrie in a ainplified fom aa a 
vector of Boolean variabiea c • (cx»C2» • • • tC^c) • 
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Fig. 1 FuDCtioBAl tubsyttwt and lafoiBatlon flow for 
iatoractlvo eontrol of oforational procoMoo 
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Fig* 2 CoBoral conflguracion of an 
iBtaractlva control ayacoB 

which rapraaanca 2*^^ poaalbla Inatructlona ragarding 
tha procaaa by which taaka ara to ba axacutad. In 
Fig 2 tha ayabola which rapraaant inputa and outputa 
of tha dawicaa hava wawa aign that indicataa '*raw** 
valuaa of data. Thia aign will ba ignorad in thia 
oaction. 

PiaplaT Pavica 

Canarally, tha dlaplay davica racaivaa, froa tha 
control coBputar, all InforBatioo raguirad by tha 
huaan oparator to aaka appropriata daciaiona for ovar* 
rida control. Thia nay include video data and audio 
aignala aa wall aa force and poaition data, all of 
which ara tranafornad into a huaan oparator-convaniant 
for* for diaplay and aanaing purpoaaa. 

In aiiplar varaiona of diaplay imita, one haa 
only a fixed nuabar of diacrata aaaaagaa, such aa 
warning and alara aignala. Thaaa Maaagaa can ba 
rapraaantad by tha character anring a ■ ^1020304 . . . 
which ia ganaratac by tha control cooputar. 


(Si •S2» ••••>«•) of thi na aimltanaoua aiipiala 
lad within a particular tlaa iacarval. 


Aitmtovt 

The actuatora ara a coUaetioo of N davicaa each 
conaiating of actuator unite, aarvo unite, and intamal 
atata aanaora. Tha intamal atata ia daacribad by tha 
vector of phyaical atata variablaa ^ • (^, •*>«)• 

lach actuator axacutaa aat-point valuaa or alnply *aat- 
pointa** racaivad froa tha control conputar aa a vector 
dp ■ ^SpxfSp2*****^^* ^ phyaical-atata variablaa 

are neaaured and aant back to the control conputar 
aa a aanplad vector dp - (Sy^tS«2f«.»yg). The 
actuator can axacuta any adnlaaibla aat*point ia a fin* 
ita tine due to ita own internal feedback aarvo control. 

Control Conputar 

Tha inputa of tha control c onputar are aanaor data 
S, annual controla N, connnada c, and feedback atata 
variablaa jj. Tha outputa ara naaaagea ^ and aat-pointa 

•f- 

Tha critical control outputa ara tha aat-pointa 
Thaaa outputa ara generated iteratively by conplai algor^ 
ithna baaed on continuoua and diacrata intamal variablaa 
and functiona. Ganarally, tha output valuaa jp do not 
only depend on tha preaant valuaa of tha inputa nantionad 
above, but alao on their peat valuaa, i.a., on tha con* 
plate hiatoty of tha control procaaa. To avoid handling 
of infinitely nany current and peat valuaa of tha input 
variablaa, a well*known practice ia to introduce atata 
variablaa and tranaition functiona. Thaaa atata varia* 
bl^a can ba grouped into a nunaric atata vector £ • (rj^, 
Bacauaa tha algorithna of tha control 
conputar alao deal with logic variablaa, it ia necaaaary 
to introduce a logic atata vector ^ • • • *^L) 

which ia a vector of Boolean variablaa. Generation of 
tha iterative aaquanca of the aat*pointa {0p(h)} 

(k • 0,1,2,3,...), where k rapraaanta tha iteration 
index, can now ba daacribad by a racurranca relation in 
tha following general fom 


i — • .1 *• 

. .,„(k) ,(k) ^(k) . (k) ,(V) ,(k)v 

i i i * £p • — * • ' • 

,(k41) . . .„(k) -(k) (k) , (k) ,(k) .(k), 

— s.'!i »i •£ *2, •£ ’• 

k-0,1.2,3 


Manual Controller 

Tha manual controller conaiata of ana or more 
hiaavi-factorad davicaa by which tha hunan oparawor 
ganarataa control aignala that affect tha axacution 
of the taaka. Tha hunan oparator acta aa a aupar^ia^'^ 
by overriding and/or inatnicting tha control conputar 
aa required. Tha generated aignala ara converted into 
an appropriata atandardiiad fom and aant to tha control 
conputar aa a vector H • . . . .Hoy) of tha valuaa 

of tha nn ainultaoaoua aignala ■aaplad^within a par- 
ticular tint interval. 

Sanaora 

Tha aanaora rapraaant a aat of aanaor davicaa 
which eaaaura^tha af facta of tha actuatora on their 
anvironnant and the deviation of a control procaaa from 
a praacribad atata. Tha aanaor data ara tha aaaantial 
alanant of all autonatUally aupportad, i.a., hunan 
operator* independent oparationa, that appear aa a part 
of tha overall interactive control procaaa. Tha aanaor 
aignala ara converted into a atandardiiad fom and ara 
aant to tha control conputar in fom of a vector 


-e 


Hare r^ and t ara initial conditiona of tha nunaric 
and logic at^a vector raapactivaly , while g, and 
h ara general vector functiona that rapraaant tha con- 
trol algorithna inplanantad on tha control conputar. 
Tha laat two functiona ara often called tranaition 
functiona. 

Tha aathanatical conatruct of Eq 1 rapraaanta 
a aaquantial nachioa, tha output of which can alao ba 
rapraaantad in tha oparator fom 


»P • i. £. ip. lo* io>* 


where ^ ia a vector oparator defined by Cq 1. In tha 
following, tha oparator ^ will ba called **control 
algorithn" or ainply "controller**. 

Tha ranaindar of tha paper will ba devoted prinarily 
to tha intamal atructura and functiona of the controller 
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Tim eoDCrollar li uMally iapltBMt«d m a bobo- 
lythic aofcifara packagt, l.a.t m a aiagla eoaglox 
aaguaatlal prograa. ABichar approach to acntetorlng 
tha eootrollar la an appropriata dacoapoaitloo into a 
aac of aubalgerlthaa or aubcoatrollara ** *iii* 

Thia aaeoad approach la dlacoaaad hara. Ca aparad to 
tha aoaolythlc apptx aeh, it appaara to of far laaa coa- 
plaxlty and bactar portability t flaalbllityt aad aala- 
talaablllty. Noraovar. thia approach of far a tha 
poaaibllity to oaa aora attract Ira procaaalng altam- 
atiraa* aocb aa aultiprocaaaiag. which could ba tha 
ooly aolution whan tlaa eonacralnta bacoaa critical. 


CoBtrol Ta«to Actloo. 

To docoapoaa tha eoatroUar, it ia aaaful to con- 
aidar tha atructura of tha ayataa taaka that hara to 
ba par formed by tha actuatora. Tha actiooa of tha 
actuatora hara rarylng dagraaa of coaplaxicy. For 
aanipulator ayataaa thay hara baan grouped la laf 2 
aad 3 Into chraa catagorlaat (1) a fisita aat of prla- 
Itlra actlooa (2) coapoalta actlona conalatlag of a 
•aquanca of prlaiclra actlona* aad (3) coaplas actlona 
which are nacworka of coapoalta actlona that era oaa- 
cut ad aaquantlally and/or In parallel, la aggregate 
of coaplas actlona nakaa up a ayataa tank. Haro wa 
•hall not ba concamad with tha hierarchy of actlona. 

Tha ayataa teak will alaply ba conaldarad a network 
of actlona which* In chaaaalraa* era dafiaad by apaclflc 
control algorlthaa and which hare apaclflc I/O require- 
nanta. 

Initiation and taralnatlon of an action during 
teak asacutlon dapanda on tha atata of tha ayataa 
which la defined by tha ralatlra '*>oaltloo of the actu- 
ator In Ita anrlronaent and by other eranta that 
deacrlba tha dagrae of teak coaplatlon. For thia 
purpoaa* tha evant-atatua ract.' * • ■ (ax*a2» • • ••no) • 
haa bear Introduced. It la a Boolean rector whoae 
coaponenta Indicate occurrancaa ot abaancaa of tha 
corraapondlng eranta auch aa **end effector la wall- 
allgnad with the plane*** or **the object la graapad*** 
etc. 

Tha Initiation of actlona la a chain-raactloo-llka 
proceaa* where one action Inltlaraa the following action. 
Theae action Inltlatlona propagata through the action 
natwork until tha ayatea tank la coaplatad. Tranaltlona 
from one action to another action are defined by tha act 

of Boolean functlona Tij(c*a), 1*J - 1,2 idj* 

of tha conaand vector £ a^^the erent-atatua vector £. 
Thaaa are called tranaltlon condltlona. If an action* 

•ay 1* la In asacutlon* and If tha tranaltlon condition 
Tij bacoaaa true* then tha action J will laaadlatcly 
ba Initiated. Tha termination of an action la defined 
by a alallar Boolean function txi(c,a>,l • 1*2*... *n. 

If tha action 1 la In asecutlon and *lf bacoaaa 
falsa, this action will ba tarmlnatad laMdlataly. If 
Tij • true Inpllaa in • falsa* then tha action 1 will 
ba tarmlnatad alnultanaoualy with the Initiation of 
tha action j. Otharvlsa tha two actlona will remain 
In parallel asacutlon. 


Subc patrol lar a 

Tha aubcontrollara are daalgnad ao that each one 
correaponda to a alngla action. Tharafora* thay will 
generally ba asecuted aa asynchronous* aaquantlal pro- 
graas which can ba Independently coded* tasted* and 
Integrated with tha ayataa by aatlafylng certain aub- 
cootrollar Interface raqulraaanta to ba dlacuasad 
later. 


where tha aMsrle mA logic state rectors ia and U 
and thalr raapactlra Initial eondltlcms Ijo are 

derived by appropriata daco^aaltlon of M atata rectors 
£ and £. The rarlablaa ipj • Nj • ij t and |jj are sub- 
rectors of Ip* H* £ and Ip raapectiaely. 

Tha suhrMtors hare tha aans dimmisiona aa the 
correspond Ing rector in which tha co^onanta of no 
latareat for tha action imdar emacutlon are Ignored la 
the coeputatlona. TWs a u bm e tore are mtuaJly disjoint 
If thay do net hare noml^Mrabla covenants la eomnon. 
Therefore* tha argMMta of |j Includa only thoae com- 
ponents required for tha confutation of tha sat-poiata 
associated with tha action j* In othar words* for an 
action j tha subcontrollar |j gamarally does not need 
to coaputa tha aat-palMts for ail actuators and dees 
not need to use all seaaors* tha manual controls* 
sad the actuator feedback. 

It should ba noted* that some actions are fully 
autoaatlc and do not depend on the tnmm operator. 

Thaaa actlona do net call for nanus 1 Inputs* l.a.* 
tha subractor is empty. 

As seen* Iq 3 does not Includa eJeneats of the 
coMond coda rector e ao aa Input srgumsat. As shown 
aboru* this rector la an srgnnmit of tha tranaltlon 
conditions Txj(f.*s) which coordinate tha actirltlas 
of tha sabcontroller. 

Generally* two different subcontrollers 6^ and 
^ can use tha asne coap o nenta of the Input rectora 
IT and £* and can geaarsta tha aat-polats of the aaaa 
actuators. If ouch subcoatrollers work slaultana- 
ously* they Interfere with each othar. raarefore* 
nubcootroliers conpetlng for the asne actuatora should 
be coordinated ao that thalr outputs |pi and |pj are 
nutually disjoint subrectora. Kowerer* this can con- 
siderably reduce the n\aber of actions that can be per- 
fomad simultaneously. In order to solve this problen, 
we Introduce s coordinate transfometlon 


fl •!.(•)• 


(*) 


which nlnlmliea the coordinate Interactions for the 
aajorlty of actlooa. Mote that this trans format loo 
auat hare an Inverse 


1 ■ 


(5) 


The transformed eet-polnts £p • (Qpl »0p2* • * 
represent the new space colled ''controller apace**. 

The coordinate trans format loo (Eq 4) also slap Ilf las 
the aubcont rollers * because the controller space Is aora 
convenient for constructing the algorltbns by which the 
actions are laplenented than the previously used actuator 
•pace. For aaaaple* the controller apace la the case of 
smulpulator control ay stems* Bef 2* can be the hand 
coordinate ayctea, while the actuator apaca Is the ^*>lnt 
•pace of the aanipulator. The change of only one i uor- 
dlnate of the hand coordinate ayatea, say the front 
distance o^ the end-effector from the object being nan- 
Ipulated* can cause alaultaneous adjustments of several 
joint coordlnatea. The coordinate transformation 
(Eq 4) In tha case of nonredundant asnlpulstora Is 
defined by the nenlpulstor geometry equations. 

In order to make further simplifications of the 
•ubcontroUera* a similar transformation for sensors 
data Is introduced 


I • 


( 6 ) 


A subcontrollar associated with an action j, 
generates the aet-polnta for those actuators which used 
to be active for the aaecutlor the action. The corre- 
sponding operator aquation csi be representad as follows 


where T - (Yx*T2». • • .Tgr) Is the rector of the trans- 
formed*^ sensor data. Exaaiples are outputs of strain 
gauges placed In the wrist of tho reelpulator arm, 
which are transformed Into orthogonal forces and 


i 
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tOT^yi 4yfiiM4 la tha haad eoordiaata ayatM. Traaa- 
foraacloo (Iq I) doaa aot a acaaaa r lly aaad to hava aa 
lanraraa, l*a«, VT ^ aa* TIm aaaual eoatrola N do aot 
haaa to traaafofaad« b aoa a aa tha aaayal eoatrollar 
eaa daalpiad aa tl^t It diraatly gaaarataa aitaala 
that Qorraapoad to tha ooatroUar ayaca* laeaaaa of 
tha Itaratlaa natura of tha auhcoatroUofa aad aoaa 
foraal coaaaniaacaa* tha oatpata of tha auheoatrollara 
will ha rmaaaatad hr tha laeraaoBta of tha aat-yoiata 

Tha aaa oyaratora of tha 
aahcoatrollara can aov ba vrittaa la tha fon 

lha oyaratora ^ ara axyaetad to ha alaylar chaa 
tha oyaratora of Iq 3 ha c a u aa of traaafotBatloaa 
(Bq 4 Mid Iq 4)* Thaaa traaafoiMtlooo can ha ean- 
trally aaacutad fnr all ai^hcoatroUara tad caa ba par- 
foraad vlthia aaothar faaetioaal block to ha dlacoaaad 
latar. 

Aa aaao* cha aub^ractor la oalttad la Eq 7. 

It la aot aaaaotial to tha eoatrollara daftaad la tha 
eoatrollar ayaea. Aa vlll ha ahowa latar* 9j la aaad 
vlthia tha actuator block (Iqa 13 thru 17). la fact* 
cha could ha §aaarally ratalaad aa aa argiaaat of 
^* hacauaa aoaa algorlthaa of ouhcoatrollara might ha 
baaad oa tha laforaatloa of tha actuator atataa* hut it 
vill ha omittad hara for tha aaka of tha alaylicity. 

Channala 

A trlpla of layut aad output auhvactora (M**Y** 

Dj) ia aaaociatad with ovary aiibcaotrollar u the 

folloviag dlacuaaioaa thla triple will ha callad 
'*chaimal'*. channala (3x*Xit5x^ 

ara mutually nonin tar farina If: 

1) Inyuta Mx do not affact tha outputa U 2 * and 
tha Inyuta M 2 do not affact tha oucyutu 

2) Qiangaa of tha anvironmant dua to tha axacu- 

tion of Ux cauaa changaa of Tx* and 

changaa dua to tha ajMcutlon of U 2 will not 
cauaa changaa of Tx* 

It follova that tha nacaaaary condition that two chan- 
nala ara mutually nonintarfarlng la that tha corraa- 
poodlng output vac tore Ux and U 2 ara mutually disjoint 
•uhvactora. 

If tvo auhcontrollara work on mutually Intarfaring 
chanaaia, chay can Intarfara with aacb other* iacauaa 
they uork aaynchronoualy* l.a.* time Indapandancly * It 
la practically Impoaaibla to dacarmina and to control 
tha Intarfaranca baevaan cha auhcontrollara. Thalr 
•Igorlthma muat tharafora ba daaignad laidar tha aaaiay- 
tlon that they ara complataly Indayandant of each 
other* Conaaquantly* only for auhcontrollara working 
•imultanaoualy on mutually noninc erf aria g channala can 
stability ha guarancaad* l.a.* tha convarganca of cha 
ccrraapondlng algorithms. Tharafora* tha ay scam must 
provide tha channal managamanc that will anahla tha 
coordinated allocation and deallocation of mutually 
Interfering channals to alaultanaoualy active 
auhcontrollara. 

Iacauaa aoma actions can ba more urgent than 
othiirs* aa It ia with actions Invoked In asMrgancy 
situations* prlorltlaa should ha aaalgnad to cha 
acta.ona. Thaaa priorities will rasolva conflicts whan 
two or mors suheoncrollars ara competing for mutually 
interfering channals. Tha matter of priority assign- 
manes vill ba diactiasad later. 

Having dona tha dacomyoaltlon of tha controller* 
tha fiAccional block diagram of tha control computer 
caa ba rapraaantad as shown la Fig 3. Thera ara six 
functional blocks* five of them having assigned I/O 
functions and ccordinata traaa format Iona and one with 



Fig. 3 Oiatribmtiom of fumcclona of tha 
control comyutar 


control functions* lha control block eonaiata of n 
auhcontrollara and tha nacaaaary intarfacaa which 
enable data coamunication ha t waan si^omtrollars 
iatamally and hatwaan auhcontrollara and pariyharal 
blocks axtamally. Tha iatarfacaa also hava tha 
raayonaibility of channal aaaagamant. Thaaa fiaic- 
tlonal hlocka and iatarfacaa ara diacuaaad In tha 
next two aactlona. 

MAUI S7STIM COMFOHPrrS 

Tha functional blocks of Pig 3 will ba Impla- 
aantad as parallel and asynchronous procaaaaa. In 
accordance with tha ayatam structure diacuaaad so far* 
wa define five peripheral procaaaaa and n procaaaaa 
which support subconc rollers Tha latter ara 

callad action procaaaaa* Than^rlpharal procaaaaa 
support tha corraapondlng pariyharal functional blocks 
aad ara callad conaand procaaa* display process* annual 
controller* aaaaor procaaa* and actuator process. Thaaa 
procassaa will now ba diacuaaad in greater detail. Aa 
vill ba seen latar* thaa^i procaaaaa ara not tha only 
systam coaponanta. There, ara additional conpenants 
which appear aa a result of tha decomposition* aad 
which support tha communication hatwaan tha procaaaaa. 
Thaaa components will ba considered in tha next section. 

A more detailed and formal description of tha 
systam Is given in tha Appendix* using tha notation 
developed by Dljkstra (Haf 4) and Hoara (laf 4). Tha 
reader not familiar with thla notation can skip tha 
Appendix without loss of tha Ideas praaantad hare. 

rnini frnrm 

Asaumlng* for tha aaka of simplicity* that tha 
command device la an array of nc on/off swlcchas* than 
tha command procaaa cyclically lutarrogataa thaaa 
switches* whose poaltlona ara defined by tha raw command 
vector F • * • * »^nc^ * It transforms this vector In- 

to tha command vector £ and checks If there la any change 
In tha coMands. Checking la dona by comparing £ with 
its value £^ 2 ^ from tha previous cyelr. If £ d c^xd* 
than tha new value £ will he sent to the action pro- 
cesses for further conputationa. 

Display Process 

Ha asaums again tha simpler case that tha dicplay 
device la Intended for discrete massages such as warn- 
ing and alarm signals. In thla case* the display pro- 
caaa racalvas* from tha control block* a uaasaga code 

in tha form of a Boolean vector a • (ax *02 

whenever It has a new value* It produces* than* tha 
character string o • 0 x 020 ^ 0 ^..** and displays it to tha 
human operator. Tha cnaractar string Is generated by 
the program C OW TgT. which for each state of a responds 
by a fixed* predefined verbal mat saga. Tha dacaila of 
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i^ich 4 Ucumw4 U Om ^rtirlova MttlM. flMlly 
it. teflSM tiM tomlmm vector tmmif ■ OrM^yxtrMiy}* 
thm eo^OMitt of oliieh 4oflM tho Mt- 
:»olBta of tiM octootor yroeoM tiMC Imoo boon alroo4y 
otcooitoi, l.o.» 


Thio yrocoM la yorfoimlat tivoo oforatlooo. Fint 

S t qrclicAlly rood# tbo row vaXuo of tho ■■ wi l l eootroXo 
__ froa tho A/D eoBoortora ohich oro ot t o ch od to tho cca- 
?rollor dovleoo. Socond. It colihrotoo roif ooluoo !■ 
occordooco with tho rolotioa 


ohoro tho ^etor fiactlMi |n royrooooto oil oyorotlono 
Bocoooory to tranofoni H into Ptoadordlood oIimIa* 
thuo proTldias • WOod latorfoco ho tw o on tho hoMO oyor 
otor Ml tho control eovyutor (llnoorlootion, conyon- 
ootloB for uoMBtod offocto duo to nonltnoor ehnmctor 
iotlco of tho controUor dowlcoo, introduction of 
ooturotion and dond-hond» otc.)« Third* tho yrocooo 
chocho if tho now unimol controlo M oro difforont fron 
thooo lonorotod in tho yrowiouo tyclo IJqXA ^ ehochini 
tho inoquolity 


idiom in tho mctor of tho oetunl oot-yointo* and 
^ io 0 giwm toloronco mctor. Tho mctor roo^ in 
oout to tho action yrocooooo in ooch yrocooo cyclo. 

In tho outynt ynrt tho octuntor yrocooo rootimo 
tho total incmnont (if any) of tho oot-yeinta (thia 
mctor ia eooyoood of tho ynrtinl incroaanta Ui gonor- 
atod hy tho difforont oobcontrollora uorking alanltan- 
ooualy) and foma tho now aot-yointa tiainc tho enrront 
yooitiona ^ 


ia a glwoB toloraaco mctor. If tq 9 dooo 
tho now mluo H will bo oont to tho action 


Tho oot-yointa in controllor oyoca oro now tranafomnd 
into tho actuator ayoco by tho laworoo noyying 


ard convortod into tho raw output woluo 


Tho mctor fiactlon oyoratoa aiadlar to but can 
bo uoro conyllcatod duo to tho coayloxlty of aono aonaor 
ayatono. For oxaoylo* proxiaity aonaora with fibor 
optica Buat bo calibratod for onrlroanontal conditiona* 
auch aa color and raxturo of tho nonlpulatod objocta* 
otc. Thorafora* tho calibration yrocoaa nay nood a 
Ornanic doflnltlon of tho callbratloo paranotora. Tho 
fubctloB nay alao includo a filtorlng capability to 
roduca nolao. 

Aftar calibration, the procooa choclrn if tha now 
mluo % ia changod with roayoct to tho mluo ^xd 
tho prorloua cycle by tho inoquallty 


Thia ia than wrltton on tho D/A conmrtora of tho 
phyalcal actuator. 

Csocution of tho input and output port of tbo 
actuator yrocoaa ia not ooonntlally aynchronoua and 
altomatim* but tho aoftwnro inplonantation nuo** 
onmiro faimoao to both porta. Alao, to amid dalay 
affocta ia tho control looyo aonaor -control lor- 
actuator. tho cyclo poriod of tho actuator procooa 
Buat bo long enough in conparlaon with tho cyclo poriod 
of tho aonaor procooa. 


Tho action procoaooa yorfom tha actiona. l.o., 
thoy oxocuta tha aubcontrollara ij (j ■ I*2*...»n) and 
prowldo tho nocaaaary adnlnlatratlon. Thoy all ham 
a ainilar atructuro, with tha oxcaytiou that thoy 
aaiploy a difforont algcrlthna For tha aubcontrollara. 
rhoroforo* tha action procoaaao will bo rayroaontad aa 
an array of procoaaoa with tho oubacrlptad naaaa Ad). 

A(2) A(n). Bwory ar«;loo procooa can ham a bloclrod 

and an actlm atato. In tho blocked atato, tho pro- 
coaa ia waiting to ba actimtod by aono othar action 
procaaa. In tha actlwa atata. it cyclically chocka 
tho ralomnt tranaltion cooditlona and oxocutaa tho 
aubcontrollor. if it haa tho right for furthor 
axlatonca. 

Chocking tho tranaltion cooditlona lor aa action 
procooa. aay A(j). oaana ovalutiag and chocking tha 

ioolaan functiona T4 x(£.o) (i • 1.2 i d j). 

idilch wara dlacuaaoa la tho prowloua aoctlon. If aona 
of thoao axyroaaiona turn out to bo true, tho p'oeoaa 
A(J) will notify tho ayaton which will then aond rwaka 
algnala to tho corroayocdlng action procoaaoa that 
ouat bo actimtod. If Y4j(c,o) boconaa falao. tha 
yrocoaa A(J) will iaB»dlatoT/""yuc itaolf la the blcckcd 
atata. Othorvlao, it will axocuta tho aubcontrcllor , 
1.0. . it will invoke tho procoduro with tha naoa fUl- 
COimOLLZIj which ia ayoclfic to tho corraayondlag 


where C5 ia a given tolaranca vector. If Eq 11 ia 
true, tho procaaa will read now raw aonaor data fron tho 
A/D coov^^nom. OtSorwiao. it will flrat porfom tho 
coordlnata tranafomotion 


which haa boon dlacuaaod in tho provloua aoction, and 
it will than road now raw aonaor data. Tha tranafomad 
aonaor data X ^ ^ action procoaaoa. 


Actuator Frocoaa 

Thia procaaa haa aa input and output part. Tha 
input part cyclically roada raw feedback valuoo of tho 
actuator atata mriabloo and porfoma tha data 
raiibratlon 


It than tranafoma tho mctor fy fron tho actuator to 
tho controllor apace by a napping 
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•cclM. WmmtU. tmd logic otcco ooctoro n ood ^ 
of tbo cubcoocroilor on loccl vorlokioo ^ tlio 
grocooo A(J). 

Cowormiat tU octloii grocoooM end tiMir ooA co a - 
crollors* tlwro «ro coo lagortaac icoooo Booclo”iod la 
tbo proolooo ooccloa. Tbooo cro ehMMol cllocotioo and 
eoatroUar aya chr oalaatloa > 

Wkea too oi^OBcrollora era eoapoclag for thm •amm 
chaaaal or for diaaaola that era aataally iatarfarlag* 
oaly oao oaheoatroUar cm ha aaalfaad to tha chaaaal. 
The other av^t aait oatil the firat unHaataa aad 
rolaaaao Ita chaaaol. Xaaadiataly after halag acU- 
iratad* m ectioa procaaa ragaaata froa tha ayataa tha 
chinaal raguirad hy ita aabcootrollar aad thaa aaita. 

Aa aooa aa tha raguaatad chaaaal ha c oa a a arailahla, tha 
actloa procaaa will ha aotifiad aad activated. It will 
thaa itarativaly aaacuta tha aaheoatroUar laitil tha 
corraagoadiag actiaa haa haan coaglatad. Whoa tha aah- 
eoatroUar taraiaataa« tha actioa yrocaaa aotifiaa tha 
ayataa that the chaaaal cu^ bo ralaaaad . Tha chaaaal 
caa aov ha naad hy aaothar aaitla^ actioa yrocaaa. 

It haa baaa autad that aoM actioaa aa> ha aor^ 
orgeat thaa othara. For iaatfficaai jaargaacy actiona 
aoat ha yarfoiaod ianadiataly to pravaat ayataa fail* 
ura» iacorract parforaaaca* or colllaioa. la thia 
caaa» tha m mrgmc j action ahould ba ahla to get tha 
chaaaal avan if it haa baijn allocatad to oaothar 
action procaaa. In ordar to adalnlatar thia, tha pri- 
ority vector £^(pitP2f*tPn) ^ iatrodacad, whom tha 
prloritlaa pj (J*l,2,...n) era aaaociatad aith tha 
procaaaaa A(J). If tha procaaa A(k) haa a higher 
priority tl^aa A(j) ahich ia la prograaa, (p^ ^ Pj)» 
than A(k) aill ba crantad laBadiata accaaa to tha 
channal Hhila A(j) aill ba raapo«’anly blocked. 

Bara, va only conaidar atatlc prlorltiaa. Horn gan- 
aral capacilitiaa of the ayctm oould ba with dynaalc 
prlorltiaa which can ba changed during taa^ axacrtlon 
by tha Ynmmn op^ %cor or by toaic algorithi. 

Tha aynchre liaatlon of a a^«bcontroUar eonaiata 
of tha folloalng. Tha osacutlon of tha tubcontroXlar 
^ wltain tha action piocaaa A(j) ia dona ia aa itaia- 
ti/a aonaar, l.a.. tha aubcontrollor ganarataa tha 
aaguanca of cha aat-polnt iaeraaMnta {Uj(k)} 

(k«0,:,2,3,. . .). Ha.*a, any lacr^Mat ;^(k^ cannot 
ba tent to tha actuator procaaa bafora tM prartoua 
Incranant haa bean uMcutad. In other uorda, 

whan vanta to ganara‘;a and «and c m it 

Buat wait until tha aetjator prccaaa firat asacutaa 

Ana than noflfiaa A(J^ by tha algnal **go** that 
It can proceed. Only than can A(j) aak for, and 
racalvt, rha new valuaa M(a) and T^k) to coaputa the 
DOW raJua JJ](k) and than"*sanc it to the actuator pro- 
caaa. Thia aynchronlaatlon nachanion. and tha pravl- 
oualy dijcuaaad channal allocation aachanlma, are 
axtamal to tha action procaaaaa. They will ba dU- 
cuaaad in uha next aactlon In aora detail. 

crmFACCs and ststdi ikteciatiob 

In tha pravloua aactlon, tea aaln coapenanta ware 
dlccuaaad. In ordar •‘o provide a alapla and raliab.a 
^laaan^atlcB, tha procaaaaa defined ^heuid not cc.#- 
mmlcara directly with ^h ether, but the: r hoy Id 
coanunlcata through aona Intar^acaa. F tr thia pjr^ 
poaa wa Introduce two addltioril aofrwart c<>npooantc 
whlc> are a atawua aonltor and a channal oo^itor. Thaaa 
nenitora can ba built ac paaalwa rnaponantw w^th a data 
atructura /H aaaociatad oparatlcnr (nooitor pro- 
caduraa) aa waa pToacaad In tha pravloua wotk (laf 3) 
which waa baaed on aonltor cnoc^pt by Irlnch-Baaaan 
(Faf 7) nd Boara (laf I). Tha IxMt^cioe of thia 
approach ia tUrt tha pmcraariu accaaaing Ltm aonltor 
nuat ahara coonjn jmorr apace. Bowavar, a aora 
attractive poaaibility ^a if wha procaaaaa are aada t/% 
ba Btttually aiajoint la aaBory apaca. Than, ihb ayataa 


can ba lapIiMacad on a c oape t er aatwork with diatrih- 
atad etorage. Solutiona to thia problon have bean alao 
propoaad hy Boare (Baf 4) and by Brineh-Banaaa (Baf S>, 
wham tha noaitom axe inpleneatad aa proceaaaa which 
noonuaicata with ether praeeaaea by I/O. Thaae idaaa 
are applied in the following neragmphe. 

itatua Hottitor 

The Itatua noaitor eepporta (a) data coMunica- 
tion hatwaen all parlpheral procaaa a a (except the 
actuator proar'a) end tha action proceaaaa A(j) 

<b> evalMtina of the event-etatue vector 
a and the neeaag i c ode vector mi (c) activation of 
Tha action procaaaae. 

Tbo data rneMMicatita ij organ 1 cod ao tint tha 
etatua nonitor ewandnae other pcocajaaa to datemina 
if they ara ready to racaiva a. to aand data. T.v 
eonanid pracavia, nniial eoatroUar, aad aeaeor promaa 
eand dowo to the etatua nonitot whe n ever they am maiy, 
aad the etatua noalior will mcalvo theoa Mata end •r'U 
copy then into ita local vaciablao. Whan tha action 
procoe )aa naad thaaa data, they aak for chan by aaading 
tha appropriate aignalA to the atatua nonitor (the 
algnal **atatuo ’ for c and a, mA tha aigLol **xnpu^'* Tor 
N aad J)* Tha etatua aoniTor wUi ^.nva di ttaly mapond 
^ aondlug tha cormapondiag data to the eorreapeading 
prucaaa. The actioa procacaaj neat aak toe tha data, 
aacauaa the atatua arnitor aarva aom than oac procaaa, 
and it nuat ioapona to their fragtaont roguaato 
thai afore ahou. r' jot ba dolayad for longor porloJo. 

Uhon tho Itatua noaitor : icolvw eanaor data T it 
will alvayj invoke the two proeadtioa tVCfTTFSATF aad 
KBSSAdOODi. Tha fiiat procedure deflnee end updecee 
the now v.ilun of the event-etarur vnetor e« uelng chn 
neet mcent ecoeor iate end poaalbly their '%et 
valuaa. Peviag T, tad thn procndu.4 HESSAdCOOF 
daflaae aad updaTaa^the aaoaage code a. If them ie a 
change of thia vnet^ with raapact to the old vaiuo, 
i.t. , if £ p loU* •Tatna noaitor will aand tha nou 
value a to tho Replay procooa. 

Thi no "haul iF of activating aa cctirvn jrocaaa, 
aay A(k), by aootaa.* actiou procaaa, aay A(j), ia por- 
fomad aa follows. If A(j) fiaft that the traaaltlon 
condxtlon T 4 |^(£, g) ia true. It wlU aeno tha Inragor 
k to tho cijtua aonltor. This latogor Idontlfloa tbo 
rction procart that nuat bo actlvatod. Aa soon aa 
statue an Itor mcalvaa k. It wiU eonJ the "wfaka'* sig- 
nal to tho i^rocma A(k). This procooa. If in Jio blocked 
ante, l/“ waiting for tb^ •'awake** signal fr-n the ...atue 
nen.tor and wiU change ita atata to tbo activo orvi. 
r I **cwake** rignal aaa affect If tLa procaaa Av't) ia 
almni*y in the active etacc. 

C hiToal honltc^ 

The chewnel nonitor aup:«arta throo functlonet 
(a) drea cjoBunicar ton botvoeu action procnotva and 
actuator procooa, (h) channel allo' tion, end (cf ches* 
nel eynchronixeclon. 

The flmt fiBctlon la baaed on the aene piln- 
clplea aa for tho atatua nonitor. Tho nonitor rmolvas 
froa tha oct^or. pricooa tho Sooloan vocto«* roady and 
tends thn vnetor of eot-polnt Incrtnente j/tet tc tho 
actuator procooa. Tha vector £tot U dlacueao^ below. 

Tho cheanel allocation, l.o., oorvlclng too 
rosuoat/raloaao roqulranonta of tho action ymcaatM 
la diwia by jwr p.ocaduraat ISQUCST a^d BZLCAIF. Tha 
fl^at procedure la invoked whenever a "regueat** alpial 
la sent f*on tha action proraaa. Tha pxranetar 

of tha procodura la «n Intagar j t7it Ida riflaa tho 
rriuaoring procooo. Tho procoduroe haa eccaaa tv the 
Bxa Boolaea natrlx a • [n^j] (1^1,2,... ,Bt 
j«l,2,...,n) which laflnna tha channel fTeuimneacs 
ftr »11 ec** ^ pweeeee. Fir exMpie, nij • true 
Indlcetee that the 1-th input of tho act^tor procaaa 
la rjgwimd by tho awhcoutrollor Aj. The j-th eolian 
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- (•! 4.^1 HO ^ Mtnx ■ tlM 

riitinl ra^ttirMBC of thm oo^eoocrollor Aj. 

TIm rfoeo4«ro oloo 1 m occom to tlM priority roctor 2- 
Tho BOtrls a Md tiM roctor 2 coaotooto local to 
tlM cliiMilVoaitor vtiicb arc dofiaod at thm tlao of 
oyotoa laatallatloa. Whoa a rotaaat la ao4o hy thm 
A(j), tho procodoro UQUIS7 flrat ehocka If tho trao 
pattorm of tho aj aotchoo tho portom of tho ehooaola 
alroaiy oiloeotM. If oot» tho choaaol vlll ho laao- 
dlatoly graatod to tho A(j) oad tho alp^l '*allocatod** 
irlll ho ooat to it. If yoo, tho procoduro will pot Cho 
pToeooo iados j into tho ooltlat foouo, which la alao a 
data atractara local to tho choaaol aoaltor. If tho 
roqoootlaf procooa haa a priority grootor thoa tho pro- 
cooaoo holdlag tho choaaol, thoo tho aoaltor procoduro 
will porfora choaaol roollocatloa. 

Tho hrecodoro IILIASX la lowohod whoaowor a 
**roloooo** algaal la ooat to tho rhaanol aoaltor. Tho 
input paraaotor of thla procoduro la alao tho latogor j, 
which Idaatlflaa cho procooa roloaolag tho choaaol. 
r^o procoduro oraaino tho waiting quoua for a procooa 
waiting for tho choaaol. If auch procooa la fouad (prl- 
orlrloo aro cakya Into accoioit), aa "allocatod'* algaal 
.nil bo aoat to It, oad Ita Idoaclflcrtloa ladoa will 
ho rourwod froa tho waiting quouo. 

Tba third functloo of tho aoaltor, tkio choaaol 
aynchrc^laatioo, la dooo hy tho procoduro GnCKOUIOIEL. 
ChxM procoduro la Inwokod whoaawor tho choanal aoaltor 
rocalvoo tba woccci raady froa tho actuator procooa. 

It will thoo choc! for cnaaaola which aro atlU huoy, 
or which art r^ady to acc^t now aot-polata, Tho chock 
la aoda hy coaporing tho woctor ready with all coliaana 
cf tho aatria a, wulch corrorpooJ to tho actlwa chan- 
nola. Activa channola which aaot tho truo pactam of 
tho vector reedy ^ro riady to accept oar data. Tboae 
channola ara dafiaou hy the index aot Jg * (jlvi'l 
ready ^ Oj and jeJa) whora la »:ho iodox aot ^ all 
active chanaola. Tba procod««ra CHECKCHAWtEL will aond 
a **fo** algaal to all procaaaoa A(i), jeJg* 

Wo now turn hark to the vector L Tho channel 
ocoitor playa tho rolv. of a huff or collocta cllffor- 

ojt portlona of U which art aonr aa ouhvoctcra gj 
f.on tho dif foro^^^rocoaaca a(j). The colloctiot. of 
thoao portii*na li doiined by: 


U 


tot 


£ s’.i • 


( 20 ) 


whor^ y* are **clo'^r»d** 
dof Inad^by 


cor a loan of tho auhvoctor*. 
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Syatfu D ata Tlow Diagran 

lEo ayaton conpononna rjn can ho connected into a 
ayaton ahown hy .ho data flow ulagran la Tig 4> Tho 
procaaaoa are roproc anted hy circle* and tiai y«-riphoral 
davl.'oa *jf aquarod he aa. Tho data and CM aignala aro 
ropriaontof' ar^ owe indicating tholr aourcoa and 
their deot^ not Iowa. To indicata tUo difforonco borwere 
data and jipiala^ tho la:tor ara ahown ty daahed arrowo 

Aj aeon, tho diefkan in 7^g 4 hoa tho fore of t 
double Otar. Vth star* hm. too mnltor procaaaoa in 
^ho cortex, nia a-tloo ar^ceoata conatituta the c o oron 
hranchea of tha atara. Th i frew htaar.hea holon* to 
pariphoral proceaioa which have acceaa to tho axtomal 
pOyaical devicjo. 
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Fig. 4 Data flow diagran of tho 
intogretad ayat« 


OOKLDSIOM 

Tha propoaed aoftworo atructuro of tho intoractivo 
:ontrol ayotan la highly nodulariaod and la oultahlo to 
ho inplonontod on a eonputor network with diatrihutad 
atorage. The nodular loot Ion lo dona in accovdanco with 
tho varied I/O and control functiona of tho ayaton. 
iecaueo of thla, tho ayaton can he inplonooted ixedually 
and with tho poaalhillty of auhaequont roflnoBMct and 
inprovonent. All noduloa can ho doaigned, coded, and 
taated aoparatoly and alnoat indapondantly of each 
other. Tho aecond character la tic ia attractive hecauao 
of tho tine conatrainta and tha ayatan rollahllity 
which hecono nora and nora inportant factora of auch 
ayatano. Alao, tho trend of eonputor technology juati- 
fiwa auch an orientation. 

Tho intention Loro la not to give a doflnlto and 
dotallod iaplanontatlon of tho aoftwaro, but to Iden- 
tify certain atructural propertieo of tho intoractiva 
control ayatano and to given hlnta for tho application 
of nodom concopta of roal-tlno progroHilng propcaod by 
Dljkatra, Hoaro, Brlnch-Honaon and othora. 

For tha aaka of ainpliclty wo ignored. In tnla 
paper: (1) tha paranotorliatlon of tho coMonda c and 
noacegoa a, (2) tho dependency of tho aubc out rol lore 
(Eq 7 ) and aenaor trana format Iona (iq 12) on tho 
actuator-atato feedback jry, and (3) the dynanlc 
aaaigiaoot of channel prlorltloa to tho aubcontrolora. 
Thoao Jaauea will bo conaldorod in tho further work 
and through practical realiaationa. 
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Tho folloinag otatooonco ohculd not bo conoldorod 
oo a progroo, but oo a pooudocolo uooo to loocrlbo 
algorltboo. Thorofort, tho oothooo^'LCol oyobolo for 
voriobloi and fimctiono aro ujoo rathor choo ayobolic 
fioooa. Tho dotollod asploootioo for tho oyotia of tho 
aotatlofi lo giron In lof 4. Ir thla aoctlot. tho fol- 
loirlng oonotandord data typoo aro uaodt '’Icoloaa' 
voctor**, '*roaIuoctor** and **chartctorotrlng**. Thoao 
aro arraya of corrooponding atandard data typoo, uhort 
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ABSTRACT 


The CURV Arm Control System (CACS) is a computer aided control 
system for Interactive computer-aided control of the six-degree-of-freedom 
manipulator of the JFL teleoperator laboratory. The manipulator is 
equipped with proximity and force-torque sensors. In order to perform 
complex tasks like tracking, capturing and stopping of slowly moving 
heavy objects, the human operator commands are supported by automatic 
control algorithms based on sensory feedback data. The general objective 
of this development project is to evaluate the performance benefits of 
sensor-referenced and computer-aided control of manipulators in a complex 
environment. This progress report represents the first phase of the 
CACS software development, and gives the basic features of the control 
algorithms and their software implementation. The control structure 
development is based on three concepts: incremental motion synthesis, 

basic control routines, and parallelism of algorithms. Incremental 
motion synthesis consists of generating a series of motion Increments 
instead of generating endpoint values. This enables a unified handling 
of position and rate control modes of the manipulator, and uses simpler 
coordinate transformations based on linearization. The basic control 
routines represent the set of elementary algorithms for generating 
different kinds of motion increments shared by all algorithms of higher 
levels of the control hierarch. The parallelism of algorithms is a 
natural consequence of considering manipulator activities as integral 
components of complex manipulator tasks. Because of the relatively 
complex control structure and its Inherent parallelism, special attention 
has been paid to its software Implementation. Therefore, modem concepts 
of monitors and concurrent processes are applied In this work. 
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SECTION I 


FUNCTIONAL DESCRIPTION 


A. INTRODUCTION 

The CURV Arm Control System (CAfS) is a computer aided system 
which provides interactive human operator and computer control of the 
manipulator. The system generally consists of two parts: hardware com'- 

ponents and software. The hardware components include computer hard- 
ware • the JPL/CURV arm equipped with sensor systems and an operator 
control console designed as a universal control panel. These compo- 
nents are part of the JPL Teleoperator Laboratory and are described in 
Refs. 1-3. 

The CACS software is a new system compcneac Intended to support 
a class of real-time manipulator control actlvltloa such as tracking, 
grasping, and stopping of slowly moving objects. 


B. OBJECTIVES 

The objectives which are pursued in this project are given in the 
proposal "Develop Experimental Modeling and Evaluation of Sensor-Aided 
Manipulator Control" which was submitted to JPL in November 1978. The 
objectives will be briefly reviewed here. 

The general objective of the project is to demonstrate and deter- 
mine experimentally the Impact of sensor and computer aided manipulator 
control on overall task performance. The experimental nature of this 
general objective implies two major points: 

(1) The sensor/computer aids are tools in the hand of a human 
operator, and consequently, a major concern is to provide a 
proper interface between sensor/computer tools and the 
human operator . 

(2) The sensor/computer aids are real-time tools, and conse- 
quently, their performance properties should match the 
versatility of a real-time control environment. 

The specific objective is to develop real-time computer control 
programs for the JPL/CURV manipulator referenced to proximity and force- 
torque sensors. 
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The specific objective includes also that a well-designed software 
basis shall be provided for further research on the following issues: 

( 1 ) Development, testing and improvement of new control 
strategies and algorithms. 

(2) Investigation of the iispact of design parameters including 
the physical characteristics of proximity and force-torque 
sensors and other relevant system components on overall 
control capabilities and on overall system performance. 

(3) Investigation of the istpact of real environment (irregular- 
ities of objects, noise, component impe'*f actions, etc.) on 
control capabilities. 

(4) Study of software approaches to the solution of manipulator 
control problems. 

This report is the first quarterly progress report in which the 
basic concepts of control algorithms and their software Implementation 
are considered. The programs developed thus far are described in the 
program documentation given in Appendix B of this report, which has been 
Issued as a separate volume. 

Before starting the system description, the functional require- 
ments of the CACS software will be reviewed. 
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SECTION II 


FUNCTIONAL REQUIREMENTS 


Th« class of manipulacor tasks, which the CACS software has to 
support, can be defined by the following control functions: 

(1) Unconstrained control of the arm, hand and jaw in 
Cartftslan coordinates (world space). 

(2) Constrained control of the manlpalator, i.e., moving the 
hand with the erd-effector (EE) at a constant distance 
above a fixed tracking plane (TP). 

(3) Tracking a target object (TO) which is slot^ly moving on the 
TP with constant speed, with arbitrary orientation and with 
straight line trajectory. 

(4) Grasping and stopping the TO ^Ich is stationary or moving 
on the TP. 

These functions must be performed interactively from the opera- 
tor's control console (universal control p^uiel) which is specially 
designed for this purpose. 

To facilitate the interactive manipulator control, the following 
automatic operations should be supported by CACS software. These opera- 
,tlons are as follows: 

(1) Roll and pitch alignment of the EE to TP. 

(2) Tracking (identification) of TO speed. 

(3) Yaw alignment of the EE to the TO. 

(4) Centering the EE to the TO for best grasping. 

(5) Grasping and stopping of the TO with complian''.e to its 
motion dynamics. 

All manipulator control functions should be capable of being per- 
formed Independently, or to be Imbedded in one continuous sequence of 
operations. Such a sequence is shown in Fig. 2-1, where an idealized 
ordering of system states and the corresponding transition o leratlons 
are depicted. This sequence pattern exists if all operatloi.^ are 
successfully accomplished and all system states are stably attained, 
if some algorithm fails or system disturbances occur, the corresponding 
transition operation oust be repeated, and the sequence pattern becomes 
more complicated. 
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SYSTEM STATE 



rig. 2-1. Manual/Autoaatic Tracking/Crasping Operations Sequence 


















Am •<»n In Fif. nil of-nratlona nr« ••qunntini. Thin ftouM 

due they are axecueed one after anocuar, '^ringing the ayatam from 
ataca to acata. taara la a aeed for aoma operationa that can 

auv^nort Che atane Cranaltion proceaa and maintain the atablllti of 
attained atacaa. Thaae operationa ar« not ^hovn in Fig. 2>1, but ahculd 
be perfonafti almultMeoualy with the aequenMal operations. Sxaa^lea 
are: 


(1) Maintenance of conacant dlatanre of EE from TP. 

(2) Malntanance of conacant distance cf EZ from TO. 

The flrat operai.ion 3tarta ioBUdiaLaly aa EE io roll and pitch 
aligned tj TP, i.u. , whan “.he manipulator begloe the conatrained moclun 
over the TP, aud tamlnacea tdten TU ia stepped, ‘^e aecona operation 
atarta when the apecd cf TO ia Idantlflad, and CazmJntiCaa wben EE atarta 
the praaplng procedure. 

Automatic operationa ulp “.he operator to carry out manipulator 
control caaka. However, air jatlona can occur when actomatlc operationa 
require operator'a aaaiatance. For example, the yew ill^ancpc algorithm 
can have difficultiea rauaed b/ an inconvenient angle of TO tiajectory 
and/or apeed. In that caae, it ahould be allowed for the operator to 
control Che hand anglea directly t;om hia control panel. The: it are many 
almilar examplea. Therefore, to faclllvate or apeed up the Swtomacic 
operationa, two new operationa are added to the oparacion list given 
above: 

(1) Operator'a manual caaiatence. 

(2) Operator'a emergency atop. 

The flrat operation enablea the operator to iaaue to the manipulator 
conmanda which con be erecuted alnultnneoualy with the ongoing auto- 
matic operationa, without interrupting them. The ascond operation allows 
the operator to put the system in a hold state for a certain time period 
to do off-line interventiona without aborting the idiole control process. 
Both operatloiu are parallel to other operations. 

The operationa dlscuaaed thus far define the basic functional 
requirements of the CACS software related to the basic interactive con- 
trol of the manipulator. Howevet, th>«re are also other capabilities 
which are provided through this project. With reference to the objec- 
tives given in S^ictlon I Introduction, the CACS is conaidwred a research 
and development syetern rathsr than a final product. This i]q>oses addi- 
tional requirements which extend the control features. The moat impor- 
tant adblMonal capablliciao are the following: 

(1) Centralized paramator maintenance. 

(2) Automated eyatem taatiiig. 
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Th« first cspAblllcy providss s fc.sc, sasy and rallsbls sscclng and 
aodlflr.ation of syscaa paraaacars, whlla cha second provides a fasc and 
coop rehens ivA syscesi ceecing after hardware /sofevare changes. Both 
feaciuea are in^lenented as an incegral pare oC cha CACS sofeware. 
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SKCTZON HI 

OPEBAZZOMAL OESCRZPIIOIt 
A. GENERAL SYSTI:M DESC&XPTIOK 

Ao ovsrvlcv of ch« CACS !• jlvaa In Fig. S'*!. As ssen, tits systca 
sr.tlojs Uts following hssdwsrs ualcs: 

(1) tfsiil|>ulstcr. 

(2) Operator control consols. 

(3) Control cocputsr. 

(4) CoaqmCsr consols. 

C5) Alsra d'lsplsy. 

( 6 ) Cooputsr Intsrfscs. 

A brlsf dssc..’iptlon of ;^sse u ilts Is givsc In ths following six 
psrsgrspbs. 


a 


a 


( I 


1 . Msnlpulstor 

Ths ssnipulstor conprlsst ths JURV llnksgt srw with E£, sct\istors, 
ssrvD-potsntloMtsrs scd ssnsor systsns. Ths Isttsr srs forcs*-tory)us 
snd proximity ssnsort supportsd by corrsspondlng slsctrcnlcs. Dstslls 
of thsss coo^onsnts srs g.<ven In Rsfa. 1-3. 

Ths outputs of ths nsnlpulstor csn bs divided into two groups of 
dsts: ssnsor dsts sad Jclnt position fssdbsck dsts. The sensor dsts 
srs represented by tbs vectors w ■ • • • <vg) snd I “ (C|.S 2 . > . . .24) 

which correspond to force-torque sno proximity date respectively. The 
wsvy sign over the variables denotes their “rsw” values which must be 
converted Into their corresponding asthaastical values by s cslibrstlon 
procedure. Ordering of proximity sensor date Is as follows: 

S]^ - front left 

02 - froat right 

03 - lower left 

S 4 - lower right 
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CONSOLE 



Fig. 3-1. Data Flow Diagram of CACS 






The joint position feedback vector 6p ■ (Qpi »^F 2 » • • • »®F 7 ) represents 
the servo-potentiometer readings of joint angles: 

- arm azimuth 

62 - arm elevation 

63 - arm extension 

6^ - hand azimuth 

6^ - hand elevation 

0g - hand twist 

Qy - gripper opening 

2 . Operator Control Console 

The operator control console is an Interface between operator and 
control computer, which enables the interactive control of the manipula- 
tor. The main features of the control console are two joysticks for arm 
and hand control, potentiometer for EE control (opening/closing the jaw), 
and a group of functional switches and pushbutcons. 

The joysticks provide easy position or rate control of three inde- 
pendent space coordinates where the magnitude of posltlon/rate Is pro- 
portional to stick deflection. Details of joysticks, potentiometer for 
EE control and other features of the console are described in Refs. 1 - 3 . 

Joysticks and potentiometer for EE control generate an output 
vector u which has two versions: 


u = (Ax ,Ay ,Az ,Ao ,AS »Av »Ag ) - for position control 

^ c cccc c c 

u a (ij.,yg,z^.,a^,^g,Yc»gc) - for rate control 

where AX(,,Ayg,AZc, and Aaj,,A6c,AYc ®re commanded translational dis- 
placements of the arm and rotational displacements of the h^d respec- 
tively, both in world coordinates. Variables Xj,,yj.,Zg and aj.,Sg,Yc 
aie translational and rotational speeds of the arm and hand in the same 
coordinate system. Values Ag«. and represent commanded increments 
and commanded speeds of the jaw opening, respectively. 


Functional switches and pushbuttons are new features which are 
peculiar to CACS, and which have to be added to the operator console. 
During the development phase, these switches will be Implemented on the 
computer panel. The list of all switches is given In Table 3-1. The 
following comments will help the understanding of this table. 

If the switch SMC Is In the "off" position, the operator has full 
manual control of the manipulator. The commanded position/ rate values 
will result In arm/hand movement and jaw opening/ closing, which Is 
unconstrained within the isanlpulator's motion envelope. Choice of posi- 
tion or rate mode of control is made by switch SRT. 

If the switch SMC Is In tne "on" position, the manual control of 
the manipulator becomes bounded to a fixed plane which has been defined 
previously (default value Is x-y plane In world coordinate system). In 
this case, the control system will Ignore the commanded values Az^ or 
z^, and the resulting motion of the arm will be constrained to the 
given plane. For Instance, the distance of the arm from a given plane 
will automatically be maintained at a constant value. Commanded values 
Ax^, Aye or xc, ye will be taken as displacements or speeds In the new 
coordinate system defined by the plane. 

Switches SAT, SAY, SAC and SAG define the corresponding automatic 
operations If the required conditions are met. For exaotple. If the 
switch SAT Is set, the control system will automatically take the con- 
trol over from the operator when the front proximity sensors register 
a preselected "proxl^ty distance" from the TO. 

By switch SSC the operator can start an automatic search for the 
TO by scanning with the EE in the work space. Parameters of this 
operation must be Inputted previously from the computer console noting 
the parameter editing procedure. 

Switch SOS enables operator interference during automatic opera- 
tions. For example, the operator can adjust the hand angles simultane- 
ously with automatic tracking of TO to accelerate this action or, to 
support this action conq>letely. Without setting this switch, no joy- 
stick commands will be acknowledged during automatic actions. 

By switch SIM the operator can reinitialize the system, l.e., to 
put the system In the Initial state as It was in the beginning of the 
manipulator operation sequence. After, reinitialization, the system 
starts from the beginning. 

Setting the switch SES, the operator Immediately stops the 
movement of the arm and puts It In the hold state. The arm can again 
be released only by resetting this switch. 


Table Functional Switches 


Symbolic 

Name 

Functional Description 

SIN 

Initialization of the system 

SMC 

Manual constrained control (searching for TO 
in tracking plane) 

SAX 

Autonatlc tracking 

SAY 

Automatic yaw alignment 

SAC 

automatic centering 

SAG 

Automatic grasping and stopping 

SRI 

Rate control 

SSC 

Automatic searching by scanning 

SOS 

Operator's manual assistance 

SES 

Operator's emergency stop 


The list of pishbutcons is given in Table 3-2. 

By pressing tht pushbutton TSE, the operator can release the arm 
from the blocked posltiuu v,"hold state"). Namely, during manual control 
of the manipulator, the computer control system automatically puts the 
arm in the hold state when the proximity sensors detect the proximity 
of the TP or the TO. In this case, an alarm signal will warn the opera- 
tor who must conmtand an appropriate motion to the arm to avoid collision 
with TP and/or TO. Before Issuing corrective commands, the operator 
must release the arm by pressing pushbutton TSE. 

If the pushbutton TPC Is pressed, the computer control system 
will Immediately memorize the current values of the hand coordinates In 
order to define the orientation of the tracking plane for subsequent 
constrained control. This action will always be taken when TPC Is 
pressed. This means that the TP coordinate setting can be done more 
than once. 

By pressing the pushbutton TCA, the operator can turn off all 
alarms displayed at that time. 

By pressing pushbutton TPA, the operator invokes a parameter 
updating procedure. Namely, all parameters are stored In a particular 
redundant storage area which Is accessible by the parameter editor sub- 
system EDIT (see Subsection C paragraph 1) concurrently with manipulator 
operation. After updating the parameters by EDIT, the new parameter 
values must be passed to the corresponding parameter locations of the 
control subsystem (OPER) . This Is automaclcally performed by special 
transfer procedures which are executed Immediately after pressing the 
TPA pushbutton. 


3. Control Computer 

The control computer Is the heart of the CACS. It consists of 
an INTEKDATA M70 minicomputer and the CACS software package which runs 
under OS/16 MI2 real-time operating system. The CACS software contains 
three subsyscems: testing subsystem (TEST), parameter editing sub- 

system (EDIT), and manipulator control subsystem (OPER). The former 
two subsystems are described In Subsection C paragraph 1 of this report. 
The latter subsystem Is the main part of Che CACS software. All three 
subsystems are Implemented as different tasks. 


4. Computer Console 

Communication with the operating system and with the CACS soft- 
ware Is through the computer console. It Is a teletype unit, but can be 
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Table 3**2. Pushbuttons 


' 1 V 


Symbolic 

Name 


Functional Description 


any other suitable peripheral device. Connunlcatlon with the TEST and 
EDIT subsystems Is only possible through the computer console, while 
communication with the OPER subsystem Is through the operator augmented 
control console. 


5. Alarm Display 

The purpose of the alarm display la to warn the operator in case 
of system Irregularities, system abort, or whenever his assistance Is 
needed. The alarm display can be Implemented as a mosaic field of 
labeled lights or as alphanumeric messages on the CRT display. During 
the development phase the latter will be used. The complete list of 
alarms Is given In Table 3-3. 


6 . Computer Interface 

Computer interface consists of analog to digital (A/D) and digital 
to analog (D/A) converters. The Input to the A/D converter Is the com- 
pound vector dp, u, w, £)^ while the Input to the D/A converter Is 
the joint position vector^ » (§j^, 02 ,..., § 7 ). 
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Table 3-3 . Alarna 


Symbolic 

Name 

Cause/ Action to be Taken 
by Operator 

ALIM 

Arm/hand on the boundary envelope. 

AIFP 

EE In the proximity of TP. System la In HOLD 
state and operator has to release It by pressing 
THE In order to move arm sway from TP. (Alarm 
Ignored If SMC set.) 

ATOP 

EE In proximity of TO. System Is In W)LD state 
and operator has to release It by pressing TEE 
In order to move EE away from TO. (Alarm 
ignored If SAT set, or Sm reset.) 

AYAH 

Angle between motion vectors of EE and TO Is 
greater than 90 degrees. Automatic yaw align- 
ment not possible. Operator assistance Is 
needed. 

ATOL 

TO lost, control Is given back to manual con- 
strained searching. 

ARPL 

Roll or pitch alignment lost. Operator assistance 
required . 

ATOG 

TO too large, cannot be grasped. 

ASTP 

TO grasped, but cannot be stopped by given force- 
torque limit. 

AGRT 

Grasping terminated. Procedure must be repeated. 

ASCI 

Search by scanning has been terminated. 

ASCC 

Search by scanning has been completed 
successfully. 

AGOL 

EE in collision with TP or TO. Automatic action 
of moving EE one step from TP or TO has been 
taken. Operator's assistance may be needed. 
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CONTROL ALGORITHMS 


In order to develop the system described In Section III with the 
requirements In Section II, a variety of control algorithms has been 
and will be developed. The complexity of the system requires a struc- 
tural approach to the development of the control algorithms to provide 
conditions for an easy functional decomposition of the system and to 
assure easy and stable system Integration. This Is especially Important 
In a laboratory development environment In which a stepwise development/ 
refinement Is required, as Is the case with CACS. 

As guidelines for a structural approach to the control algorithms 
development, the following three concepts are Introduced: 

(1) Incremental motion synthesis. 

(2) Basic control routines. 

(3) Parallel process concept. 

All these guidelines and their application will be considered In the 
next six subsections. 


1. Incremental Motion Synthesis 

The Idea of Incremental motion synthesis requires the production 
of Increments of motion rather Chan "endpoint" values In each system 
iteration cycle. The principle Is depicted In Fig. 3-2 (graphical 
symbols used In this picture are explained In appendix A) . Vectors dx 
and V represent incremental translatory displacement and tracslacory 
speed which must be performed by the manipulator EE^ These quantities 
are effected by the control algorithms to perform Che required automatic 
operations Initiated by manually Issued comnands. The position Incre- 
ment Is added to the content of the buffer TRP (translatory-posltlon) . 
The speed veccor Is converted Co a corresponding rate Increment Ax^, by 
multiplying with the clock interval At, and stored In the buffer TRR 
(translation-rate). In each cycle the contents of both buffers are read 
and summed. The resulting increment Ax^^^ represents a composed motion 
increment which must be transformed InTS a correspcndlng Increment 
in the joint space coordinate system . This value will be added to 

joint position vector which has been generated in the previous 

iteration cycle. The new vector will be used in the current 

iteration cycle. The essential point here is that the position buffer 
TRP must be reset after reading, while the content of the r-'te buffer 
TRR remains unchanged. This will cause arm/hand movement with constant 
speed V and. simultaneous arm/hand displacement Ax. The latter will not 
occur in the next iteration cycle, if the position buffer Is not filled 



Fig. 3-li. Principle of Incremenf:al Motion Synthesis 





again. In ordar to axacuca -Isplacamaot Lx in ona itaration cycla, cha 
nagnlcuda of this vacto** muat ba laaa than tha flxad valua datarminad 
by tha itaration cycle and by the oanipulator dynamic parfontanca. Ir 
this is satisfied, the position displacement of the manipulator can ba 
cousldered as in inatantanaoua displacemant If the position increment 
doer not satisfy this raqulramant, tha corresponding operation muat ba 
delayed. This will bo dirrusaed latter. 

The Incremectax approach enaoles easy coordinate trana format Ion. 
Ujually Che eq-Mtlon of joint variables are given in the following form; 


X “ 1(®) 


( 1 , 


•#here XCfi) represents a vector mnpp'ng of Joint coordinates into tha 
world (Cartesian) coordinate system. Using the incr^ntal form of 
variables, (1) can be written: 

For sufficiently small increments , the Taylor expansion of (2) 
gives ; 


^x(k) 3j g(h) ^0(k) 


(3) 


or; 




(4) 


where: 




36 


e»e 


(k) 


( 5 ) 


represents the Jacobian of the transformation (1) 


Having («iia In Bind, aa wall as cha da^crlption glvan abova, cha 
advancagaa of ta Ircraaancal motion ayncnasla principle can be acs- 
aarlzad as followt: 

(a) It tics tha nature of cyclic procasacs and, tharafora, 
anatlas a more convaniant ImpleBancatlon chan cha endpoint 
appioach. 

(b) It a low>. unified handling of position and rata control 
modao. 

(c) It enables easy ayntUa->ls of composed aocions of both 

types: tosition-poalcion and positlon-raca. (An example 

of a combined ss>ticn of tha posltlon-ratr. type is given in 
F^g. .^-3, where the EE which Cracks tha target object with 
Che seae spaed nnisc corract its relative position hn>) 

(d) It provides easy coordinate transformations. 


2. Aprli>ation of the Incremental Motion Synthesis in the CICS 

The prin<:iplc of incrementnl motion synthesis described in che 
previous subsection ie applied to the CACS. The main features of the 
application are given in Fig. S-A. This diagrar is derived from a 
basic property of the CURV arm: hand orientation is independent from 

arm elevation and extension. This is due to the double parallelogram, 
mechanism added to che linkage. Therefore, che coordinates of the arm 
<nd h'.nd can be handled reparately, l.e. , tne position displacement and 
speed composition shown in Fig. 3-2 is now directly applied to the 
arm position coordinates x “ (x,y,z) and lo the hand orientation angles 
■ (a,S,Y)« For that reasou two new buffers arc introduced: ROP 

'Irotat ion-posit ion) and ROR (roratlon-rafe) . The jew operation is 
considered only as a matter of position control. Therefore, only one 
buffer JAW is introduced for this coordinate. 

Fig. 3-4 will now be explained step by step. Raw values read 
from A/D buffer are grouped into four vectors which are already 
described In Ssetion IIIA. All these variables are calibrated separ- 
ately. Calibration oi joint position feedback is done Ly che followlrg 
linear equations: 


* ®Fi * ^i ” ^i’ ^ " i«2,...,7, (6) 

where a^;^ are scale factors of dimension rsd/Volc, and bi are zero offsets. 
These coefficients will be considered as system parameters. 


3-13 



SA'JTCH/POSHWJnON 



Fig. 3-4, General Data Flow Dlagraa of CACS 








Ig. 3-4. General Data Flow Diagram of CACS (Continued) 
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A similar calibration procedure is applied to vector u. Sensor 
data are calibrateo by special table look-up procedures which will not be 
considered here. 

Calibrated variables are used by control algorithms to produce 
displacement increments Axp ■ (Ax_,Ayp,AZp) , Aop ■ (AUp.ASp.Ayp) and 
Ag, and velocities v ■ (u,v,w) ana u • l.e.. the correspond- 

ing rate Increments > (Ax^^Ay^ and Ao^ « (Aaj.,AS^,AYr) . These 
valiies are stored in the buffers TRP through JAH. 

Because of the finite dynamic response characteristics of the 
manipulator, the rate Increments must be limited. Therefore, the fol- 
lowing value will be stored in the buffer TRR: 



if 


max 


|Ax^,Ay^,Az^| < 


max 



otherwise 


( 7 ) 


where: 




and r . is a given parameter and represents the maximum possible value 
of anf^isplacemcnt that can be achieved iu one clock cycle. This is 
the condition for the rate buffer input operation, which prevents getting 
the manipulator into unpredictable working conditions. Similarly, the 
rotation rate Increueuts are also limited. However, they are limited 
component by component, because in the case of rotation it is not Intpor- 
tant to keep the proportloncilty of the vector components. 

In each iteration cycle ail buffers are read, the composed values 
Ax^ and Ac^ are formed and brought to the coordinate transformation 
block. 


i 




I 


i 


1 






c 


i 

i 


; 


r 


Coordinate transformations are given by the following equations: 


Ae . = T, (0 . ) • Ax 
-A A --A -IT' 




( 9 ) 


t 


j 



r 


1 


- kj • is . 


where: ^ ■ (01f02»*^3^ ^ * (<54»05»96) hand joint vari- 

ables, respectively. 

Tiansformation matrices Ta end Tg are Inverse Jacobians of joint 
to world space transformations given in Ref. 3. These matrices have 
the following elements: 


where: P » a^^ + r3S2, t3 » 2DS3 + d, - sin 63, c, = cos 03, 

S2 * sin 02» ^2 “ ®3 “ (63/2), C3 ■ cos (03/2), D, d, kj 

cad 33 are mechanical parameters of the manipulator, and I is the unit 
matrix. The most recent feedback values can be used as the current 
values of joint variables 63,62 and 63. 

The r.rm and hand coordinates are not completely Independent. 

It can be seen from the equation for the hand azimuth (Ref. 3): 


»-«l + «4-2 • 


Therefore, the following correction must be done for A6^, after the 
coordinate transformation has been completed: 


A0. = A0' - A0- . 

4 4 1 


( 
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Finally, the transformed and corrected joint variable Increments 
are joined into one vector A 6 (^), where the index k denotes the current 
iteration cycle. This vector is added to the joint position vector 
0 (k-l) from the last Iteration cycle to form the new joint position 
vector 9 ^^). This vector is then converted to ‘'raw** values and is 

brought to the D/A buffer DACBUF. The conversion to raw values is an 
inverse to the calibration procedure defined by equation ( 6 ). 


3. Basic Control Routines 

The algorithms which perform the operations of roll/pitch/yaw 
alignments, tracking, centering, grasping and stopping of TO as well as 
other interactive operations specified in Section II, are based on a 
variety of specific motions of the arm, hand and/or jaw. A look-ahead 
study of the whole CACS has shown that all foreseen operations can 
successfully be carried out by a unique and finite set of elementary 
manipulator actions. These actions will be implemented as a set of com- 
mon routines which can be used by all CAC^ a - /"'rithms. The procedures 
related to these actions are listed in Table 3-4. Their explanations 
are given in Figs. 3-5 - 3-9. The following comments will supplement 
the explanations. 

Rotational motions (changing the yaw, pitch and roll angle) are 
effected by three Independent routines, YAW, PITCH and ROLL. This is 
for the sake of simplicity of algorithm implementation, since different 
kinds of angle changes are used in the various algorithms. 

Motions referenced to the hand coordinate system (Fig. 3-7) 
require additional coordinate transformations. The rotation matrix 
which transforms the hand coordinate system into the world coordinate 
system can be represented in a block form as follows: 

i 4 j 4l 
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Table 3*4. Basic Control Routines 


Routine 

Identifier 

Description of the Operation 

SHIi'T 

Translatory displacement of the arm one step 
In the world coordinate system hx ■ (Ax^Ay^Az). 

SHIFTC 

Constrained translatory displacement of the arm 

one step in the tracking plane ■ (Ax^,Ay^). 

"“t t t 

YAW 

Rotational displacement of the hand one step in 
the world coordinate system Aa (changing the 
yaw angle) . 

PITCH 

Rotational displacement of the hand one step 
in the world coordinate system A8 (changing 
the pitch angle). 

ROLL 

Rotational displacement of the hand one step 
in the world coordinate system Ay (changing 
the roll angle). 

MOVER 

Translational displacement of the arm one 
step in the hand coordinate system AC (longi- 
tudinal motion of the hand) . 

MOVEE 

Translational displacement of the arm one 
step in the hand coordinate system Ari (lateral 
hand motion of the hand) • 

MOVEZ 

Translational displacement of the arm one step 
in the hand coordinate system At (vertical 
motion of the hand) • 

EXPND 

Expansion (opening) of the jaw for one incre- 
ment Ag« 

CONTR 

Contraction (closing) of the jaw for one incre- 
ment Ag. 

RUN 

Motion of the am with the constant speed 
V » (u,v,w) in world coordinate system. 

RUNG 

Motion of the arm with constant speed 
» (u^, v^) in tracking plane.. 

ROT 

Rotational motion of the hand with constant 
angular speed ^ « (a,3,y) in the world 
coordinate system. 


Note: All translational motion of the arm assumes unchanged hand angles » 

and all rotational motion of the hand assumes unchanged arm 
coordinates. 
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Pig. 3-5. Translatory Motion of the Arm in World Coordinate System 








Fig. 3-6. Rotational Motion of the Hand in World Coordinate 






Translatory Motion of 









Constrained Translatory Motion of the Ara In Tracking Plane 







e 


where the submarrlcee (1 ■ 1»2,3) are the followlag column vectors: 


C 


-2 


cos6 cosa 

-slny sin6 cosa *f cusy slna 

cosy sln0 coco + slny alno_ 
-cos(t slna 

ciny sinB fllaa + cosy coso 
-cosy slnfi 8ina_*f slni^cosa 


(15) 


-3 


-slnfi 


-slny coffi 


cosy cosfi 


The angles a, fi and y are the 'iurrent values of the hand angles In the 
world coordinate system expressed by Joint coordinates: 


I 


i 


« “ - f 


6 ■ 0 


5 


IT 

2 


Y ■ 0^ . (16) 

The most recent feedback values can be used as actual values of Joint 
coordinates. 


C 


Con«cr»lned cranslatory motion in tha cracking pl^na also naada 
additional coordinate ccan8fonDacloaa4 Hara» the ttanaJomatlon matrix 
has tha form: 


^t * -1 
-2t “ -2 


where r\ - (ac*dt>Yc^ repreaen«.a tracking planr. angles, i.e., the hand 
angles taken In tae TP-coordlnace setup procedure (invoked by presaliig 
pushbutton TPC). As will be shown in Section IIIO, elements of 
matrix are calculated during the coordinate r tcup procedure and 
ntemorlzed In the special buffer (AT) cor subsequent ure. 

Using the Basic Cont*'ol Routines, an extended sec of n^'f routines 
can be derived for toore complex manipulator operations. For exnaple, 
the centering algorithm can be simplified considerably If asymmetric 
extension/contraction jaw operations are ir croduced Therefo e, four 
new routines are added to the list In Table 3-4: 

EXFNUR - Expansion of the jaw for ^ue increment Ag to Che right. 

CONTRR - Ccntractlon of Che jav for one increment Ag from the 
right. 

EXFNDL - rxpanslon of Che jaw fcr one increment Ag Co the left. 

CONTPL ' CortracClon of the jaw one Increment Ag from Che left. 

r..e first two routines assume opening/closing of the jaw with 
fixed left aid . and right slda mowing Ag out/in. The lest two routines 
are opt-oslte to the operations implied by Che first two routines. These 
ro<\cines are shown in Figs. 3-10 - 3-11. 


4. Dynamic Response of Manipulator and Delay 

Incremental manipulator JisplacemenCs, commar.ded by CAC£ software, 
are r<dalized by relatively complex electro-hydraulic servo systems. 
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Fig. 3-10. Asyanetrlc Expansion and Contraction of the Jaw (Fixed Left Side) 











c 


c 


Because of limited dynamic perfomances of these servo systems, the 
aajustments ot joints, l.e., the realization of corresponding displace- 
ments In the vrorld coordinate system, requires some finite time. There- 
fore, the displacement magnitudes must be In accordance with the manipu- 
lator servo time constants and with the Iteration frequency of the 
computer system. This means that the commanded displacements must be 
sufficiently small so that they can be realized In one Iteration cycle. 
Let ^t be the time Interval of one Iteration cycle, and Ax mav and Aamjy 
the maximal values of translator^ and rotational displacements of the 
manipulator which can be realized within the time ht. Then all compo- 
nents of the vectors * (AX(.,hyc,AZ(,) and • (AO(.,AB(.*^Yc) must 
be less than or equal to the values Axq^x ^°hsax> respectively. 

Of course, this Is a rough consideration since displacements executed In 
a given time Interval can differ from case to case, depending on the 
current state of the manipulator and on the attached load. But It Is 
acceptable here as an Initial step which can be refined later. 

As will be seen. In some cases It will be necessary to realize 
displacements greater than the limiting values defined above. In such 
cases, delay mechanism must be Introduced. If, for example, the value 
Ax > AXg^]( has been Issued by one algorithm. It mast wait for a next 
decision until the displacement Ax Is completely realized by the mani- 
pulator. This situation is Illustrated In Fig. 3-12. As seen, the 
displacement will be completed during the third Iteration counting from 
the moment the command has been Issued. It means that the algorithm 
must be delayed two Iteration cycles. In general, the number of delayed 
cycles n]) can be determined by the following approximation: 




max {Ax ,Ay ,Az } 
c c c 


Ax 

max 


1 


(19) 


where fal denotes upper integer value of **a'*, and is a given 

parameter and represents the maximum possible value of any displacement 
that can be achieved in one clock cycle. The same formula is used for 
the rotational displacement A^ » (Aa^, ACc>^Yc) • This formula is derived 
under the assumption of a linear shape of the dynaLiic response of the 
manipulator. Implementation of the delay mechanism will be explained in 
Section IIIC* 


5. Parallel Processing Concept 

As shown in Section the entire process of searching, tracking, 
grasping and stopping of the target object can be broken down into a 
sequence of small units of activities or operations. Some of the unit 
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activities are carried out simultaneoualy, i.e., they are executed in 
parallel and asynchronously. This is the case with the following 
activities: operator manual assistance, operator emergency stop, 

input-output operation, automatic maintenance of the distance between 
EE and TP, and automatic maintenance of the distance between EE and TO. 
These activities are performed independently from and simultaneously 
with other activities which are being executed sequentially. 

The structural approach to the algorithm design and the software 
implementation of the CACS requires that all activities be considered 
as separate programming modules. Furthermore, the parallelism of the 
activities requires the concurrent execution of these programming mod- 
ules. An additional benefit of the concurrent programming approach is 
that more than one processor can be allocated to the system to increase 
execution speed and reliability. This is especially important in the 
case of real-time environments expected for future versions of manipu- 
lator control systems. Of course, the concurrent programming system 
can be realized with only one processor which is multiplexed among con- 
current programs. The scheduling of one or more processors for con- 
current programs will be achieved through system software not discussed 
here. The structure of the control algorithms and their software 
implementation considered here are invariant to the number of processors 
which can be Ignored in further considerations. 

The program which is being developed to accomplish the specified 
activities will be called process . A process can generally have two 
states: active and blocked . Active state means that the process is 

in execution. (The active state has two substates: running when the 

process is being executed by the processor; and ready when the process 
is waiting to he allocated to the processor. But, this Is not important 
for the following discussions.) In the blocked state the process is 
waiting to be resumed, or it is waiting for a specified time interval 
by which it has been delayed. Putting the process in active or blocked 
state is done by special synchronization signals. By means of these 
signals a process can be waiting or delayed, or it can resume another 
activity if some specific conditions occur. 

For a concurrent programming system it is important to describe 
the interprocess relations in terms of synchronization signals. This 
can be shown graphically by process precedence charts which are intro- 
duced here to simplify the discussion. Cn these charts. Fig. 3-13, 
processes are depicted by circles which are connected by arrows represent- 
ing synchronization signals (no data flow!). Besides the arrows, the 
synchronization conditions are indicated as boolean variables or 
expressions. In the case of regular arrows, signals will occur if the 
condition indicated is true. In the case of dotted arrows, synchroniza- 
tion sign'ils will occur if the condition is not satisfied. Dotted arrows 
will be used to emphasize "backward signalling", that is, when the 





H 
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i 


; i 
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process loses the active state because the required condition Is not 
satisfied anymore. In order to preserve system reliability, the follow- 
ing general rule will be accepted here: A process can be made to wait or 

delayed only by Itself . Therefore, only process resumption Is presented 
by arrows. 


At the right side In Fig. 3-13 an example of a process precedence 


chart Is given. It shows the process being resumed first. If during 


Its execution condition C]^ occurs, the process ?£ will be resumed and 
will be waiting Immediately. If this condition Is changed to false. 


process ?2 will be waiting and Px resumed again. If ?2 Is active and 


condition C 2 occurs, two processes will be resumed at the same time: P 3 
and P 4 . These processes are executed concurrently, together with 
process P 2 , until the condition C3 holds. If condition C 3 changes to 
false, P 3 will be waiting and Px will be resumed again. P 4 , once 
resumed, remains active until C 4 becomes false. 


Besides the synchronization signals, data Is another reason for 
Interprocess communication. Some processes produce data while others 
use data. The data used (accessed) by several processes are called 
shared data . In order to preserve system reliability and to facilitate 


data communication, special attention must be given to this question. 
Today's software techniques offer special mechanisms, called monitors, 
for handling both, shared data and synchronization signals (Refs. 4-7). 
By means of such monitors It Is possible to control the access to shared 
data and to enforce various access right policies. Monitors and their 
Implementation, as well as process scheduling, will be discussed In 
Section me paragraph 6 . 


6 . 


Application of the Parallel Processing Concept In the CACS 


The list of all processes of the manipulator operating subsystem 
OPER Is given In Table 3-5. It has been made using Fig. 2-1, Section II. 
As seen, five new processes AMP, AMO, lOP, OMA and OES have been added. 
The first three processes are unconditionally active when the OPER sub- 
system Is active. Therefore, two groups of processes can be dis- 
tinguished: permanent active , and temporary active . A new rule will 

be added to the one mentioned In the preceding paragraph: permanent 

active processes cannot be waiting or delayed . This rule will be 
enforced by the scheduler as will be shown In Section IIIC paragraph 7. 


The process synchronization conditions will be handled by two sets 
of boolean variables: functional switches and event flags . The former 

is already described In Section IIIA paragraph 2, while the latter is 
given in Table 3-6. The event flags describe the state of the manipu- 
lator. They are updated by an input-output process which is a permanent 
active process. More details about event flag generation are given In 
the program documentation In Appendix 3. 







Table 3-S. CACS Processes 


Identifier 

Description 

lOP* 

Input-output process 

OMA* 

Operator manual assistance 

OES* 

Operator emergency stop 


Unconstrained manual control 

ARA 

Automatic roll alignment 

APA 

Automatic pitch alignment 

MCS 

Constrained manual control (search for TO) 

AMP 

Automatic maintenance of distance from TP 

ATS 

Automatic tracking of TO speed 

AMO 

Automatic maintenance of distance from TO 

AYA 

Automatic yaw alignment 

AGO 

Automatic centering of EE on TO 

AGO 

Automatic grasping 

ASO 

Automatic stopping of TO 


“k 

permanent active processes. 


The general outline of all temporary active processes Is given 
In Fig. 3-14. As seen, the processes are cyclic. If a process Is 
active, the switches and activity conditions are checked In every 
Iteration cycle. These conditions are represented by boolean expressions 
of functional switches and event flags, depending on the particular 
process. Synchronization signalling Is performed by two primitive pro- 
cedures: signal and wait. (Signalling procedures put In parentheses are 

not applied In all processes.) Completion of the activity to which the 
process Is dedicated Is also determined by event flags. If the activity 
Is completed, special state variables must be updated. These variables 
will be discussed later, together with the algorithms implemented by 
the processes. 

The Interprocess relation Is described by the process precedence 
chart In Fig. 3-15 which comprises all temporary active processes. 

This chart Is self explanatory by looking up the condition variables 
listed In Tables 3-1 and 3-6. Hyphenated event variables represent a 
logical condition or operation. For example, ef-t replaces the logical 
expression eflt o£ efrt. 
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Fig. 3-15. Process Precedence Chart of the OPER Subsystem 




Table 3-6. Event Flags 


Identifier 

* 

Description 

CALM 

Arm has reached boundary envelope (joint 
variables ^3 ^ position). 

EHLH 

Hand has reached boundary envelope (Joint vari- 
ables and/or dg In end position). 

EFLP/ETRP 

Front-left/right proximity sensor Indicates 
proximity distance. 

EFLT/EFRT 

Front-left/rlght proximity sensor Indicates 
tracking distance. 

EFLC/EFRC 

Front-left/right proximity sensor Indicates 
collision distance. 

ELLF/EL8P 

Lower-left/right proximity sensor Indicates 
proximity distance 

ELLT/ELRT 

Lower-left/right proximity sensor Indicates 
tracking distance. 

CLLC/ELRC 

Lower-left/right proximity sensor Indicates 
collision distance. 

ESOL 

EE Is roll aligned to TP. 

EPIC 

EE Is pitch aligned tu TP. 

ESPD 

Speed of TO Is attained. 

EYAW 

EE Is yaw aligned to TO. 

ECNT 

EE Is centered on TO. 

EGSA 

TO Is grasped. 

ESTP 

TO is stopped. 

EJCT 

Jaw closed to the cracking aperture. 

EJCL 

Jaw closed. 

EJOP 

Jaw open. 


description holds for flag true. 
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SOrmRE IMPLEMENTATION 


In prnvloun subsections the CACS bes been described from func* 
tlonel, hardware and algorl'ihmlc points of view. Now, the software 
Implementation of the system la discussed. The general software archl- 
cecturo, the basic software components and the scheduling of the proc- 
esses are considered In the first seven subsections. The processes 
themselves are not designed at the time of writing this report and will 
not be considered here. The description of the program documentation is 
given In paragraph 8 of this section. 


1. General Software Architecture 

As elready mentioned In Section III A paragraph 4, the CACS soft- 
ware consists of three parts: 

(a) Manipulator control subsystem (OPER). 

(b) Parameter editing subsystem (EDIT). 

(c) Manipulator testing subsystem (TEST) . 

Thase three subsystems will be implemented as three independent tasks 
under the OS/16 MI2 real-time operating system. 

The first subsystem Is the main part of the CaCS software ^Ich 
supports interactive control of the manipulator from the operator con- 
sole. The preceding s>sctlons as well as the major part of this section 
are devoted to this subsystem. The later two subsysteois are not part of 
this report and will only bclefly be discussed In the next two 
subsections . 


2. Parameter Editor Subsystem 

The Ei)IT subsystem Is Intended to rupport fast and easy editing 
of all CACS parameters and constants. This Is extremely Important in a 
laboratory environment and for the purpose of experimentation as empha- 
sized in Section 2. The relatively large number of parameters and the 
frequent need for their readjustments make this problem nontrivial. The 
following needs must be taken care of: 

(a) Provide facilities for easy parameter changes before or 
during manipulator operation. 

(b) Preserve system Integrity and reliability. 

Both problems can be solved simultaneously by means of a centralized 
and redundant parameter file ^Ich is maintained independently of the 
manipulator control software. This is shown in Fig. 3-16, where four 
routines are outlined. 
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Fli^. 3-16. Syshea Block Dla^raa of Paraaeter Editor Subsyat 



Th« first rojclns psrsMCtt editing provides en on-line updece 
end display of the user (source) peresscsr file. This file contains 
one record for esch parasictet, where each record contains at least the 
followlog data fields: 


(a) 

Paramater identifier (Its symbolic noma 
programs) . 

used In OPER 

(b) 

Index nuad>er for case of parameter arrays (every vector 
coB^onMit represents one parameter entity). 

(ft) 

Descrrlptlon of the parameter limited to 
characters. 

fixed number of 

(d) 

Unit of moasuremsQt. 


(e) 

rarameter value. 


(f) 

Date of the last parameter update. 



The updating procedures should cooprlse addition, deletion and 
o»dlflca*;lon of en**lre rscords and of particular record fields. To 
facilitate the experlnenls, two values can be associated with each 
parameter: the "actual value" and the "try value". The former 

represents the value established In previous experiments, while the 
latter rapr«.sents the new value ready for use In the current experi- 
ment. The try value can alwoys be returned to the actual valu* If the 
user Is not satisfied with Its effect on the system, or It can be fixed 
as a new actual value If It gives better results. In the last case the 
da**e of the last parameter update must also be changed. The user param- 
eter tile can be displayed by a CRT or by console printouts as shown by 
example In Fig. 3-17. 

The parameter conversion routine automatically generates the 
parameter table after the editing procedure Is complete. This cable 
contains paraiseters and their combinations in Che "object" form used 
by the OPER procedures. These values are not necessarily redundant, 
and cannot be accessed Ly the user. Examples of parameter conversion 
are the following expressions: 

AR - (H + DLOWER(l)) * SIN(GAMINC) 

BR - -DBS/2 * (1 - COS (GAMING)) 

CR - DBS/2 * SIN (GAMING) 

’)R • (H -f DLOWER(D) * (1 - COS (GAMING)) 

where AR, BR, CR and DR are "obj .stive" parameter values used directly 
in the roll alignment algorithm, while H and DBS are constructive con- 
stants of the CURV arm, and DLOWER(l) and GAMING are algorlthadc param- 
eterr, all representing redundant source parameters in user form and 
maintained by the user. 
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Pig. 3-17. Exaaple of Paraaeter Display 
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The file creation routine Is used only once. Immediately after 
establishing the EDIT task. This routine Is trivial If the sequential 
organization of the user parameter file Is accepted. 

The save routine Is a system utility routine which provides 
copying of the user tile onto the output media as punched cards, paper 
tape, or magnetic tape. This enables the keeping and maintaining of 
more than one set of paraipeters. For example, each set of parameters 
can be associated with different sensor calibrations or with different 
experimental runs. Before starting the manipulator operations, a par- 
ticular set of parameters must be chosen and the corresponding save 
file Bust be loaded Into the user file. I'he user file can be a direct 
access auxilllary storage (disk for example) or it can be directly In 
the core. The storage required for a user file Is about 12 K bytes, 
and the storage required for the EDIT code Is not expected to be large, 
so that both can be core resident during the execution of the EDIT task. 
The parameter table, which Is an output of the parameter conversion 
routine, contains only one set of parameter values requiring about 3C0 
bytes of computer storage. This table must be transferred from EDIT 
to OPER task, where It will be placed In a shared data buffer. This 
Is discussed In paragraph 6. 


3. Manipulator Testing Subsystem 

The TEST subsystem Is Intended as a maintenance tool for checking 
the OPER subsystem after every hardware change and/or after system 
reinstallation caused by changes In the software. It Is obvious that 
the dynamic environment of the CACf will require frequert modifications 
of the software to iorprove the existing pregrams or to extend the system 
by new features and capabilities. This means that the CACS software 
must be considered as a "living" part of the system, subject to changes 
and continuous growth. To maintain system integrity and efficiency of 
maintenance. It Is strongly recommended that considerable attention he 
paid to the manipulator testing subsystem. It must be developed as a 
sequence of testing procedures which will autoniatlcaily and syste-iftlcaJly 
check all modules of the subsystem, their mutual interactions and their 
interactions with the I/O devices. The sequencial order must be 
designed In a way so that all possible errors in the hardware and soft- 
ware can be located fast and, eventually, without operator assistance. 

The major part of the TEST subsystem can probably be synthesized from 
testing programs which are an outcome of the overall program development 
process. This must be kept In mind when creating the testing programs 
which must be flexible and suitable for future retailorlug. 


4. Monitor Concept 

The basic software components of the OPER subsystem are processes 
and monitors. Processes are already discussed in the preceding Section 
(paragraphs 5 and 6), and their implementation will be considered In 
paragraph 6 of this Section. Monitors are introduced by Brlnch Hansen, 
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Refs. 6-7, and by Hoare, Bef. 8, as a new concept for a hierarchical 
structuring of nmltlprograimiilng operating systems. This concept is 
adopted to develop the CACS software. It will be discussed briefly in 
this subsection to complete the discussion of Section III B paragraph b. 

A more detailed description of monitors and their implementation is given 
in Refs. 6, 7, 9, 11 and elsewhere. 

Processes communicate with each other by sending or receiving 
data. As already mentioned before, these data are called shared data. 

To preserve system integrity, the processes should not operate directly 
on these data. An example of direct operations on data is common data 
structures which enable unrestricted data accessibility by minimal sys- 
tem overhead. However, the reliability of complex systems, especially 
systems with concurrent programs, is highly sensitive to Incorrect data 
communication. Therefore, the accessibility of data, l.e., the access 
rights of the processes, should be carefully controlled. This can be 
done by allowing access to data only by special procedures. Procedures 
provide much safer interface mechanisms than common data structures. 

This is the basic assumption underlying the modem approach to struc- 
tural design of concurrent programming. Kef. 11. It will also be 
employed in the design of the CACS software. 

The shared data and the procedures which can access them, are 
grouped within an abstract object called monitor . These objects are 
usually high level language constructs designed for concurrent pro- 
gramming, Ref. 9-10. Tney also can be used as a powerful concept for 
building hierarchically structured systems in sequential languages. 

This is the case for the CACS software. 

Regardless of the interpretation of monitors, they define the 
following four entities: 

(a) shared data structure, 

(b) monitor procedures, 

(c) Initial operations, 

(d) access rights 

The shared data structure Includes data transmitted among proc- 
esses, as well as other data necessary for the correct functioning of 
''he monitor in the context of a multiprogramming environment. Monitor 
procedures explicitly define all operations which the processes can 
perform on shared data. Initial operations define all operations which 
must be executed in time of creation of the monitor. Finally, access 
rights define all connections of the monitor to the rest of the system 
hierarchy. This is achieved by giving an explicit list of all processes 
or other monitors which can enter the monitor. 


Besides Che control of operations on shared data and Chelr 
accessibility, another feature concerning the monitors Is essential: 
mutual exclusion . It must be ensured that only one process can operate 
on shared data at the same time. This can be accomplished by allowing 
Chat only one procedure within the same monitor can be called by proc- 
esses or other monitors at the same time. If some process Is executing 
a monitor procedure and another process Is trying to call this or 
another procedure within the same monitor, the latter must be delayed 
until Che first process leaves the monitor. Simultaneous monitor calls 
are scheduled outside the monitor by special procedures grouped within 
one program called kernel . The kernel Is a basic part of the system 
software which Implements exclusive access rights of processes and 
scheduling of CPU's and other physical resources among concurrent proc- 
esses. For Its work, the kernel uses Interrupt mechanisms and low 
level communication facilities which are Implemented In hardware or 
lowest level machine software. Therefore, the kernel Is highly machine 
dependent, l.e.. It can be considered as an extension of the machine 
which hides Its details from the user who builds the concurrent pro- 
gramming system. This Is an essential point of the monitor concept. 

As mentioned before, monitors are not only Intended for trans- 
mission of shared data. They are also used for process synchronization 
and for scheduling of physical resources. Although, monitors have a 
structure essentially Identical for all purposes, two general types of 
monitors can be distinguished: buffer monitors and resource monitors . 

The buffer monitors are designated for shared data transmission 
among concurrent processes. Their data structure Includes three parts: 
shared data buffer, full-empty Indicator, and two single process queues. 
Buffers can be designed to handle one (single buffers) or more (multiple 
buffers) data portions, ^klltlple buffers are usually designed as linked 
circular lists, which Impose three more data items in the monitor data 
structure: current buffer length, and two pointers, one for the head 

and one for the tall of the list. (In CACS the simpler single buffers 
are used.) The sizes of the buffers correspond to the amount of data 
to be handled. The data can be arrays, records or sets. (In CACS only 
data arrays will be used.) 

A full-empty Indicator is a boolean variable which tells the 
processes whether or not the buffer contains data. If the buffer Is 
full no data can be transmitted to the buffer, but the data can be 
taken from the buffer. If the buffer is empty, data can not be taken 
from the buffer but can be transmitted to the buffer. Single process 
queues are usually Integers which represent process waiting to send 
data to Che buffer and process waiting to receive data from the buffer. 
These Integers are basic elements of process synchronization by the 
kernel which performs the blocking (delaying) or the activating of the 
process. 
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A typical Slagle buffer monitor Is shown In Fig. 3-18. (In this 
section the Pascal language. Ref. 14, Is used as a more precise and more 
comfortable means for program description than program flowcharts.) As 
seen, two monitor procedtires are associated with the monitor: SEND for 

transmitting data to the buffer and RECEIVE for taking the data from 
the buffer. Procedures BLOCK and CCNTINUE are kernel procedures which 
block and activate the processes In a pairwise manner. 

Resource monitors do not have buffers. They use multiprocess 
queues to schedule the resource among processes. This kind of monitor 
will not be used In CACS In Its usual form, and It will not be con- 
sidered here. 


5. Basic Prerequisites and Assumptions 

The Implementation of monitors and processes depends on com- 
puting capabilities assigned to the project. Other limiting factors 
are the attitude of the research personnel about the utilization of 
capabilities and directives derived from the experience during the 
foregoing work. In the case of the present project, the essential pre- 
requisites and assumptions can be summarized as follows: 

(a) Only one processor Is assigned to the CACS project. 

(b) The frfiole OPER subsystem should be Implemented as a single 
partition user program that Is Imbedded In the OS/16 

^f^2 real-time operating system. 

(c) The clocked I/O version for CACS software-to-manlpulator 
Interface is suggested. Ref. 12. 

The first prerequisite is due to current limited hardware 
resources. However, the eventual possibility of utilizing more than one 
processor should not be rejected. Current trends in microprocessor 
technology development make this perspective realistic, and specific 
suggestions have already been made. Ref. 12. Hence, software design, 
compatible with a multiprocessor environment. Is advocated. 

The second prerequisite is made for two reasons. First, the 
characteristics of the MT2 real-time operating system are not Ideally 
adjusted for this kind of real-time environment. It Is designed to 
fit more on-line Interactive Information systems than systems such as 
manipulator control. Second, the top-level operation of MI2 requires 
a fair amount of knowledge about its operation. Only with a single 
partition user program can the operator be unburdened from many details 
required by the operation of the MT2, Ref. 13. It should be mentioned 
that no suitable high level languages, especially for concurrent pro- 
gramming, are available at the present time for this machine. 


typ» vector : qrroy [l. .bufdlmj^f int«g«r; 

vor buffid: vector; | shored data. | 

vor full : booleon; {Buffer full*etnpty indicator.} 

vor sender, receiver: integer; | process queues*} 


procedure SEND (orgvec: vector); 
begin 

if full then BLOCK(sender) else 
begin 

I put orgvec into buffer } 
full:=true; 

CONTI NUE(receiver) 
end 
end; 


procedure RECEiVE (vor oi^vec : vector); 
begin 

if not full then BLOCK(receiver) else 
begin 

{ get content from the buffer into orgvec } 
full:*faise; 

CONTINUE(sender) 

end 

end; 


Fig. 3-18. General Outline of Buffer Monitors 
(in Pascal language) 
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Finally, the clocked I/O is the most suitable vehicle for 
communicating between the digital computer as the controller and the 
continuous dynamic system as the controlled object. 

6. Implementation of Monitors and Processes 

Taking into account the assumptions made in the preceeding 
subsection, it is obvious that the problems of mutual exclusion, proc- 
essor switching, and scheduling of other physical resources do not 
exist any more. Therefore, the implementation of monitors and proc- 
esses becomes trivial. The processes can be developed as siiq>le 
subroutines and their switching to the processor as simple subroutine 
calls. Also, the producer-consumer relations between processes are 
not asynchronous due to periodic I/O operations. It means that in 
every iteration cycle, all input buffers are first filled by an input 
process (DOIO), and then they are read by the internal processes. The 
same is true with the output buffers and the other internal buffers 
for shared data handling. In addition, most of the internal buffers 
are not of the producer-consumer type. As shown in Subsection 4.2, 
the condition for putting a new value into the buffer is not its empti- 
ness, and the condition for getting the content of the buffer out is 
not its fullness. This considerably simplifies the monitor procedures 
and their data structures, because the buffer full-empty indicator as 
well as the receiver and sender queues are not needed anymore. There- 
fore, the data structures in the buffer monitors are reduced to the 
buffer content itself. Consequently, the kernel is dramatically sim- 
plified and reduces to a single monitor of the synchronization type and 
to one short program. 

Now, the justification of the monitor concept in this project 
might be appropriate based on the following reasons: 

(a) The complexity of CACS software and its dynamic environment 
demand a highly structured organization. Concurrent proc- 
esses and the monitor concepts are powerful methods for 
such software structuring, even if the processes and monitor 
are simulated. The natural parallelism of manipulator con- 
trol activities reinforces this reason considerably. 

(b) The simplified versions of the monitors used out of a 
language context can dramatically reduce the hazard of 
destroying the system integrity. The system overhead 
(in time and space) , which can be considered a procedural 
drawback for data communication, is still negligible in 
comparison to other computations in the system. 

(c) Because the monitor calls and process switchings are 
reduced to simple subroutine calls, the kernel run-time 
overhead is practically reduced to zero, and therefore, 
cannot be used as an argtiment against the monitor concept. 
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(a) If the system is built based on the monitor concept it can 
be easily transferred to a multiprocessor environment. 

As mentioned earlier, multiprocessor control of a manip- 
ulator has a real future. 

(e) Even in the single processor case more complex scheduling 

strategies must be considered if the execution time of CACS 
programs exceeds the maximum al lovable value of the cycle 
time. Preemptive scheduling strategies and processor 
switching will not disturb the general system structure if 
it is based on the monitor concept. Only the monitor pro- 
cedures must be extended by queuing features, and a new 
version of the multiprocessor case. 

Let us now consider the implementation of the monitors in the 
OFER subsystem. The list of all OPER monitors is given in Table 3-7. 

All monitors listed, except the last three, are of the buffer 
type. However, the majority of the Input/output buffer operations, as 
required by the control algorithms, are not of the simple send/receive 
type. They are dependent on some specific conditions, or they have 
to provide some nontrivial data conversion. For example, calibration 
and decalibration (MPXS, MJIN, MJOU) , computation of trigonometric 
functions (MSCT), computation of matrix coefficients, and coordinate 
transformations (MSCH, MAT) . These operations could be done outside 
of the monitors, but for structural reasons they are kept inside as an 
integral part of the data transmission. A detailed specification and 
description of the monitor procedures is given in the program dociimen- 
tation (Appendix B). 

Monitors MilEL and MHLD are used as locking flags for controlling 
some arm manipulations, such as putting it in the hold state or releas- 
ing it from the hold state by operator intervention. These monitors 
and related locking mechanisms will be explained together with the 
processes that use them. Monitor MPRQ is Intended for process scheduling 
and synchronization. It is discussed in the next subsection. 

The implementation of Che monitors and processes will be done 
in an identical manner in accordance with the features of the INTERDATA 
Common Assembler Language (CAL). Thus, every monitor or process will 
be coded as separate assembly block called module . The structure of a 
module is given in Fig. 3-19. As seen, it contains three parts: list 

of entry-point and external symbols, program code (procedures), and 
data declarations. Programs within one module can be hierarchically 
structured, and they cannot be accessed by procedures defined in other 
modules unless they are listed in the entry-point list. The same is 
true with data. If their identifiers are not included in the entry- 
point list, they remain private variables of Che monitor and cannot 
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Table 3-7. CACS Monitrrs 





ENTRY-POINT AND 
EXTERNAL SYMBOL LIST 


PROCEDURES 


DATA DECLARATIONS 


3-19« General Structure of the Module 
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be accessed from other modules. The monitor procedures called by a 
process must be In the external list of the process. The same holds 
for nested monitor calls. If the monitor procedure calls a procedure 
from another monitor, the later must be in the external list of the 
former module. The entry-point and external symbol features (lENTRY and 
EXTRN pseudo-instructions) are the only means for access right control 
in the CACS . 

It is accepted here ttiat all shared variables are liLplemented as 
private variables of the modules, i.e., they must not appear in any 
entry/external list. Furthermore, all private variables of processes 
which must be initialized are implemented as global variables declared 
as entry-points. 

In order to describe the hierarchical structure of the system, 
the graphical technique called access graph has been used. Ref. 9. By 
access graph the monitors and processes are depicted as circles and the 
access right by arrows (the latter should not be Interpreted as data 
flow!). The access graph of the OPER subsystem is given in Fig. 3>20. 
In this figure all internal processes are represented by oae circle to 
simplify the diagram. 

A more precise definition of the access rights can be done by a 
cross-reference table of the form shown in Fig. 3-21. The columns of 
this table represent processes and monitors, while the rows represent 
the monitor procedures which are called by processes and/or monitors 
indicated at the head of the columns. The access right is checked by 
"X". For example, "X" in the first column and the second row of the 
table means that process Pi has access to procedure S 12 of the monitor 
Ml. As a consequence, the module representing the process Pi must have 
a pseudo- Instruction EXTRN S 12 , while the module representing monitor 
Ml must have a pseudo- instruction ENTRY Si2. 


7 . Scheduler 

Scheduling of I/O ooerations and Internal processes is provided 
by the program called scheduler (SCHED) and by one monitor (MPRQ) called 
process queue. These twc system components will be discussed in the 
present subsection. The scheduler is the main program of the subsystem 
OPER, whicn is a clock interrupt driven cyclic program. The Pascal 
program of the scheduler is given in Figs. 3-22 - 3'*23. 

I/O operations and cycle control is done by subroutine DOIO, 
which reads A/b converter rnd swltches/pushbuttons and writes D/A con- 
verter and alarm display. These operationa are performed through spec- 
ial clock interrupt service routines. After executing service routines, 
the control is given to the first Instruction following DOIO call. 

The rest ct the main progra.i is executed within a clock i.iterval. After 
execution, the processor is trapped in a “busy wait" loop (1 goto 1), 
until the next clock interrupt. 
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Fig. 3-21. Access Right Cross-Reference Table 












jiniHolizotion twirch.} 


\S*art clock, porfoim I/O oporatiem on 
ADC, OAC, twitch«/puihDuttont ond 
olom dliplo^. Rotum oddroas of Clock 
Infomipf Son^ico follo^Ml this stotomonU | 

iNITSW(tin); | Got sin from SWTBUF. { 

If not sin thon flrst:«truo olio 

bogin 

If first thon 
boflin 

first iso; 

INIT; jlnltiolizo OP8R.I 

ond 

MUPLEX; {Multipiox intomol procossos.| 


program SCHEO; 
vor first : beoioon, 
sin : booioon; 
lobol 1; 
bogin 

first :«truo; 
DOiO; 


I goto I 
ortd; 


]Busy wait until clock intorrupt,! 


Flg« 3"22< Scheduler - Mein Prograa 
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PfOc<duf MUPLEX; 
vor pindpx : intagtr, 
procno : integer, 
nc : integer; 


pindex:»0 

whiie pindex < piocno do 

^22i2. 

pindex:*pindex *■ 1; 
RESUME(pindex, nc); 

if nc ^ 0 then 


{Process ID index.} 

{Number of processes (porometer).} 

{Number of skipped dock cycles ond 
resume indicator. | 


{Check if process is for resumption 
(if nc ^ 0 ).{ 


begin {Branch to process with index pindex} end 


Fig. 3-23. Scheduler - Process Multiplexing 
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In every iteration cycle the switch SIN is interrogated. If SIN 
is set, the subroutine INIT is called. This subroutine performs ini- 
tialization of all pertinent shared variables (for example: process 

queue, alarm status table, event status table, etc.). Initializations 
are done by executing corresponding monitor procedures. It is obvious 
that Initialization of the system is possible at any moment without 
terminating the program SCHED. This feature is helpful for experimental 
work with a manipulator. 

Scheduling of internal processes is accomplished by subroutine 
MUPLEX, which is shown in Fig. 3-23. This subroutine scans during every 
iteration cycle all internal processes represented by their index (pin- 
dex), and it examines if they are for execution or not. This examina- 
tion is done b> subroutine RESUME which is a monitor procedure of the 
synchronization monitor, MPRQ. The scanned internal process will be 
resumed, l.e., the corresponding subroutine will be called if the out- 
put variable nc of Che RESUME is nonnegative. If it is negative, the 
next process will be examined . If it is positive, the variable nc 
represents the number of skipped iteration cycles in the case of delay- 
ing the process. This number is important Information for the process 
which must know Che exact time passed since the last execution in order 
to provide time dependent computations. The variable nc is set to a 
negative value if the process is blocked or delayed. 

Figure 3-24 shows the monitor MPRQ which is the basic synchroni- 
zation component of subsystem OPER. The data structure of the. monitor 
consists of five vectors: pqueue, pqlnit, permpr, delays, and dcount. 

Vector pqueue is a process queue consisting of a bit map pointing to all 
processes waiting for execution. The process identification is pro- 
fi .ed by a vector index, so that pqueue (plndexj » true, means that the 
process indexed by pindex is ready for execution. The vector pqlnit is 
an Initial value of the pqueue, while the vector permpr is a bit map of 
the same size as pqueue and pqlnit, defining permanently active proc- 
esses. Permanently active processes cannot be delayed. Integer vector 
delays contain information about delay, i.e., the number of iteration 
cycles which have to be skipped by the scheduler. For example, delays 
[pindex] > 6 means that the process pindex will not be executed for the 
next six iteration cyules. During every iteration, this number is 
decremented by one until delays [pindex] becomes zero. The vector 
dcount is of the same size as delays, and it contains information about 
the number of iteration cycles being skipped since the last delay opera- 
tion. It is used to form the output variable nc of the procedure 
RESUME. 

Monitor MPRQ contains six procedures: INITPQ, INITDS, SIGNAL, 

WAIT, DELAY and RESUME. The first two procedures are called by sub- 
routine INIT, Fig 3-22, in order to initialize the vectors pqueue, 
delays, and dcount. Procedures SIGNAL and WAIT are used oy all internal 
processes for synchronization purposes. The first procedure makes the 
process active, while the second one blocks the process. As seen, it is 
ensured that permanent active processes cannot be blocked. These 
procedures are tools for implementing the process precedence chart 
given in Fig. 3-15. 


3-58 


oortst mctxpno - 32; | Moximal number of processes*} 

type integv : ofroy [l«. moxpno] integer; 

type boolve : orroy [l . • moxpno] boolean; 


vor pqueue : boolve; 
vor pqinit : boolve; 
vor permpr : boolve; 
vor delays : integv; 
vor dcount : integv; 

procedure INTPQ; 
begin 

pqueuer^init 

end; 

procedure INITOS; 
begin 

deloys:^; 

dcount:^ 

end; 

procedure SIGNAL (pindex : integer); 
begin 

pqueue (pindexj :^rue 
end; 


I Process queue*} 

{Initio! value of process queue*} 

{ Def inition of peimonent active processes* } 
{Delay semaphore*} 

{Oeloy count*} 

{ Initioiixe process queue*} 


{Initialize delay semaphore*} 


{Signal process pindex*} 


{Wait process pindex*} 


procedure WAIT (pindex : integer); 
bjg|n_ 

ifjnot^p*rmpr [pindex] th*n pqueue [pindex] :«faise 
end; 


Fig. 3-‘?4. Process Synchronization Monitor MPRO 
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proctdurt DELAY (pind«x^ nd : integer); jOeloy process pindex nd dock cycles* { 
begin 

if not permpr [pindex] then 
begin 

if deloys [pindex] < nd then delays [pindex] :^d 
end 


procedure RESUME (pindex : integer; vor nc: integer); 

{Examine if pioces pindex is for resumption*} 


begin 


nc:^ *1; ] Initializotion of skipped clock cycle number nc«| 

I (if nc < 0 process is not for resumption)*} 

if not pqueue [pindex] then 
begin 

delays [pindex] 0; 
dcount [pindex] 0 
end else 
begin 

i^ delays [pindex] > 0 then 
begin 

delays [pindex] delays [pindex] - 1; 
dcount [pindex] doount [pindex] + 1 
end else 
begin 

nc;- dcount [pindex]; 
dcount [pindex ] :- 0 
end 
end 
end; 


Fig. 3-24. Process Synchronization Monitor MPRQ (Continued) 
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The DELAY procedure causes blocking of process for a limited 
period of time, defined by the number of Iteration cycles to be skipped 
by the scheduler. If during the same Iteration cycle and within the 
same process two delay operations are called, the one with the greater 
delay will have priority. The permanently active processes also cannot 
be delayed. 

The procedure RESUME is called by the scheduler. Its purpose is 
to determine the variable nc based on the state of the vector pqueues 
and delays. It updates also the vector delays and dcount. 


8. Program Documentation 

The programs under design at the time of writing this report will 
be presented In Appendix B. To facilitate the program design phase and 
future program maintenance, the program documentation Is elaborated 
systematically by the use of forms which are specifically developed for 
this purpose. There are three types of forms with the following titles: 
"Module Definition,” "Procedure Definition," and "Data Definition." All 
of these forms have Identical upper right corners and bottom parts. The 
upper right comer contains the identifier of the object (module, pro- 
cedure, or data entity) . The bottom part Includes Identifiers of sys- 
tem, subsystem, designers, date, and page number for the case when the 
form is continued. This uniquely defines each form regardless of the 
time, system or project. For easy look-up, these forms can be sorted 
alphabetically by the object identifier, where different object types 
a.'e separated. Examples of compiled forms are given in Figs. 3-25 - 
3-27. Entries in the forms are self-explanatory through the corre- 
sponding headers, and the following comments will help their 
understanding . 

The procedure hierarchy of the module is given by the line 
indentation, rather than by a hierarchical diagram. Procedures are 
represented by their identifiers together with identifiers of their 
I/O parameters which are enclosed in parentheses. 

The dimensions of data arrays (number of data items in the case 
of records) are enclosed in square brackets. The procedure identifiers 
which appear in the entry-point list are underlined, while the identi- 
fiers of procedures which belong to another module, l.e., which appear 
in the external symbol list, are enclosed in the brackets. The same 
is true for data identifiers. Explanatory comments are given in nar- 
rative form on the right side. Every part of the form has a check-box 
which has to be checked by an "X" if that part of the form must be con- 
tinued. For this purpose, a new form must be attached to the first one 
and page numbers must be given. 

A procedure definition form is created as HIPO (Hierarchical 
Input-Process-Output) form. Column "I/O Method" defines the way of 


MODULE DEFINITION 


IDENTIFIER: 

AAJOUT 


FURroSE/OESCRIPTION: 

monitor, supports access to calibrated output joint variables. 

c 

PROCEDURE HIERARCHY (ENTRY-POINT SYMBOLS UNOERUNED, EXTERNAL SYMBOLS 
ENCLOSED IN BRACKETS): 

UPOTEJ 


[G£THL0(HLDF[1])] 

AAonitor MHLD 

DJARM(VPC3]) 

[GEnRPCVPCa] )] 

Monitor MTRP 

[GErTRR (VR[3])3 

Monitor MTRR 

e3ETSa(SICOTFll2l)] 

Monitor MSCTF 

DJHND(VR[3])) 

CGCTR0P(VP[3])] 

Monitor MROP 

[GETROR(VR[3])] 

Monitor MROR 

DJJAW(VF[1] )) 

[GETJAW(VJ[1])] 

Monitor MJAW 


CPJRAW (oraument[7]) 


DATA STRUCTURE (ENTRY-POINT SYMBOLS UNOERUNED. EXTERNAL SYMBOLS 
ENCLOSED IN BRACKETS): 


hldf[i] 

Hold floQ (HLD) 

VP [3] 

Auxilloiy voriablo 

VR[3] 

Auxiliaiy voriablo 

VJ[l3 

Auxiliary voriablo 

SIC0TF[12] 

Trigon. fun« of joint* var* 

JOUT[7] 

Colib* joint voTe (monitor) 

[Z0TETA[7]] 

Zoro offset! (poromotor) 

CsftetaC/]] 

Scaling factor (poromotor) 


□ 


SYSTEM: 

SUBSYSTEM: 

DESIGNER: 

DATE: 

PAGE; 

OF: 

CACS 

OPER 

M. V. 

3/15/79 

1 

1 



PROCEDURE DEFINITION 


IDENTIFIER: 

DJARM 


IDENTIFIER I DIM. TYPE 


DESCRIPTION 


I/O METHOD 


incr«m«nts of orrifjoinf variabiM 

A6^, Afg) 


od - link 2 


1. Tok« from bu/Fof TRP by GETTRP. 

2. Tokw^frenbufforTRRby GETTRR. 

3. Comp>w«> ^ 5 ^ » 

4. Tokm iln ( 6 jp) and ooi ( 5 Ip) (I = 1 , 2, 3) from bufrar SGT by GETSCT. 

5. Trarafenm ooordinatai from world to joint space: 

^A*(®1P ^2F' ®3F^* 


SYSTEM: 

SUBSYSTEM: 

DESIGNER: 

DATE: 

PAGE: 

OPr 

CACS 

OPER 

M, V. 

3/15/79 

1 

1 


Fig. 3-26 








































passing Che parameters co/from the procedure. In case of CAL Che 
following notation will be used for different methods: 


(a) ad “ link + m Argument passage by argument address 

located Immediately after calling Instruc- 
tion (m Is relative position from return 
address) . 


(b) ad « 

(c) 

(d) local 


Argument passage by argument address which 
Is contained In Che register R^. 

Direct argument passage, as a content of 
Che register R^. 

Argument passage through local variable of 
Che module. 


(e) 


No arguments are passed to/from the 
procedure. 


The description of procedure functions can be done In a narrative 
form or by using flowcharts, pseudo-codes or higher level languages 
(Pascal, for example). 

Data definition forms are Intended for simple data structures, 
such as arrays or simple records. 
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SECTION IV 


CONCLUSIONS AND PLANS 


The funcClonal and operational description of the CACS software 
is given in this quarterly progress report. The operational description 
comprises the basic description of CACS hardware, the structure of con- 
trol algorithms, and the principles of the system software implementation. 
In this document the project objectives and framework have been clarified, 
the method of approach has been established, and the design of basic 
system software components, monitors and basic control routines, have 
been completed. 

The next phase of development will consist of three activities: 
software implementation, additional control algorithms development and 
real-time experimentation. 

The software implementation will be performed in two steps: 
implementation of the monitors and implementation of the processes. 

Since the monitors are the basic system components, containing procedures 
which are common to all processes, they will be implemented separately, 
out of the top-down development line. Therefore, the set of special 
off-line testing programs will be elaborated. These programs will be 
Implemented in a high level language, in this case Fortran V. Their 
purpose is to test extensively all monitor procedures, in order to 
provide a reliable and secuie programming base for further process imple- 
mentation. This approach of program development is usually called 
"bottom-up development . " 

The implementation of the processes will be carried out by top- 
down step-by-step refinements. The development will start with the main 
program of OPER subsystem, the scheduler, and will continue with the 
development of the processes. First the permanent active processes (lOP, 
OES, OMA) will be implemented, and then other Internal processes, starting 
with MUC and ending with ASO, according to the natural sequence of manip- 
ulator operations. The processes will also be Implemented by a top-down 
technique, as far as it offers practical benefits. As a consequence of 
this approach, no special testing programs are needed, except the proce- 
dure stubs for simulation of as yet unimplemented lower l<'.vel procedures. 

The control algorithms will be developed within the frame of the 
proposed principles, and the basic control routines and other relevant 
monitor procedures designed in this report. It should be noted that 
almost all CACS algorithms have been analyzed and evaluated by a look- 
ahead design which has preceded this report. The final version of the 
algorithms and related parameters will be established by experimental 
procedures . 
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The real-Cloe' experloents jlll be carried ont using the full 
C4.paclty of the CACS hardware under re^<^ oporatlonal condiclona. The 
purpose of this activity Is to check out the system Interface and tc 
adjust and review algorithmic parameters, calibration constants, and 
other pertinent data. Tn addition, the system performance and limita- 
tions as well as the Impact of the real operational conditions on system 
functioning will oe examined systematically In oiaer to be able to 
Improve control performance through algorithmic modifications. 

The u«xr progress report will describe the development activities 
outlined above. The first part of the report will be devoted to testing 
problmns. The second part will consider processes, l«e., the related 
algorithms and their aoftwaje implementation. 


(^Kt.CCD!NG PAGE BLANl^ NOT FILMED 
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APPENDIX A 

GRAPHICAL REPRESENTATION OF DATA AND STANDARD OPERATION 


A-1 


In order co eimplify the syitem deecrlpclon, a graphical 
technique ia employed in thla report. This technique conaiata of 
graphical symbola which express both the data and the operations on 
Che data. 

Two categories of data will be considered here: scalar and 

vector data. The former represent single data values with which a 
symbolic name is associated and indicated on the diagram. Vectors 
represent one-dimensional data arrays, which are ordered secs of scalar 
data under one symbolic name. 

Operations on data can also be considered as two general cate- 
gories: elemetary algebraic operations and data buffering . Elementary 

algebraic operations include summation and multiplication of scalar 
and/or vector data, as well as elementary arrangements/rearrangemencs 
of data arrays (Join and disjoin operations). For operations wtiich 
are more general special symbols are used. This will require additional 
explanations. 

Data buffering are specific kinds of stor«ge operations which 
play an important role in CACS software. Data buffers are memory 
locations with dimensions corresponding to the dimensions of the data 
arrays indicated on the diagram. These locations are accessible through 
special procedures which ensure mutual exclusiveness and access right 
control. The length of all buffers will be the same. This means that 
all buffers can accept only one scalar/vector data entity. 









Algebraic tummoHen 


x*x - y 





store X In buffer ff It It ampty. If ft 
Is tiill, the store operotfon will be 
deloy^ until the ^ffer becomes empty. 









OUTPUT BUFFER OPERATION 
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ABSTRACT 


This study investigates a flow shop scheduling problem 
which is defined as the linear flow shop problem (LFP) . 

The time to perform the jobs in the shop is a linear 
function of the batch size. The results of the research 
show that the shortest processing time (SPT) sequencing is 
the optimal solution for the mean lateness, m-'?an flow time 
and waiting time. A typical scheduling problem that can 
fit the definition of the linear flow shop problem is an 
assembly line. 

The main objective of this research is to find an 
optimal solution for a production scheduling problem that 
will be easy to implement in practical situations. The 
results show that the optimal solution is very easy to 
implement in real life problems as the jobs have to be 
arranged according to a monotonic increasing processing 
times. 

The proofs for mean lateness and mean flow time use 
the idea of pairwise interchange of adjacent jobs 
in the sequence. The method of mathematical induction is 
used to prove that the theorems hold for the case of m jobs 
in the sequence. A numerical example is presented to 


i 





explain how the theory is implemented. A sensitivity test 
is conducted on the numerical results to show how the 
shortest processing time sequencing yields the optimal 
solutions. 
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INTRODUCTION TO SEQUENCING AND SCHEDULING 

1.1 Definition and Cla»$ification of > Sequancing/ 
Scheduling Problen 

A sequencing problem is defined as one that determines 
the relative order or sequence in which given jobs or tasks 
are to be performed by the available facilities or resources. 
A scheduling problem is defined as one that assigns the 
actual starting time for each task on various resources. 

! In general, there are two kinds of scheduling problems, 
frhe static problem handles a set of tasks that is available 
for scheduling before the scheduling process starts. In 
this case the set of jobs is fixed and does not change over | 
time. The dynamic problem deals with a scheduling process 

where in addition to a given fixed set of tasks, new jots | 

i 

arrive to the system after the scheduling process starts. i 

i 

[The dynamic problem presents a real life production process 

i 

where the set of tasks to be considered for scheduling 
contains two subsets of tasks. The first subset includes 
tasks taht the system has from orders on hand and the second 
subsystem includes tasks from orders that will arrive to the| 
system in the future. 
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1.2 Types of Scheduling Models 


To classify the major scheduling models it is necessar) 
to characterize the configuration of resources and the 
behavior of the tasks. Each model can be used in a static 
or a dynamic approach 

(a) Sxngle Machine: One machine is continuously 

available and is never kept idle while work is 
waiting 

(b) Parallel Machines: Several identical machines 

are available for scheduling the tasks. Tasks 
are performed starting at tin’.e t»0, consecutively, 
so that as soon as a task is completed another 
task can put on the machine that is freed. 

(c) Flow Shop: Jobs to be scheduled follow a fixed 

routing and the routing is the same for all jobs. 

(d) Job Shop: Routing for all jobs is fixed but each 

job may have a different routing. 

The solutions to the various models are achieved mainly 
by two kinds of techniques; optimal techniques and heuristic 
techniques. The optimal techniques use mathematical 
programming such as linear programming, integer programming, 
dynamic programming and enumeration techniques like branch 
and bound. Mathematical approaches become very complex for 
practical cases as computational requirements will be severe 
for large problems. Even for relatively small problems. 
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there is no guarantee that the solution can be obtained 
quickly. Heuristic techniques usually assign a priority, 
or set of priority rules to sequence the jobs that are 
ready to be performed. Computational problems are avoided 
and the solutions to large problems can be obtained with 
limited computational effort. The problem with heuristic 
methods is that they do not guarantee optimality. The goal 
in practical problems is to get a **good" (not optimal) 
solution by using a simple sequencing rule, in a relatively 
short time. In special cases the results of using heuristic 
priority rules will result in an optimal or near optimal 
solution. 

1.3 Performance Measurements of a Production Schedule 

Quantitative measures for evaluating schedules are very 
important to determine whether a specific schedule is effi- 
cient. Performance measurements are functions of variables 
that define a scheduling process. The basic variables are: 

Ready Time (r^). The point in time at which job i is 
available for procesi^ing. In a static model where 
all the jobs are available for processing at time 
zero, rj^«0. 

Processing Time (t^) . The amount of processing 
required by job i 

Due Date (d^) . The point in time at which the proces- 
sing of job i is due to be completed. 
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Conpletion Time (C^). The time at which the processing 
of job i is finished. 

The major quantitative measures for evaluation are: 

Flow Time (Fj^) , The amount of time job i spends in the! 
system: " *“1 

Fj^ • Cj if r^^ « 0 . 

Makespan (MS) . The amount of time from zero until all 
jobs are completed or the completion time of the last 
job. 

Lateness (L^) . The amount of time by which the comple-j 
tion time of job i exceeds its due date: L. » 

Tardiness (TR) . The lateness of job i if it fails to 
meet its due date, or zero otherwise: TR^ « max{ 0 ,L^}i 

Idle Time (ITj). The amount of time that machine j 
is not productive. 

Schedules are generally evaluated by aggregate 
quantities that involve information about all jobs , resulting 
in one dimensional performance measures. Suppose that m 
jobs are to be scheduled on n machines then aggregate 
performance measures are: 


Mean Flow Time; 


Maximum flow time: 


1 Z F. 
m ^ 1 


i*l 


F_„ = max {F. } 
max 1 
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f 



m 


Mean Lateness: 




Maximum Tardiness: ■ max{TR. } 

max .gj 1 


1.4 Research Objectives 

This research investigates a linear flow shop model 
where the time to perform the jobs on the machines is a 
linear function of the batch sizes. The results of the 
research show that the shortest processing time (SPT) 
sequencing is the optimal solution for the mean lateness, 
mean flow time and waiting time. 

The main objective of this investigation is to find 
an optimal solution for a production scheduling problem 
that will be easy to implement in practical situations. 

The results show that the optimal solution is very easy to 
implement in real life problem as the jobs have to be 
arranged according to a monotonic increasing processing 
times. The results can be implemented in processes as 
assembly lines where the jobs are following the same 
routing. 

The proof for the minimum mean lateness and the minimum 
mean flow time is carried out by using mathematical induc- 
tion for the case of m jobs in the system. An analysis of 
a numerical example shows how the results are implemented 
in a specific case. 


5 



CHAPTER 2 


FLOW SHOP SCHEDULING > BASIC CONCEPTS 


2.1 Definitions and Terminology 

A flow shop scheduling system is defined as a process 
where the jobs to be scheduled follow a fixed routing 
and the routing is the same for all jobs. 

The shop contains n different machines and each job 
consists of n operations, one for each machine as illu- 
strated in Figure 2.1 for a ^’Pure” n machine flow shop. 

The machines in a flow shop are numbered so that the j^^ 
operation of any job precedes its operation, then the 
machine required by the j^^ operation has a lower number 
than the machine required by the operation. The 
machines in a flow shop are numbered 1, 2,...,n and the 
operations of job i are numbered (i,l), (i,2) , . . . , (i,n) . 

It is not required that every job have an operation on each 
machine in the shop as in the case of a general flow shop 
illustrated in Figure 2.2. Jobs must not enter the shop on 
a single machine, or leave from a single machine. In a 
general flow shop, each job is treated as if it had exactly 
n operations, for in cases where fewer operations exist, 
the corresponding processing times are taken to be zero. 
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The only requirement is that all movements between machines 
within the shop be in uniform direction from machine j to 
j-*>l and from machine i*l to machine j*2 etc. 

An example of such a shop is an assembly line, where 
the workers or work stations represent the machines. 

However, a group of machines served by a unidirectional, 
noncyclic conveyor would be considered a flow shop. 

2.2 Characteristics of a Flow Shop 

A basic flow shop problem is characterized by these 
conditions . 

(a) A set of m multiple-operation jobs is available for 
processing at time zero. 

(b) n different machines are continuously available, with- 
out consideration of temporary unavailability for 
causes such as breakdown or maintenance. 

(c) Each operation can be performed by only one machine 
in the shop. 

(d) Setup times for the operations are sequence -independent 
and are included in processing times. 

(e) Job descriptors are known in advance. 

(f) Individual operations are not preemptable, once an 
operation is started on a machine, it must be completec 
before another operation can begin on that machine. 
Only a single interval (b,c) is to be assigned to each 
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operation with (c*b) equal to the processing time of 
the operation. 


! 


(g) Each machine can handle at most one operation at a 


time. Consider the interval (b^, c^) , the assignment 


of operation x to a particular machine. For every 


other assignment (b^, Cy) to that machine, either 


“•x - 'y " *=X i “y 


n 




CHAPTER 3 


REVIEW OF PAST WORK 


The literature on flow shop scheduling contains many 
papers that introduce general ideas, algorithms to solve 
specific problems and analysis of different performance 
measures. This literature survey presents the tyipcal 
papers in the area, papers that present general ideas or 
general techniques to solve the problem. 

The earliest optimal results for the two machine flow 
shop problem were obtained by Johnson (1954). The objectivr 
of Johnson was to minimize makespan or minimize the maximum 
flow time. Johnson proved that in an optimal sequence, 
job i precedes job j if: MIN{t^j,tj 2 ^ ~ MIN{t^ 2 »^jl^ 
where t^^ is the time to perform job i on machine 1 etc. 
Johnson extended this algorithm for the case of three 
machines. The problem loses some of the nice s'^ructure of 
the two machine case. The problem is formulated, however, 
and for the special cases where MIN{tj^^} ^ MAX{tj 2 ) or 
MINftij) > MAX{tj 2 ) the complete solution is found 
analogously to the two machine problem. 

Many papers have been written about other performance 
measurements but no one has found an optimal algorithm, 
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other than enumeration, for the general flow shop problem. 
Conway, Maxwell and Miller (1967) state in their book. 

Theory of Scheduling . "Even for the two machine flow shop 
the optimization of mean flow time is a very difficult 
problem... Johnson's procedure is not optimal with respect 
to this criterion and, in general, it is not even very good"i 
(p. 89). 

I 

Ignal andScharge (1965) applied a branch and bound 
technique to the three machine flow shop problem. They 
observe that the three machine maximiam flow time problem is 
easier for the "branch and Bound" procedure than the two 
pachine flow tims problem, in that a higher proportion of 
jthe job sets were solved with the minimum number of nodes. 

Jackson (1956) considers a case in which the m jobs 
■have a common machine for their first operation and a common 

Imachine for their last (third) operation, but in which the 

! 

jsecond operation of each job is performed on a different 

1 

machine. There are thus m^Z machines, m of them correspon 

I 

ding to a second machine in a flow shop which can process 
aany jobs simultaneously. The algorithm is similar to 

iJohnson's three machine method. I 

1 ■ 

Dudek and Teuton (1964) have proposed an algorithm to 
minimize the maximum flow time in a three machine flow shop, 
which is also applicable to larger shops if one arbitrarily 
limits consideration to permutation schedules. The algo- 
rithm suggests a method for selecting the job to be placed 

11 
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o 
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first in the sequence, the job to be selected from the 
remaining m>l to be placed second, etc. 

Wagner and Story (1963) have used integer programming 
to formulate and solve the problem of three machine flow 
shop, to minimize maximum flow time. 

Gupta and Dudek (1971) examine various optimization 
criteria and investigate the interaction of several co.*>t 
factors on optimal schedules. Based on the results of a 
sensitivity analysis performed to study the interactions of 
several cost factors on the optimal schedule, they suggest 
the adoption of minimization of total opportunity cost as 
the optimization criterion for the flow shop schedules. 

Akers (1956) has a graphical solution for a special floi 
shop problem. There are only two jobs to be scheduled 
through a flow shop of any number of machines. Suppose that 
the machine are numbered l,2,...,n in the order in which 
they process the jobs, so that the processing times are 
given by: 

^2,l*^2,2****»^2,n 

These times can be marked off on axes for the two jobs as 
shotm in Figure 3.1. A scliedule can be represented by any 
line. 
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1. From (0,0) to (2t, St, J; 

j«l •‘»J j«i 

2. That is composed of horizontal (work on job 1 only), 
vertical (work on job 2 only) and 45^ (work on both 
jobs) line segments; 

3. That does not enter the interior of any of the shaded 
regions (which would imply one machine working on both 
jobs simultaneously) . 

Line SL in Figure 3.1 is such a "schedule line." 
Mathematical approaches such as linear programming, 
integer programming, dynamic programming and branch and 
bound to the flow shop schedule become very coiq>lex for 
practical problems. The formulation of specific problems 
in mathematical models is time consuming and need experts 
that are not available in every production facility. In 
practical problems many variables are involved in the model 
and conq>utational requirements will be severe. Even for 
relatively small problems, there is no guarantee that the 
solution can be obtained quickly. In dynamic organizations 
where schedules change every day for example, solutions that 
are the output of a ^.mathematical model cannot be used 
efficiently as it takes some times hours to get the results, 
even if one uses a fast computer. There are cases where in 
the production plant there are not people who can read and 
understand the computer output and then implement it in the 
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shop. Using heuristic algorithms, computational problems 
are avoided. These techniques obtain solutions to large 
problems with limited computational effort. Experts are 
not needed to develop mathematical models and to understand 
the output of the models. Heuristic methods usually assign 
a priority, or set of priority rules to sequence the jobs 
that are ready to be performed in the shop. Priority rules 
as first coa» first serve, shortest processing time and 
minimum slack time are easy to inclement in the shop floor 
and it is very easy to run a production facility following 
these simple priority rules. In many cases priority rules 
are used in production plants as a result of experience of 
years and they turn to yield good results. The problem 
with heuristic methods is that they do not guarantee 
optimality. However, the goal in practical problems is to 
get a **good" (not optimal) solution by using a simple 
sequencing rule, in a relatively short time without the 
need of experts. In special cases the results of using 
heuristic priority rules will result in an optimal or near 
optimal solutions. Heuristic rules of assigning jobs in a 
scheduling process can be modified after using it in the 
shop and getting some results about the level of efficiency 
Palmer (1965) developed a heuristic algorithm for the 
makespan problem. The algorithm gives priority to jobs 
having the strongest tendency to progress from short tines 
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to long times in the sequence of operations. He proposes 
the calculation of a slope index for each job. 






-(n-3)ti 2- 


where : 

n is the number of machines in the shop, 
ti is the time to process job i on machine n. 

Then a permution schedule is constructed using the job 
ordering ^ $2 ^ ... ^ 
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CHAPTER 4 

FORMULATION OF THE PROBLEM 


The machine scheduling problem considered for this 
research is for a flow shop facility engaged in the batch 
production of parts. Each job consists of a known number 
of parts (n) . The parts are produced, in order, on a 
number of macnines (j). The goal is. to sequence all of 
the parts production jobs (i) so that the entire production 
schedule is completed in the minimum amount of time. The 
feature that makes the problem both attractive and practical 
is that the processing time on each machine is a linear 
function of the number of parts in that job. 

A common representation of machine time is by 


tj. - bjij 


where 


tij • the time of job i on machine j, 

bj ■ processing time per Part of machine j , 

n^ » the number of parts in job i. 


This problem is a typical example of an assembly line 
problem where one has to produce different types or models. 
In an automobile assembly for example, you will produce n^ 
cars of Model A, n£ cars of Model B, etc. 
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The value of bj is either computed from historical data 

or estimated from machine operati;>^ characteristics. 

Periodically the value of is updated in order to keep it 

as a good representative of the real machine performance. 

The conditions that characterize a linear flow shop 

problem (LFP) are: 

a. A set of m multiple -operation jobs is available for 
processing at time zero. 

b. Each job requires n operations and each operation 
requires a different machine. 

c. Set up times for the operations are sequence independent 
and are included in processing times. 

d. Job descriptors are known in advance. 

e. n different machines are continuously available. 

f. Individual operations are not preeniptable. 





Substituting equation (4.2) in equation (4.1) yields 
the following linear relations between the time to perform 
a batch on a specific machine and the times to perform the 
same batch on all other machines. 

^ij * j,k t J, i e I (4.3) 

4.2 Time to Perform All Jobs Through All Machines 


Let Tj^ be the time to perform job i (i e I) over all 

n machines that the shop contains. 

n 

T; - (4.4) 

j»i J 

Substituting equation (4.3) into (4.4) yields: 
n 
T^ 




^ik ^ 


-r- i e I, k e J (4.5) 

°k j»l J 


Let T be the total performance time of all m jobs over 
all n machine‘s, then: 


m 

T - STj 
i«l ^ 


in in n ^ 

I)t. • E[-^ SbJ , k 
i»l ^ i-lL'^k j»l 


e J 


m 


m n 


t., 


T • Stj ■ E E ^ 

i-I ^ i-1 j-1 “k 


■BT ‘’j’ 


k e J 


(4.6) 


(4.7) 


(4.8) 
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Let Mj be the time that machine j is busy, then: 






4.4 Idle Time Involved in Linear Flow Shot 


Idle time is defined as the time in which machines or 
work centers are not producing any parts. In any productio 
plant the purpose is to minimize the idle time of machines 
in order to increase the total efficiency of the process. 

In case this idle time of the machines can be controlled 
and it is known in advance when the machines are not busy, 
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it is possible to take advantage of this fact and do some 
preventive maintenance or any other planned maintenance to 
keep these machines running. Let's define as the idle 
time of machine j, immediately before the i^^ job starts on 
machine j. There is no idle time on the first machine. 


4.4.1 Two Machine LFP 

Figure 4.1 illustrates the case of a two machine 
linear flow shop. In developing the terms for the idle 
time of machines in the shop, the linear relations 
constituted by the definition of the LFP problem are used. 
The idle time immediately before job 1 starts on machine 2 


^12 ■ hi 


(4.20) 


IZ2 - max{tj3^ ♦ tzi - t^^ - 


(4.21) 


The sum of (4.20) and (4.21) yields 


^12 * ^22 * ^21 ‘ 5^ ^ll'^ll^ (4.22) 

3 2 2 

I,, ■ max(St., • St., - Sli?,0} (4.23) 

il i.l 11 »>i i«i 11 iz 


The sum of (4.22) and (4.23) yields: 







3 

Et, 


hz hz * hz • 


3 3 

■ Sli 7 “ raax{ 23 1 
i-1 i-1 

b. 


Et,,, E 

i.i 11 i«i 


^il 


^11* ^11^ 


(4.24) 


The idle time immediately before the last job starts on 
machine 2 is: 


ni^ b l in* 1 

IjnZ * ®a3c(Et., - E til ■ S 

mz 11 bj 11 i2 


(4.25) 


The total idle time on machine 2 is: 


m 

m 


m-1 

m-1 

u 

m-2 

^hz 

i-1 

= maxlEt., 
i-1 

■ El 

E t-i, 
i-1 

E tjj 

i-1 

®2 

■ET 

E 

i-1 


til, 


.... E 


i.i ii ' B7 hi- hi* 


(4.26) 


In general 


m 


E 1^7* K 

i-1 

l<u<m 


(4.27) 


where 


u 


Et. 


“ “i 


'’2 y ' . 
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4.4.2 Three Machine LFP 


Figure 4.2 illustrates the use of a three machine 
linear flow shop. The idle time immediately before job 1 
starts on machine 3 is: 


11 * *12 *[l * 5j]*ll ■ S *1J 


*13 ■ * 


lot • max{ ♦ 2 1., - I 


(4.28) 


23 


i-1 ^ M3 - FJ ^11 


t,,.0} (4.29) 


The sum of (4.28) and (4.29) yields: 


^3 * ^23 • - TT til. (4.30) 


i-1 


i2 


i-1 i2 ■ ^ jt-i Ij 


j j 5 ^ ^ 

^33 * max(2t42 * " 2^i3.0) (4.31) 

^ i-1 i-1 *>1 i«i 11 i.i 13 


The sum of (4.30) and (4.31) yields: 


2 lit • max{2ti2 ♦ 2li2 “ F^ ^t.,, 2lit] 
i-1 i-1 i-1 *>1 i»i 11 fTi 13 


max{ 


Et., ♦ Ei,, - h E*<i. St,, * 


i-1 


•i2 ' jTfi* ft‘i2 

2 


Sli2 ■ 5^ tjj. Etj,) 

i»i 12 *>1 11 j.i iJ 


(4.32) 


The idle time immediately before the last job starts on 
machine 3 is: 
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n-1 


J.J ■ ^ S*i2 - 2 ‘ll - IlJ.OU^-SS) 


The total idle time on machine 3 is: 


m mm. m^ m-1 

2 1^3 - max{Et^2 * ^1^2 ' ^ ^il* ^ 

i-1 i-1 i-1 °1 i-1 i-1 


^i2 ^ 


m-1 


mr2 


^i2 * TT” ^ii » • * * • 2 * 2 " Tl^^n » 

i-1 °1 i-1 i-1 i-1 ®1 


St j) 

j-1 


m-1 


m-1 


max{^ S^il " ^ ^il 23li2» ^ ^il " 

h fTi ^ i-1 W hj il 

V m-2 ©•! w 2 • 

3 y . . y T ^ Tt 3 1 ^ 

i-1 i-i i-1 ®1 


S^i 2 » S^lj^ 


i-1 j-1 


(4.34) 


In general 


m 

E >12 

i-1 


max k. 


(4.35) 


1 < V < m 


where 


v-1 


K • ^ Ht,, i: t., . Zh, 

^ ^ i-1 n i-1 >1 
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4.4.3 The General Caf 


The general term for the idle time is: 

I . ^ z , . i:, 

13 ®i 1-1 ®i t-i i«i ‘j'l 

i-i 

S '‘r“» 

The total idle time on machine j is: 


(4.36) 


I. 


max{-^ ^ ^il ■*“ £li 4 1 ^ 

^ i-1 ^ i-1 61 ^ 

V T' ^ **i-l 

ITl ■ ®1 W * i.l 

2 2 W 

S'll * S'l*’ 


(4.37) 


4.5 Makcspan (MS) in Linear Flow Shop 

I Makespan is defined as machining time of machine n 

(N^) plus the idle time that is generated on machine n 

(lin)« Using equations (4.14) and (4.37) yields: 
m 


MS . M„ ♦ • 

m-1 


m-2 


"n 52 , . '’”•3 t I! t 

®Ii.l il itl'i."-!’ "ET i4i ‘il ■ SJ *n 

V I *’n-l Vt •’n , 

ft “HT i.l “ ' ^ *ll*ft‘i.»-l’ 


n-1 


i:*u> 


t-i 


(4.33) 
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f* ( 

^ h 


c 




iJ 


Substituting equation (4.2) into (4.38) yields: 


MS - b 


m n n 

2n. ♦ max(b , £n. - b T, n. ♦ 2l. , , 

n ii .1 1 n 1 1 n i i,n-i 

«-l m-2 n-1 

■1 ft "i ‘ ‘n i?i "i * S 


n-1 


b« 1 £n. - b_n, ♦ „ i» 2 t, ,} (4.39) 

n-1 ^ i n 1 ^ i,n-l* U 


4.6 Bounds for Perfornance Tine in a Linear Flow Shop 


4.6.1 Lower Bound (L) 

The lower bound is equal to the total performance time 
of all the m jobs over all n machines, or equal to the total 
machining time of all n machines. Using equations (4.11) 
or (4.19) yields: 


Ml 

L • E (n^ 


i-1 


n 

j-1 


j 


Z^(b. 2^n.] 

j-1 J i-1 ^ 


(4.40) 


4.6.2 Upper Bound (U) 

The upper bound is the lower bound plus the idle 
time that is generated in the process. Using equations 
(4.37) and (4<40) yields: 
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n ffi m n n ^ 

U » L + 2( £li4) • SUi 2b.) ♦2(roax(4^ 

• j-2 i-1 i»l ^ j«l ^ j-2 ”1 
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CHAPTER S 

OPTIMALITY IN LINEAR FLOW SHOP SCHEDULING 

This chapter describes three performance measurements 
of the linear flow shop scheduling: (a) Mean lateness, 

(b) Mean flow time, 

(c) Waiting time. 

It is proved that the shortest processing time sequencing 
yields the optimal sequence that minimizes the given 
performance measurements. A numerical example is given to 
ishow how the shortest processing time sequencing yields the 
jminimum mean lateness. 

I 

5 

S.l Lateness in a Linear Flow Shop 

Recall that job lateness is defined as L^ = C^ - d^, 
or the discrepancy between the due date of a job and its 
completion time. The objective of every production plant is 
to minimize the lateness in completing the orders because 
it is involved with penalty. 

Theorem 1 

The mean lateness in a two machine flow shop problem 
where the job times on the machines are a linear function 

i 

J 




i 

i 

h 

i 
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cf the batch sizes (t.^ ■ t. minimized by 

shortest processing time (SPT) sequencing. 

Proof 

The proof has two steps. 

(a) Consider a sequence S that is not a SPT sequence. That 
is, somewhere in S there must exist at least one pair 
of adjacent jobs, i and i***!, with l+l following i, 
such that t^ > Now construct a new sequence, 

S', in which jobs i and i+1 are interchanged in 
sequence and all other jobs are not changed. The 
situation is depicted in Figure S.l, where T^ denotes 
the point in time at which job i begins in sequence S 
and at which job i+1 begins in Sequence S'. A, 
denotes the set of jobs that precede jobs i and i+1 
in both schedules and B, denotes the set of jobs that 
follow i and i+1 in both schedules. 

The processing times date d^^ are 

given and completion time and the lateness of the job 
i are computed. 

t ^2 * processing time of job i on machine 1, compute< 
by 

ti 2 * processing time of job i on machine 2, computed 
by tjj - bjn^ 
d^ * due date of job i 
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(b) 


Figure 5.1 A Pairwise Interchange of 
Adjacent Jobs 

(a) Sequence S 

(b) Sequence S' 







Ci 


completion time of job i 
« lateness of job i» computed by ^ 

C. « completion time of the last job in set A. 

A in 


t-. 


The proof is sho%m for 2^ L- (sum of the lateness over all 
the jobs), as the mean lateness T « 


Using equation (4.1) the processing times t^^ and t ^2 
are computed by: 


hi “ hh 


t-2 « b2n. 


(5.1) 


(5.2) 


From equation (5.1) 


« 

“i TJ" 


Substituting equation (5.3) into (5.2) yields: 


t., = 


(5.3) 


i2 FT il 


Through all the proof steps it is considered that: 


b. 


(5.4) 


and 




*i2 • 'X^il 


(S.5) 


34 



The total lateness is: 


m 

* La * Li ^ L. Ln (5*6) 

k«l K A 1 fl 

where: is the total lateness of the jobs in set A. 

Lj^ is the lateness of job i 

is the lateness of job i+1 
Lg is the total lateness of the jobs in set B. 

The completion times of jobs i and i-»'l in sequence S 
are as follows: 

Cj^(S) « max{T^ + ^i,l»^A^ * ^i,2 “ * ^i,l» ^A^ 

♦ bt. 1 (S.7) 

Cj*i(S) - miut{Tj^ ♦ Cj(S)) ♦ 

“ * *i,l * *’*i*l,l C5.8) 


Using equation (5.6): 


m 

El. (S) 
k«l ^ 


La + maxiT^ + t. C^} + bt. ^ - d- + 

maxCTA + t^ ^ + fi+i^l»CiCS) }'»-bt^^j ^ 


‘^i+l * 


(5.9) 


The completion times of jobs i-i-l and i in sequence S' are 
as follows: 
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C.^j(S-) . max{T;^ ♦ j, CJ^^ ♦ j (5.10) 

CjCS') . «ax(T^ ♦ Cj^j (S-)) + btj 1 

(S.ll) 

Using equation (5*6): 


Zw(S') 

k»l ^ 


d. + LgCS*) 

Let's define P and Q as: 


La ♦ maxdA ♦ Ca) ♦ bt^^j ^ - d.^j * 

"“"a * ‘i»l,l ♦ ‘i.l- Ci.l(S')) ♦ bt^ j - 

di ♦ Lj(S') (5.12) 


* *1,1- V 

(5.13) 

* *i*l,l’ V 

(5.14) 


Then 

P > Q as t._j > ^ 

Substituting equation (5.13) into (5.7) yields: 

C.(S) • P + bt.^^ 

Substituting equation (5.14) into (5.10) yields 

Ci*i(S') - Q * bt_ , 


(5.14) 


(5.15) 








Then 


(S) > P ^ Q end ht^ 2 ^ i* ^ ^ ^ 

Let's define R and V as: 

R • max{T^ + t. ^ + ^i+i,i» C^CS)} ( 5 . 16 ) 

V . max{Tj^ ♦ j * tj j. q^iCS-)) (5.17) 

It is clear that R ^ V as C^CS) > C^+jCS*) 

Substituting equations ( 5 . 13 ) and ( 5 . 16 ) into ( 5 . 9 ) yields: 

m 

£lj^(S) « + P ♦ bt^ 2 - d^ + R + bt^+j 2 - d .^2 j 

+ LgCS) (5.18) 

Substituting equations ( 5 . 14 ) and ( 5 . 17 ) into ( 5 . 12 ) yields: 
m 

Slj^(S') « ♦ Q + bt2+2,l ■ ‘^i+l * ^ * ^^i,l ■ ^i 

+ Lg(S') (5.19) 


The result of comparing equations ( 5 . 18 ) and ( 5 . 19 ) is: 


m m 

Zwcs) - LlvCS') 

k«l * k -1 


^ * t’^i,! - <ii * R * bVi,l 


■‘^i+l * Q ■ **^i+l,l 


+ d 


i +1 


V * bt. T + di 

1 y 1 1 


Lb(S') 


P ♦ R - Q - V ♦ Lj(S) - Lb(S') 


'B' 


( 5 . 20 ) 



Equation (5.20) can be divided into two sums as follows: 
m m 

Sl^(S) - - (P ♦ R - Q - V) ♦ (LjCS) - 

LgCS')! (5.21) 

The first sum (P+R-Q-V) represents the difference in 
lateness between sequence S and S* for the set of jobs 
which are in set A and jobs i and i-*>l. The second sum 
[L£(S) • Lg(S*)] represents the difference in lateness 
between sequence S and S' for the jobs that are in set B. 
The first sum is independent of the second sum but the 
second sum is dependent on the first sum. 

Observing the first sum shows that the total lateness 
of sequence S is greater or equal to the total lateness of 
sequence S as P ^ Q and R ^ V, it means that the first 
job in Set B in sequence S will start at the same time 
or later than the same job in sequence S'. Concludes that 
Lg(S) > L^(S') as the jobs in set B are the same for 
sequence S and S' but the starting point is different. 

It is shown that 

m m 

Sw(S) > Sli.(S') as P > Q, R > V, L»(S) > U(S') 
k«l ^ ~ k-1 ” n - n 

and P > 0, R> 0, Q>0, V^O. 
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(b) In part (a) of the proof it is shown that any inter- 
change of a pair of adjacent jobs where a shorter job 




is following a longer job, can improve the total 
lateness or leave it unchanged. 

The second part of the proof shows that in the last 
interchange of a pair of adjacent jobs (going from a 
non-SPT sequence to a SPT sequence) a strict improve- 
ment can be made in the total lateness, concluding 
that a SPT sequence produces the minimum total late- 
ness. Part b of the proof has two phases. The first 
phase uses contradiction. Let us define a sequence S 
as a finite set of non-negative numbers such as 

^i,l' ^i+l,l» ^i+2,1- 

Sequence S is a non-SPT sequence in such a way that 
there is somewhere in the sequence only one pair of 
adjacent jobs i and i+1, with i+1 following i, such 
that tj ^ j^. Assume that sequence S is an 

"optimal" sequence. Construct a new sequence S*, in 
which jobs i and i-*-! are interchanged in sequence 
(so that you get a SPT sequence) and show that a 
strict improvement can be made in this "optimal" 
sequence. Therefore, the conclusion is (based on part 
a of the proof) that it is impossible for a non-SPT 
sequence to be optimal. 


s 
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The second phase involves a constructive proof that 
uses the method of mathematical induction to show that 
a SPT sequence minimizes the mean lateness. 

In this part of the proof the jobs which precede the 
pair of jobs i and i+1 that are being interchanged, 
are not considered, as the total lateness of the jobs 
in sequence S is the same as the total lateness of 
the job in sequence S'. The proof is shown for the 
case of two jobs where only jobs i and i-*-! are in 
the process, for three jobs; i, i**-!, i+2 and for four 
jobs; i, i+ 1 , i*»- 2 , i+ 3 . It is assumed it is true for 
m jobs i, i+1, i+2, . . . ; then it is proved it is 

true for m+1 jobs. 

Two Jobs Sequence 

The situation is depicted in Figure 5.2 where and 
C2 are the completion time of jobs 1 and 2 respectively. 

The terms for and C2 in sequence S are: 

€2(8) « t2j^ ♦ t22 ■ t2]^ + bt2j^ ■ (l+b)t2j^ ( 5 . 22 ) 

Cj(S) • max{C2» t2j + t^^} + t^2 • maxi (l+bjtj^* 

t2j^ + **’ ( 5 . 23 ) 

Using equations ( 5 . 22 ) and ( 5 . 23 ), the lateness is: 

l2(S) • (l*b)t2x - d2 ( 5 . 24 ) 


Lx(S) ■ maxi (l+b)t2x, t2i ♦ t^^} + bt^j - dj 


( 5 . 25 ) 








MACHINE 1 


MACHINE 2 



MACHINE 1 


MACHINE 2 


Figure 5.2. A Two Job Linear Flow Shop 

(a) Non-SPT Sequence 

(b) SPT Sequence 
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Using equations (5.24) and (5.25), the total lateness of 
sequence S is: 



(l^b)t2j • d2 ^ uiax( (l^b) t22 » ^21 ^ ^11^ 

♦ bt^i “ dj (5.26) 


The conqiletion times in sequence S' are: 

CjCS') ■ ^ ^12 * ^11 ^ ^^ll " (l^b)tj2 (5.27) 

C2(S') - max(C^, t^^ ♦ t2i) ♦ t2£ ■ max{ (1+b) t^j , 

tj2 ^ ^21^ ^ ^^21 (5.28) 

Using equations (5.27) and (5.28), the lateness is: 

L2 (S') ■ Cj - dj • (l+b)t2j - dj (5.29) 

L2(S') • C2 - 62 • max{ (l*b)t22, tj^ ♦ t2j) ♦ bt2i - d2 

(5.30) 

Using equations (5.29) and (5.30), the total lateness 
of sequence S' is: 

2 

S Lj^(S') • (l+b)t22 ’ **■ »ax{ (l^b)tj^j, tjj ♦ ^21^* 

bt2i “ ^2 (5.31) 
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Conq»aring equations (5.26) and (5.31) yields: 


2 2 

(Sj • (S ' ) • (l+b)t 22 ' ^2 ^ roax( (l^b)t 2 j » ^21 

♦ tjj) ♦ btj^ - dj - (l+b)t2i * dj - max{ (l<»‘b)t22, 

^11 ^ ^ 21 ^ " ***21 * **2 

■ t 2 i - ♦ max( (l+b)t 22 » t 2 ^ ♦ t^j) - max{ (l^bjt^i* 


Observing the terms max{ (l<«-b)t22 • ^2^ and 

max( (l'<-b)t22» + t2i) and given that t2j > 

(by definition) it follows that max{ (l‘»>b)t2i, t2x * tjj^) 
> max{ (l+b)tj2, tj^j ♦ t--} as (l+b)t2j > (I'^b)!^^ and 

(l+b)t2i ^ ^21 * *11’ ** ~ ^ 

2 2 

it follows that 2 l.(S) > JZl.(S') 

i -1 ^ ^ 

It is shown that the SPT sequence gives the minimuin 
mean lateness . 


Three Jobs Setiuence 

The situation is depicted in Figure 5.3. The compXetioi 
I times of the jobs in sequence S are: 


C £( S ) ■ t 22 ♦ t 22 “ ^21 * ***21 * (^^*^^^21 




(b) 


Figure 5.3. A Three Job Linear Flow Shop 

(a) Non-SPT Sequence 

(b) SPT Sequence 




Cj^(S) • max{C 2 , ^2l * ^11^ * ^12 * max{ (l+b)t 2 i» t 2 i * 

til} + btii (5.33) 

Cj(S) ■ max{Ci, t 2 i + tn + tji) + tj 2 * max{max[ (1+b) 

^ 21 * ^21 ^ ^21 ”*■ ^11 "*■ ^ 31 ^ 

btji (5.34) 

Using equations (5.32), (5.33) and (5.34) the lateness is: 

L 2 (S) * C 2 - d 2 (5.35) 

Li(S) • Cl - di (5.36) 

13 ( 8 ) = C 3 - dj (5.37) 

Using equations (5.35), (5.36) and (5.37) the total late- 
ness of sequence S is: 

3 

^ Li(S) * ~ ^2 * max{ (l+b)t2i, t2i + tn} + 

i“l 

btii ■ ^i niax{max[ (l+b)t 2 i, tii + tn] 
btii , ^21 ^ll ^ ^31^ ^ ^^31 ” ^3 (5.38) 

When observing equation (5.38): 
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max{(l+b)t 2 i, t 2 j + ■ (l+b)t 2 j^ as t 2 j > and 


b > 1 

(5.39) 

Substituting equation (S.39) into (5.38) yields: 


3 



L L. 

i -1 ^ 

(S) » (l+b)t 2 i - ^2 ♦ (l+b)t 2 i ♦ bt^j ' dj ♦ 



max{(l+b)t 22 ♦ ^t^^, * t^^ + ^ 32 )+ 



^^31 " ^3 

(5.40) 

The completion times of the jobs in sequence S' are: 


C^CS') 

' ^11 * ^12 • ^11 ^ ^*^11 * (l+b)t 2 i 

(S.41) 

CjCS'J 

» max{Cj, tj^i + t 2 X> + ^22 * max{ (1+b) t^^j. 



til + t 2 i) + bt 2 i 

(5.42) 

CjCS') 

= max{C 2 » t^j^ + ^21 ^ 3 l^ ^32 “ n»ax{max[ ( 1 +b) 


^ 11 * ^11 * ^ 21 ^ ^^ 21 ’ ^11 ^21 * ^31^ 



+ bt32 

(5.43) 

Using equations (5.41) , (5.42) and (5.43) the lateness 

is: 


Lj(S') » C, - dj 

(5.44) 


L 2 (S') . Cj - dj 

(5.45) 


LjCS') - C 3 - dj 

(5.46) 
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Using equations (S.44), (5.4S) and (5.46) the total lateness 
of sequence S' is: 


Sl^CS') ■ (l+b)tj^j^ “ ^1 ™a3c{ (l+b)t^j^, tj^j + t2j^}+ 

btzi - d2 + max{max[(l+b)tj^j^, tj^j^ + t2j^l + 

^11 * ^21 ^31^ ^ ^^31 ‘ ^3 (5*47) 

Comparing equations (S. 40) and (5.47) yields: 

3 3 

2^(5) - SLi(S') » (l+b)t2i - d2 + d+b)t2i ♦ bt^^ - 
+ max{ (l+b)t2j^ + ^^11 » ^21 ^11 * ^31^ 

bt3i - dj - (l+b)tjj^ ♦ d^ - max{ (l+b)tj^j^, t^^^ 
+t2i) + bt2j + ^2 - max{max[ (l+b)tj^j^, tj^j^ + 
^^21* ^11 * ^21 * ^31^ " ^f-11 


31 “3 


^21 ”‘ax{ (l+b)t2j^ + ^^11» ^21 ■*■ 

‘ll * *31> • *11 ■ ">ax{(l*b)tii, tjj ♦ - 

max{max((l+b)tjj^, t^^j^ + t 22 l + ^^21* ^11 * 


^21 * ^31^ 


(5.48) 


jlet'^ define X, Y, Z as follows 


X « max{ (l+b)tj^j^, tj^j^ + t 2 j^} 


4 


Y » max{X + bt2i* * t 2 i * (5. SO) 

Z « max{ (l+b)t2i + bt^^, t2^ + + tj^} (S.51) 

Substituting equations (5.49) , (5.50) and (5.51) into (5.48) 

yields : 

3 3 

S.Li(S) - XrfLj^(S') » t2i + (l+b)t2i + Z - tjj^ - X - Y 

(5.52) 

It is clear that Y > X (5.53) 

The relations between Y and Z are: 


Y « max{X + bt2j» + t£l ^31^ 

Z » max{(l+b)t2i + bt^j^, * tjj^} 

Then compare X + bt2j^ and (l+b)t2j^ + ^^11 after reducing 
the terms compare X and t2j + Two cases have to be 

checked: 

I. X = (l+b)tjj 

It follows that (l+b)tj^j^ < t2^ + btj^j^ as tj^j^ < t2j^ 

II. X » tj^j^ ♦ t2i 

It follows that tj^j^ + t2i £ t 2 i + btj^j^ as 
^11 ^ ^21 » b ^ 1 

It is found that Z ^ Y (5.54) 

3 3 

Sl.(S) - 5Il. (S') = t,, + (l+b)t,i - t.. + Z - X-Y 
i-1 ^ i=l ^ ^ 
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Using equation (S.54) 


Two cases have to be checked for X 

I. X • (l+b)t^^ 

3 3 

2l.(S) - Ll. (S') • t,i + (l+b)t2i - t,i ♦ 2 - Y - 
i»l i«l 

(l+b)tjl 

3 3 

It follows that Sl^(S) > Sl. (S') as t,, < t^, , 

i-1 ^ i»l ^ ^ 

2 > Y 

II. X » tjj ♦ t2i 

3 3 

Sl.(S) - El. (S') = t2i + (l-''b)t2i - t,i + 2 - Y - 
i«l ^ i=l ^ ^ 

^11 ■ ^21 
3 3 

It follows that El. (S) > El. (S') as t,, < t-,, 

i-1 ^ i»l ^ 

2 > Y, b > 1 

It is shown that the SPT sequence gives the minimum 
mean lateness. 

Four Job Sequence 

The case is illustrated in Figure 5.4. The completion 
times of the jobs in sequence S are: 
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^2liy.^31 I ^41 1 


•22 


El 


•32 I ^42 I 


t ! 

2 n ^-3 ^4 

(a) 

t ^ 1 ^ <• T i j 


EIT 5 I 


I hr^21*^31*"41 


jil_ 


MACHINE 1 


MACHINE 2 


MACHINE 1 


Til I ^22 




4 4 


Cl C2 C3 


±L 


(b) 


MACHINE 2 


Figure 5.4. A Four Job Linear Flow Shop 

(a) Non-SPT Sequence 

(b) SPT Sequence 




CjCS) 

“ hi * ^22 ■ ^21^^*^) 

( 5 . 55 ) 

Cj(S) 

• max{C2, t2i ♦ ♦ t^2 



■ max{ (l+b)t2j, t2^ ♦ 

♦ btj2 


■ (l+b)t22 ♦ *s 

^21 ^ ^ 11 * ^ ^ ( 5 . 56 ) 

Cj(S) 

• maxiCj, tj2 ♦ t2i ♦ tj^} 

♦ btji 


■ max{(l+b)t2i ♦ bt^^* 

* ^21 ^Sl^”*" ^*31 


- Z * btjj 

( 5 . 57 ) 

C4(S) 

• nax(Cj. tji ♦ tjj ♦ tji 

♦ t4j} ♦ 1.4^ 


■ max{Z ♦ btjj* ^21 



Using equations (5. 55) through (5.58) the lateness is: 


LzCS) » C 2 - d 2 (5.59) 
LjO) “ Cj - dj (5.60) 
Lj(S) « C 3 - dj (5.61) 
14 ( 8 ) - C 4 - d 4 (5.62) 


Using equations(5.51) and (5.59) through (5.62), the total 
lateness of sequence S is: 



(l+b)t22 + (l+b)t22 ♦ btj^2 ♦ 2 + btj^ ♦ 

niax(Z ^ ^^31* ^11 * ^21 * ^31 ^ ^41^ * ^^41 " 
4 

S d. (5.63) 

i«l ^ 
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he completion times of the jobs in sequence S' are: 


Cj(S') “ tj^j^ ♦ ^l2 " ^11 * ^^11 " 

€2(8') ■ max(Cj^, tj2 ♦ t2j} ♦ t22 

« max{ (1+b) tjj, tj ^2 ♦ t 22 ) ♦ bt 2 j^ 


(S.64) 


(5.6S) 


Cj(S') * max ( max [ C ^ ( ^11 ^ ^21^ ^^21* ^11 ^ ^21 ^ ^31^ 

♦ btji (S.67) 

C^(S') » max{max(roax(Cj , tj ^2 ♦ t 2 j) ♦ ^^21* ^11 ^21 

+ tjj] btjj, tjj + t2i ^ t^i * t^i> ♦ bt^j 

(S.68) 


sing equations (5.64) through (5.68) the lateness of the 


obs in sequence S' is: 


Li(S') 

L2(S') 

13 ( 8 ') 

L4(S') 


(5.69) 

(5.70) 

(5.71) 

(5.72) 


Jsing equations (5.49) (5.50) and (5.69) through (5.72), the 
total lateness of sequence S' is: 

4 

^ L^(S') - (l+b)tjj^ + X ♦ bt22 + Y ♦ btjj ♦ maxiy+btj^, 
^ll * ^21 ^31 ^41^ ^^41 ■ *^i (5.73) 


52 




Comparing equations (5*63) and (5*73) yields: 


Sl.(S) - • (l+b)t2i ♦ (l^b)t2i ♦ btji ♦ 2 ♦ 

i-1 ^ i-1 ^ i 


btjj ♦ max{2 ♦ ♦ t 2 i ♦ ♦ t^j^) ♦ bt 


- (l+b)tj^l - X - bt2j^ - Y - btjj - max{Y ♦ 

4 

^^ 31 * ^11 ■*■ ^21 ^ ^31 * ^ 41 ^ ■ **^41 


Reducing the terms brings to: 

4 4 

^L. (S) - (S') » (2+b)t,, ♦ Z + max{2 •*• bt,<i, t,i ♦ 


iTi iTi 1 


^21 ^31 ^ 41 ^ ■ ^11 ' ^ ^ 


^11 ^21 “■ ^31 ^41^ 


Using equation (5.54) there are two sets of terms to compare 


(2+b)t22 and tj^^ ♦ X 
Check the two cases of X 
1. X - (l-»-b)tj^ 

Compare the terms (2+b)t2j^ and tj^j ♦ (l+b)tjj 
It follows that (2+b)t-,i > t,, ♦ (l‘»-b)t,i as 



2. X - tji ♦ t2i 


Compare the terms (2-*^b)t22 tj^^ ^21 

It follows that (2+b)t22 > t^^ ♦ t^^ ♦ t2j 
The result is that (2+b)t22 > t^^ ♦ X 

II. Let's define and as follows: 

2j^ ■ 2 ♦ max(Z ♦ btjj, tjj^ ’ *21 ’ ‘n * *4i> 

Yj • Y ♦ max{Y ♦ btji, tjj . tjj ♦ r.jj ♦ t^j) 

I 

jusing equation (5.54) it is clear that: 


hth- 


It follows again that: 


4 4 

2l.(S) > £l.CS') 
i-1 ^ i-1 ^ 

It is shown that the SPT sequence brings to minimum the mean 
lateness . 


The Inductive Step 

After proving the theorem for the cases of two jobs, 
three jobs and four jobs, and by using the method of the 
mathematical induction: assume it is true for the case of 
m jobs and prove that it is true for the case of m<i’l jobs. 
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In general: the terms for the completion times of the jobs 

are as follows: 


C, - (l+b)t. 


“ max{Cj» ^11 ^ ^21^ ^^21 

Cj ■ max{Cj^, t^i ♦ t2j^ ♦ tj^} ♦ bt^ 


- ”«< CiK-l' * *’*' 


rhe total lateness is: 


••i «»• 

Sl. « 2c - 2d. 
i-1 ^ 1^1 ” iTl ' 


Given that < t 2 i < < t^.i < 

And proving for the cases i ■ 2, 3, 4 that 


2l^(SPT) < 2 L^(non-SPT) 


Assume that 


2‘sCSPT) < 2L.(non-SPT) 
i-1 ^ i-1 ^ 


It is necessary to check the case of ra+1 jobs. 
For the SPT sequence (tj^j < t 2 j^ < tjj^ <...< t^^j < 


f 


Vi' 


L* ♦ C* - d , 
m m+1 m+1 


for the non-SPT sequence (t2j^ < ^ ^31 ^ ^jn ^m+ 1 ^ 


Vi ■ 


^ ^Vl.l 


L + C ,-d , 

Bi+1 n in’*’l iM+1 


Observing the equations for and j 

in n I 

And adding them to 5 ^ L- (SPT) and $2 (non-SPT) respectively 

i-1 ^ i-1 ^ i 


to get the general result that: 


m+l m+1 

2 L. (SPT) < S L. (non-SPT) 
i-1 ^ i-1 ^ 


Summary 

This proof has shown that in a two machine LFP 
problem, the optimal solution is to sequence the jobs by thej 
shortest processing time on the first machine. The proof ofl 
optimality constructs the optimal sequence by mathematical 
induction. Figure S.S presents a block diagram that 
summarizes the steps of the proof. 


S6 


\ 

: ^ 


i 


A pairwise interchange of any pair of adjacent 
jobs i, and i>l where i^*! is following i and t.> 
ti^i* inproves total lateness of the sequence 6r 
iHves it the same 


The last pairwise interchange of a pair of adjacent 
jobs (going from a non*SPT to a SPT sequence) yields 
a strict improvement of the total lateness of the 
sequence 


Prove that it is true for the case of two jobs in 
the sequence 


Prove that it is true the case of three jobs in 
the sequence 


Prove that it is true for the case of four jobs in 
the sequence 


By using induct! on » assume it is true for the case 
of m jobs and then prove that it is true for the 
case of m-*^l jobs 


SPT sequence yields the minimum mean lateness in LFP 


Figure 5.5 


A Block Diagram that Summarizes 
the Steps of the Proof 
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5.2 Computational Results 


A computer program was written to simulate numerical 
examples that implement the theorem. The example that is 
presented contains four jobs to be performed on a two 
machine LFP. The information includes the amount of items 
in each job, i => 1, 2, 3, 4 and the machining factors, 
b., j « 1, 2. The specific figures are in Table S.l 

i 

Table S.l 

Data used in numerical example 


1 

2 

3 

4 

2 

3 

5 

7 

1 

2 

/ 

/ 

6 

9 

15 

21 


The input data contains 41 « 24 permutations or all 
possible sequences, and the output is the total lateness 
of each sequence. 

Analysis of the results show that the shortest 
processing time sequence, yields the minimum lateness, as 
proved by the theorem. 

A sensitivity test is conducted for each sequence, 
to find out how the total lateness of the sequence is 
changed, when a pairwise interchange of two adjacent jobs 




is done. The results of the sensitivity test show that 
each pairwise interchange of two adjacent jobs where a 
longer job preceds a shorter job, decreases the total 
lateness. Figure 5.6 describes a flow diagram that 
presents the algorithm that the 'sensitivity test follows. 
The sensitivity test is constructed as a network, where 
each block contains the sequence that is tested. The 
figures outside each block present the total lateness of 
the given sequence. There are twenty four networks. Each 
network shows how the lateness is decreasing after 
conducting a pairwise interchange of two adjacent jobs. 

The last block of each network presents the SPT sequence 
which is shown as was proven, to yield the minimum mean 
lateness. Each network is accompanied with a graph that 
shows the relations between the lateness and the number of 
pairwise interchanges of pairs of adjacent jobs. Each path 
in the network is presented as a different line in the 
graph. The curves show decreasing values of the total 
lateness for each interchange. I 









1 


2 


Number of pairwise 
interchanges 


Figure 5. 


(b) 


7 Results of Sensitivity. Test for 
Sequence 1 

(a) Network of sequence 1 

(b) Total lateness curve for 
sequence 1 




Total Lateness 


(2. 5,3. 7 M 2,3,5^_7j 
32 25 

(a) 



1 Number o£ pairwise 

Interchanges 


(b) 


Figure 5.8 Results of Sensitivity Test for 
Sequence 2 

(a) Network of sequence 2 

(b) Total lateness curve for 
sequence 2 
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Total Lateness 


i 2.3.5.7~1 SPT Sequence 
2S 


Figure 5.9 The Optimal Sequence 


I 2, 3,7,5 I -I 2.3.5.7 I 
29 25 


(a) 



1 Number of pairwise 

interchanges 


(b) 


Figure 5.10 Results of Sensitivity Test for 
Sequence 4 

(a) Network of sequence 4 

(b) Total lateness curve for 
sequence 4 
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Number of pairwise 
interchanges 
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Figure S.ll Results of Sensitivity Test for 
Sequence 5 

(a) Network of sequence 5 

(b) Total lateness curve for 
sequence 5 
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2 


3 


Number of pairwise 
interchanges 


(b) 


Figure 5«12 Results of Sensitivity Test for 
Sequence 6 

(a) Network of sequence 6 

(b) Total lateness curve for 
sequence 6 




Total Lateness 



Number of pairwise 
interchanges 


(b) 


Figure 5.13 Results of Sensitivity Test for 
Sequence 7 

(a) Network of sequence 7 

(b) Total lateness curve for 
sequence 7 







Total Lateness 
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31 25 

(a) 



1 Number of pairwise 

interchanges 
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Figure 5.15 Results of Sensitivity Test for 
Sequence 9 

(a) Network of sequence 9 

(b) Total lateness curve for 
sequence 9 
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Total Lateness 



Figure 5.16 Results of Sensitivity Test for 
Sequence 10 

(a) Network of sequence 10 

(b) Total lateness curve for 
sequence 10 
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Figure 5.17 Results of Sensitivity Test for 
Sequence 11 

(a) Network of sequence 11 

(b) Total lateness curve for 
sequence 11 
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4 Number of pairwise 
interchanges 




Figure 5.18 Results of Sensitivity Test for 
Sequence 12 

(a) Network of sequence 12 

(b) Total lateness curve for 
sequence 12 
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1 2 Number of pairwise 
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Figure S.19 Results of Sensitivity Test for 
Sequence 13 

(a) Network of sequence 13 

(b) Total lateness curve for 
sequence 13 
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Total Lateness 
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Figure 5.20 Results of Sensitivity Test for 
Sequence 14 

(a) Network of sequence 14 

(b) Total lateness curve for 
sequence 14 









Figure 5.22 Results of Sensitivity Test for 
Sequence 16 

(a) Network of sequence 16 

(b) Total lateness curve for 
sequence 16 


75 












Figure 5.23 Results of Sensitivity Test for 
Sequence 17 

(a) Network of sequence 17 

(b) Total lateness curve for 
sequence 17 
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Figure 5.25 Results of Sensitivity Tost for 
Sequence 19 

(a) Network of sequence 19 

(b) Total lateness curve for 
sequence 19 


82 




Figure S.26 Results of Sensitivity Test for 
Sequence 20 

(a) Network of sequence 20 

(b) Total lateness curve for 
sequence 20 


83 
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wise interchanges 
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Figure 5.27 Results of Sensitivity Test for 
Sequence 21 

(a) Network of sequence 21 

(b) Total lateness curve for 
sequence 21 











Figure S.28 Results of Sensitivity Test for 
Sequence 22 

(a) Network of sequence 22 

(b) Total lateness curve for 
sequence 22 


89 










Figure 5.29 Results of Sensitivity Test for 
Sequence 23 

(a) Network of sequence 23 

(b) Total lateness curve for 
sequence 23 


92 











Figure 5.30 Results of Sensitivity Test for 
Sequence 24 

(a) Network of sequence 24 

(b) Total lateness curve for 
sequence 24 











Total Lateness 





5.3 Flow Ti— In « Linear Flow Shop 

In general, flow time is defined as the amount of time 
that a job spends in the shop. 

In the case of a two machine LFP system shotm ip 
Figure S.31, the flow times Fj^, k e I are as follows: 

h • hi * h2 

Fj - maxCCj. t^^ ♦ ♦ tj2 

F3 - max{C2. ♦ tj2 


(5.74) 

(5.75) 

(5.76) 


Using equations (5.74) through (5.76), the flow time of job 
i is: 


•■i 


max{C 


i-1* 




} ♦ t 


i2 


(5.77) 


where 

t^j is the time to perform job i on machine j 
is the completion time of job i 
is the flow time of job i. 

Theorem 2 

Vhe mean flow time F in a two machine LFP where the 
job times on the richines are a linear function of the 
batch sizes (t^^^ ■ ^j°i* ^j 1 minimized by shortest 

processing time (SPT) sequencing. 
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Proof 

The proof is a mirror image of the proof of Theorem 1, 
and utilizes an adjacent pairwise interchange argument. 

Consider a sequence S that is not a SPT sequence. That 
is, somewhere in S there must exist at least one pair of 
adjacent jobs, i and i+1 with i**-! following i, such that 
Now construct a new sequence S', in which jobs 
i and i***! are interchanged in sequence and all other jobs 
are not changed. Ths situation is depicted in Figure S.l 
where Ty^ denotes the point in time at which job i begins in 
sequence S and at which job begins in sequence S'. 

"A" denotes the set of jobs that precede job i and i+1 in 
both schedules and "B" denotes the set of jobs that follow 
i and i+1 in both schedules. „ 

The proof is shown for ^ F. (sum of the flow time over 

i-1 ^ 1 ® 

all the jobs), as the mean flow time is: 7 * - z- 


I'l’k - ♦ ''i * * ''b 


( 5 . 78 ) 


ifhere 


Fy^ is the sum of the flow times of all the jobs in 


set "A" 


Fj is the flow time of job i 
Fj is the flow time of job j 

F^ is the sum of the flow times of all the jobs in 
set "B". 
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The flow tiroes of the jobs in sequence S are: 


FjCS) - maxiC^, j) ♦ j 


(S.79) 


Fi^.l(S) • roax{C^(S), ♦ t^^j ♦ ^i-t-1,1^ ”*■ ^i+1,2 


(5.80) 


sing equation (5.78), the flow tiroe of sequence S is: 


Fjj(S) “ F^ ♦ roax{C^, ♦ ^i,l^ * ^i,2 ®®3c{Cj^(S) 


* ^i,l * ^i+1,1^ ^ ^i+1,2 * 


(5.81) 


The flow tiroes of jobs i and i***! in sequence S' are: 

- lllax{CJ^, ♦ Vi,2 (5.82) 

Fj(S')- max{C.,j(S'), T/, * 1 * t. j) ♦ t. 2 (5.83) 

! 

Rising equation (5.78), the flow time of sequence S' is: 


Fj^(S') . F^ . max(C^, ♦ t.^j j) . tj^^ ^ * 
max{C.^l(S'), Ta ♦ tj^jj . t. j) 


^ *i,2 * 


(5.84) 


101 




When observing equations (5.79) through (5.84), one can 
notice that the terms are actually equal to terms of 
completion times in equations (5.7) through (5.12) that are 
developed in Theorem 1. 

The proof of the sequence that generates the minimum 
mean lateness (Theorem 1) is based on the completion times 
of the jobs, and the sum of the completion times of all the 
jobs in the sequence. Using the same steps of the proof 
can show that the mean flow time is also minimized by the 
shortest processing time sequencing. 

It is proven that the mean flow time is minimized by 
the SPT sequence. 

5.4 Waiting Time in a Linear Flow Shop 

In a flow shop scheduling problem, a waiting time can 

be defined for each jeb in the system. It is the time that 

every job must wait until it can be processed. If the jobs 

are processed in a numerical order, the waiting time for 

the i^ object is 5^ tv i • The idea is similar to the 

k-1 

problem of scheduling for single stage production, since in 
flow shop problems there is no idle time on the first 
machine. 

A performance measurement can be the sum of all 
m i-l 

waiting times 2 2 t^ . . The minimum of this expression 

i-l k»l 

occurs at the same time as the minimum of the expression 
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m i 


Y, Z) tv 1 occurs, since the two expressions differ by the 
1=1 k-1 

n 

constant Y t . . 

ft ^ 


theorem 3 


The sum of waiting times of all the jobs in a flow 
shop problem is minimized by shortest processing time (SPT) 
sequencing. 


Proof 


m 1 


In the proof , the term Y S t. , is used because of 

ft k-1 

:he reason mentioned above, as the term that represents the 


sum of waiting times. 


m 1 


t. , • 2(n+l-k)t. , 
i-1 k-1 k-1 


ccording to Hardy, Littlewood and Polya (1952), this sum is| 
least when the series are monotonic in opposite senses. ' 
herefore the sum of all waiting times is minimized, if 
the jobs are sequenced according to SPT (tj^ 1 ^2 1 

..m,j 
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CHAPTER 6 


SUMMARY AND CONCLUSIONS 


This study investigates a flow shop scheduling 
problem which is defined as the linear flow shop problem 
(LFP) . The purpose of the research was to find an optimal 
solution for real life cases that can be defined as linear 
flow shop problems. Three performance objectives were 
defined; mean lateness, mean flow time and waiting time. 

The results show that the shortest processing time sequen- 
cing minimizes the objectives that were defined. Imple- 
menting optimal solution for a practical problems is very 
easy. The jobs that have to be processed are ordered in the 

i 

system by monotonic increasing performance time. 

Chapter 1 of this dissertation introduces to the reader 
the scheduling and sequencing area. The basic concepts and j 
basic models of scheduling are defined and the main 
performance measurements of a scheduling problem are 
formulated. 

Chapter 2 defines the basic concepts of flow shop 
scheduling. Two variations of flovr shop scheduling are 
discussed, the pure flow shop and the general flow shop. 

The only requirement on the schedule is that all movements 


104 





of jobs bdtween machines within the shop be in uniform 
direction from machine j to machine j***! etc. 

Chapter 3 reviews past research in the flow shop area. 
Many studies were conducted, and many articles were 
written on flow shop problems. The studies that are 
reviewed introduce the basic ideas that were developed in 
the machine scheduling area and the special flow shop cases 
A discussion of the differences between optimal and heuris* 
tic solutions to the problem is included. Optimal solution 
are efficient in theory but difficult to implement in 
practical situations. 

The linear flow shop problem is formulated in Chapter 
4. It is a model of a practical production facility where 
the time to complete a batch on a specific machine and the 
times to complete the same batch on all other machines is a 
linear function of the batch size. The time to perform all 
jobs through all machines and the machining time required 
to perform the jobs is developed. The idle time involved 
in a linear flow shop problem is discussed in three steps. 
First the case of a two machine LFP, second a three machine 
LFP and third the consideration of the general m jobs n 
machines case. Upper and lower. bounds for performance time 
are developed where the optimal solution falls in that 
range. 

Chapter 5 describes three performance measurements of 
the linear flow shop scheduling, mean lateness, mean flow 
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time and waiting tine. Three theorems regarding the definec 
performance measurements are stated. The optimal solution 
for a two machine case is proven to be a schedule that 
follows the shortest processing time ordering. The proofs 
for mean lateness and mean flow time use the method of 
pairwise interchange of pairs or adjacent jobs in the 
sequence. The method of mathematical induction is used to 
prove that the theorems hold for the case of m jobs in the 
sequence. A nximerical example follows the proofs. A 
sensitivity test is conducted to show how the shortest 
processing time sequencing yields the optimal solutions. 

This study presents optimal solutions for the linear 
flow shop problem. It was found in many studies that 
optimal solutions are difficult to implement in practical 
cases because of the complexity of the formulation and 
computational problems. The advantage of the results found 
in this dissertation is that it is very easy to construct 
the optimal sequence according to the shortest processing 
time ordering. The results can be implemented in practical 
production situations that fit the linear flow shop 
definition. 

This research investigated the linear flow shop where 
all variables are deterministic and can be established 
according to data available from past years or experience 
of production management people. An open subject to be 
researched in the future is the case where the machines' 
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coefficients are random variables and the times to perform 


the jobs are the expected performance time. Further 
research should be conducted into production models other 
than the LFP. For example, a fixed setup time per machine 
is often necessary in parts production. Non linear 
production times often occur in chemical processes. 

Another issue in scheduling is the idle time involved in th< 
process. One can take the advantage of idle time of 
machines to conduct maintenance needed in the shop. In all 
kinds of problems created in the scheduling area one should 
keep in mind the way of implementing the results of the 
research. People in industry want easy and understandable 
procedures to make their production lines run smoothly and 
i efficiently. 
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ROBECON 

A GENERALIZED NETH00OL06Y 
for 


ASSESSING THE ECONOMIC CONSEQUENCES 
of 

ACQUIRING ROBOTS FOR REPETITIVE OPERATIONS 


ABSTRACT 


AUhough problems reldting to the engineering design of robots are 
awesome, it is the economic aspect which is fundamental to the user's 
decision to acquire robots for •'epetitive operations. Nevertheless, a 
review of the relevant literature suggests that very little exists In 
the way of providing guidance to prospective purchasers of robots as to 
the economic consequences of prospective Ov^quisltions. This conclusion 
stems from matching critiques of more than 25 published references against 
a set of explicit criteria for a generalized methodology. Recognizing 
the need for an appropriate methodology, an exhaustive set of cost elmnents 
are identified which are to be included in a comprehensive analysis. 

(This paper is intended to serve as the first in a series leading to 
a fully-develot.e<j model, ROBECON, which may be used for specifying the 
econcmilc consequences of robot systems acquisitions. The model will be 
computer based and user interactive.) 
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R06EC0N 

A GEHERALIZED METHODOLOGY FOR ASSESSING THE ECONCniC CONSEQUENCES 
OF ACQUIRING ROBOTS FOR REPETITIVE OPERATIONS 

I. INTRODUCTION 

Robots Defined 

The Robot Institute of America defines the robot as "a programnable, 
multi 'function manipulator designed to move material, parts, tools or 
specialized devices through variable programmed motions for the performance 
cf a variety of tasks.” A less precise definition has been adopted by a 
manufacturer- users group. Computer Aided Manufacturing International (CAM>I): 
"a device that performs functions ordinarily ascribed to human beings, or 
operates with what appears to be almost human Intelligence.” With either 
definition. It Is commonly understood that moder'n robots are proqramnable 
manipulators that can perform useful work automatically without hirnian 
assistance. (The term robot comes from a Czech word for forced labor; It 
was Invented for Karel Capek's 1921 melodrama, R.U.R. ). 
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Robot Installations: Substantial and Growing 

Beginning with the development In the fflld-60's of the microprocessor, 
which permitted robots to be made smaller and cheaper, and spurred by 
endemic wage Inflation, robots have been used with Increasing frequency 
In the Industrialized nations. There are differences of opinion as to the 
nundjer of robots currently in place around the world. One estimate Is 
that there are “about 7,000 working Industrial robots world-wide.” 
(Ferguson, October 12, 1980). AnoUier estimates about 15,000 robots 
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in the Western industrialized nations, with 10,000 in Japan and 3,000 
in the U.b. at the end of 1979. ( TIME , Dec. 8, 1980). Still another 
source estimates 40,000-50,000 in worldwide use, with 30,000 of these 
installations in Japan (Allan, 1979). 

Expert opinion appears unanimous that the forecasted growth of robot 
installations will be spectacular into the foreseeable future. It is esti- 
mated that installations will increase at the rate of 30%-40% over the next 
decade. (The first and largest of the robot manufactui'ers in the U.S., 
Unimation Inc. of Danbury, Connecticut, experienced a 30% per annum growth 
rate over the past seven years.) Estimates of industry salos potential range 
from $2 billion to $4 billion by 1990. (Currently, industry sales in the 
U.S. are about $90 million.) A recent forecast by the Society of Manufactur- 
ing Engineers and the University of Michigan estimates that by 1987, 15% of 
all assembly systems in the U.S. will use robot technology. 

There are several significant reasons underlying expectations for sub- 
stantial growth of robot installations in the foreseeable future. First, 
the conditions which led users to adopt robots over the past decade will 
persist, principally with respect to higher wage rates. Second, unit costs 
can be expected to decrease because robots are becoming smaller and more 
flexible and new manufacturers are entering the industry. Third, applica- 
tions will increase as the functional capabilities are expanded, especially 
with respect to the ability of robots to see properly the articles which 
they are manipulating. 

The Problem 

Th 2 engineering design aspect of robots is awesome, yet it is the economi c 
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aspect which Is fundamental to the user's decision to acquire this equipment. 
After all, robots generally perform no functions which cannot otherwise be 
performed by combinations of human workers, machines and devices. The decision 
to acquire robots Is Influenced, wholly or In part, by the economic consequences 
to be expected from that decision. A preliminary review of the literature 
suggests that this Issue has received little attention relative to the design 
and operational characteristics of robots. Certain cost estimates are widely 
quoted In the literature*, but these are generally Inadequate as a guide to 
prospective users who may be contemplating capital Investments of $5,000 to 
$150,000 per Installation. (Multiple Installations, I.e., Implementation of 
systems using two or more robots, are not uncommon. Capital Investments In 
the millions of dollars may be required In these Instances, of course.) 

Large, relatively sophisticated firms will probably have the expertise 
"In-house" to conduct appropriate economic analyses. However, as robot 
Installations become more extensive, it Is likely that smaller, less sophis- 
ticated firms will be considering the acquisitions cf robots, and they will 
need competent guidance as to the economic justification for these decisions. 

It is this issue which provides the justification for the research described 
in the following sections. 

Objective 

There are a variety of ways of describing the process by which prospective 
users arrive at the decision to acquire a specific robot or robotic system. 

*For example, Unimation Inc. reports that a robot's cost Is $4.60 hourly, 
and this has remained relatively constant since 1961. This is a rough 
estimate, however; it is based on straight line depreciation rather than 
cost of capital recovery, and taxes are ignored. 
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For our purposes, here, we may focus on three principal stages. First , the 
appropriate decision maker(s) within the firm must focus upon a limited set 
of candidates from among the much larger population of robots (and related 
auxiliary equipment and software) currently available In the marketplace. 

(It Is assumed, at this point, that the decision maker has already completed 
an analysis of the task(s) and operating enviroranent and Is reasonably 

I 

convinced that a robot system may represent an optimal solution to the roanu- i 

factoring* problem). At this stage It will be necessary to describe Important 
technical requirements for the robot(s). Including: capacity, drives and 
controls, memory, and other features such as tactile, feedback and visual 
sensors. These technical requiranents must then be matched against availability. 

Seo Exhibit 1, for example. The central feature of this first stage Is the 
Identification of a set of candidate systems with technical characteristics 
suitable to the firm's operational requirements. This Includes, In addition 
to the robots themselves, associated requlronents such as changes necessary 
to other equipment, tooling, spare parts and test equipment for maintenance, 
utilities, back-up equipment to be used If and when the robot Is down, safety 
equipment, and the like. 

The second stage is an economic analysis of the consequences, or Impacts, 
of the candidate robot systems as Identified In the first stage. This Is the 
focus of the research described here. The objective is the development of an 
evaluation methodology which will permit users to forecast, or assess, the | 

! .4 

economic consequences of acquiring one or more robots for repetitive operations . 


*Here, "manufacturing" includes fabrication, assembly, inspection, material 
handling and other tasKs associated with the production of manufactured 
goods. 
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Exhibit 1. Characteristics of Selected Robots 

^ds in Robots", article appearing in Industrial Robots. Vol. 1. William T. Tanner, Fditnr 
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Users are assumed to be any business flms (manufacturers, fabricators, 
processors, etc.) or governmental agencies who may be considering the purchase 
of robots as operational alternatives and for whom the economic consequences 
are relevant to the acquisition decision. 

The third phase In this process, as Illustrated In Exhibit 2, Is one In 
which economic consequences are considered jointly with other (non-ecoiiomic) 
consequences so as to arrive at a choice from among alternative systems. 

There are a variety of approaches to this "multiple criteria" probTen, some 
of which are relatively complex. In any event, this Is not an Issue which 
we will address further at this time. The research described In this docu- 
ment focuses only on Stage 2, the economic analysis. 

Tasks 

The tasks necessary to meet the objective described above are as follows: 

(1) Specify the criteri a which must, or should, be met by the economic 
analysis methodology. Including mathematical model(s) and associated 
procedures . 

(2) Review the literature to determine the extent to which economic 
data, models, and analytical procedures are currently available 
to prospective users. Relevant references will be critiqued In 
view of the criteria Identified In (1). 

(3) Identify the elements of total system costs , that Is, the economic 
Impacts which. In general, may result over the lifetime of the 
robotic system, from Installation to ultimate disposal. In certain 
gove«*nment applications, and In the U.S. Department of Defen;^e In 
particular, these are known as life cycle costs (LCC). (See, for 
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Exhibit 2 . A Simplified Schematic Representation 
of the Decision Process for 
Acquiring a Robotic System 



U.S. Department of Defense, 1978; U.S. Department of the 
Air Force, 1978; Graver and Jenkins-Stark, 1976; and Kolarik, 
1980 .) Total system costs Include subsequent costs as well as 
the Initial Investment. 

(4) Develop engineering cost estimates (ECEs) and/or cost estimating 
relationships (CERs) , as appropriate. These are two principal 
approaches to the estimation of future economic consequences 
which are likely to result from a current Investment decision. 
Generally, the CER approach relates system costs to a combination 

of measures of the systems (dimensions, performance characteristics, 
etc.)*. The cost estimating relationships are obtained through 
curve fitting techniques. In the ECE method, total system costs 
are broken down Into rleatively small components, or elements; 
the elements are related by ECEs which reflect the ways In which 
the system Is developed, operated and maintained. It Is expected 
that these relationships will be of sufficient generality to permit 
their use in a wide variety of analyses. 

(5) Develop computer-based models for generating Total System Costs . 
Three separate models will be developed, with Increasing complexity, 
utilizing the ECEs and/or CERs as determined In the preceding task. 
The mdoels are characterized as follows: 

(a) In the first model we assume that all economic consequences 
(the amount and timing of cash flows); as well as the planning 


*Th1s method is sometimes called parametric costing . 
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horizon and discount rate, are deterministic. It Is also 


assumed that the characteristics of the robotic system 
currently under consideration for Implementation (the 
“challenger'') are Identical to those of all future challengers. 
That Is, all future challengers are Identical to the current 
challenger . 

(b) In the second model we relax the assumption concerning the 
verisimilitude of current and future challengers. We now 
assume that future challengers are not necessarily Identical 
to the current challenger . Indeed, It Is likely that certain 
costs will decrease due to economies of scale and the ability 
of robot manufacturers to move out on the learning curve, for 
example; other cost elements, such as energy and labor, might 
be expected to Increase over time. 

(c) The third model differs from the second In that all economic 
consequences, the planning horizon and the discount rate are 
assumed to be stochastic. Unlike the prior models wherein 
all Inputs were assumed to be known with certainty, we now 
treat these Impacts as random variables. 

(6) Test the ImplementabiUty of each of the models through a series 
"of controlled" experiments. These will consist of a set of sample 
problems which will be presented to real-world decision makers for 
solution using the models. The experiences of these decision 
makers will be monitored and evaluated to determine the extent to 
which the models are useful in the capital allocation decision. 
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II. CRITERIA FOR A GENERALIZED METHODOLOGY 

Prior to developing a generalized methodology for assessing the economic 
consequences of acquiring robots for repetitive operations, it is necessary 
to make explicit the criteria by which the efficacy of the methodology will 
be measured. These same criteria can also be used systematically to critique 
the existing relevant literature. For our purposes, then, the following 
criteria will be established: 

(1) Theoretically sound — We are not interested solely in obtaining 
a solution. The solution must be internally consistent with the 
decision maker's (user's) objectives as well as the assumptions 
underlying the model. 

(2) Credible — The users must have a feeling of confidence that the 
methodology will in fact provide solutions that are useful in the 
decision making process. The methodology must be believable. 

(3) Verifiable — The user should be able to replicate, or verify, 
the results by tracing the chain of events from data input to 
ultimate solution. Verifiability is a preconditicn to credibility. 

(4) Comprehensive — The economic model (s) imbedded in the methodology 
should include all the economic impacts which can reasonably be 
expected to occur as the result of the decision. (The time interval 
over which these impacts will occur is the planning horizon .) Thus 
the methodology should include the economic consequences of the 
total system — equipment acquisition, operation and maintenance, 
taxes, and the like — throughout the entire planning horizon. 

This is the Total System Costs concept. 
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(5) Reasonable data requlreiwents — Although comprehenilblllty Is a 
desirable. If not essential, element of the assess*nent method- 
ology, It Is unrealistic to expect that the analyst will be able 
to deal exhaustively with absolutely all economic Impacts. To 
do so Is neither possible nor desirable. The data requirements 
for the economic models should be limited to only those which 
are likely to have a significant affect on the user's capital 
allocation decision. The cost of gathering Impact drta and 
exercising the models should In no case exceed the economic ad- 
vantage to be gained from the analysis. 

(6) Accuracy — The level of accuracy should not exceed that which Is 
necessary to Identify significant differences among alternatives. 

(7) Assumptions made explicit — The assumptions underlying the nethodology 
and Imbedded in the analytical models should be stated clearly. 

(8) Important factors stressed — Not all elements of the analysis 
are of equal Importance. Those which have greatest significance 
should be highlighted. 

(9) Uncertainly treated explicitly — Equipment acquisition decisions 
are properly based upon anticipated consequences expected to result 
from the various alternative courses of action. These consequences 
lie In the future, and hence are uncertain. (Some would argue that 
the more distant the event, the greater Is the uncertainty, but 
this is not necessarily so.) The extent to which this uncertainty 
affects the decision should be made explicit so that It may be 
treated by the decision maker as a separable Issue. 
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(10) Incorporates efficiencies over time — The learning curve 
(Improvement curve, progress curve, etc.) has been used for more 
than forty years to describe the relationship between productivity 
(cost/quantity) and time. During the Initial stages of production. 
In particular, productivity Is Improving as the people and machines 
In the process "learn" to operate more effectively. Economic models 
should Incorporate this effect. 

(11) Reflects real and relative price changes — Economic Impacts 
should not be expected to remain constant over time, particularly 
over a long planning horizon. In part these differences result 
from changes In the relative prices of specific goods and services, 
popularly known as Inflation . Inasmuch as relative price changes 
may be of significance to the capital allocation decision, they 
should be Incorporated Into the analysis. This Is eipeclally 
Important for those goods and services for which prices change at 
substantially different rates. 

III. LITERATURE REVIEW 

During the summer of 1981 an Intensive review of the literature was 
conducted to Identify the extent to which published material describing the 
economics of robotics Is available to prospective users. Sources for review 
Included newspapers and popular magazine articles, anthologies (especially 
W.R. Tanner's Industrial Robots) , professional conference proceedings, 
goverrnnent reports, and technical p3oers of professional societies (especially 
the Society of Manufacturing Engineers). Consultants working In this field 
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Mere also contacted for ^cads. More than 200 Individual Items were 

c 

reviewed; the references appearing In the Bibliography are representative. 

Of these, only the dozen listed In Exhibit 3 are directly related to economic 
analyses of robot Installations. 

The Accounting Method 

^ As Indicated In Exhibit 3, these references m^y be characterized by 

one or more of several analytical procedures. The accounting method describes 
economic consequences (costs and benefits) In accounting terms, that Is, the 
effect of the Installation on the firm's Income and expense accounts. Thus 
the cost of capital recovery Is defined by annual depreciation expense.* 

The principal objection to the accounting method Is that the opportunity 
^ cost Is Ignored. The opportunity cost, sometimes described as the minimum 

attractive rate of return. Is the return which would be expected from alter- 
native Investment opportunities should the specific project proposal not be 
i ( funded. As described In the literature of engineering econooiy, the concept 

j of capital recovery (CR) Incorporates the opportunity cost as follows^ 



: c 



Discounted 

Accounting Payback Cash 

Method Method Flow 



Exhibit 3. Publications Related to Economic Analysis 
of Robot Installations 




CR • (C-L)(A/P, <. N) ^ Li 
C ■ initial cost 

L ■ net ia1vag«i (residual) value at the end of N periods 
i ■ opportunity cost (discount rate) 

N ■ service life of the investment 
and (A/P, 1* N) » functional form of the algebraic expression 

liJiiil! — 

(HI)” - 1 

It may be shown* in general* that capital recovery does not yield the 
same results as those derived from the popular depreciation methods. To 
illustrate, consider straight line depreciation . The annual depreciation 
expense (0) is given by: 

0 ■ ('d - 

where « cost basis 

L^ * expected salvage value for depreciation purposes 
Nj » depreciable life 

To simplify our exmnple, let us suppose that C « L ■ and N > N^. 

It may be shown that the percent error (a) is given by 

^ ‘ ^ ‘ Cf-pi(A/Al^. U] i pi 

where p » L/C 

The percent error (a) is shown graphically in Exhibit 4 for N » 5. The 

error Increases with the discount rate and the ratio of salvage value to 

Initial cost. When 1 > 20% and p « 0 (no salvage value), for exan^le, 
the error Is approximately 40%. When 1 ■ 20%, and p « 0.50, the error Is 
about 62%. (Note that 0 ■ 0 and CR » 1 for all values of N when p » 1.00. 
Thus, In this special case, a ■ 100% for all values of N.) 



Percent Error (a) 




The Payback Method 

As illustrated in Exhibit 5, payback (or payout ) is the number of 
periods required for cumulative benefit; to exactly equal cumulative costs. 
Costs and benefits are usually expressed as cash flows, although discounted 
present values of cash flows may also be used. In either case, the payback 
method is based on the assumption that the relative merit of a proposed 
investment is measured by this statistic. The smaller the payback (period), 
the better the proposal. 

Despite the apparent fact that the payback method is widely used in 
industry, it suffers from serious theoretical deficiencies. The most im- 
portant of these is that the payback method ignores the consequences of 
the proposed investment after the period in which payback is completed. 

This may be shown with reference to Exhibit 6. "ere we have two competing 
projects. Alternatives A and B, with payback for A less than that of B. 

But it is unlikely that A would be preferred to B since the latter generates 
far greater net cash flows over the remaining periods in the planning horizon. 
With verv rare exception, payback should not be used as the sole criterion 
to measure economic efficiency. 

Discounted Cash Flow (DCF) 

There are a number of variations to the discounted cash flow method: 
present worth or present value, equivalent uniform annual cost, rate of 
return (or return on investment), benefit-cost ratio, and the like. They 
have in common recognition of the timing as well as the amounts of cash 
flows; money has value over time because of the existence of alternative 
investment opportunities. When used properly , the DCF variants lead to 
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Exhibit 6. Payback for Two Competing Investment 
Alternatives 


anmng Horizon 


consistent solutions. 


Our literature search revealed seven publications In which the DCF method 
Is applied to robotics Investments. These are: 

Abraham and Beres, 1979 . Returns on Investment (ROIs) are summarized 
for ten candidate assembly equipments and tooling as well as fifteen 
separate combinations. Uncertainties are Ignored. There Is no description 
of the process by which the ROIs are computed. 

Behunlak, 1979 . This very brief paper describes economic analyses for 
three robotic applications: swaging, die casting and painting. Cash flow 
tables are shown for 5-year planning horizons. Three after-tax figures of 
merit are shown for each application: (1) payback, (2) return on Investment, 

and (3) discounted rate of return. The formulae, or procedures, for computing 
ROI and DCRR are not Included, and the author's results cannot be verified. 
Uncertainties are Ignored. 

Behuniak, 1980 . This Is similar to the author's earlier paper, except that 
no cash flows table Is Included. Results are reported for a die casting 
robot with Initial cost of $68,000. Payback =3.4 years; ROI = 24%. 

(For the die casting robot In the author's 1979 paper, payback * 3.6 years 
and ROI = 29.2%.) 

Bublick, 1979 . This paper described an economic analysis for robots used 
In sp»*ay coating and finishing. Cost details are provided for both the 
manual and robot alternatives. The procedure for determining payback 
(1.17 years In this case) Is detailed. A cash flow table is provided for 
a 7-year planning horizon. Including an adjustment for Inflation of 6% per 
annum. The author claims that ROI = 121% for this application, although 
the formula for determining ROI Is not given. Uncertainties are Ignored. 
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Hanity and Belcher . This Is perhaps the most interesting of the papers 
reviewed. Dennis Hanify is affiliated with the ITT Research Institute; 

J.V. Belcher is with the Universal Oil Products Company. The paper 
describes the Industrial Robot Analysis (IRA) program initiated at the 
ITTRI Robot Technology Center and currently available to clients interested 
in the technical and economic effects of robot systems. An example problem 
is included — a "real world" example — but cost figures are fictitious 
since the original data are proprietary. Major costs, their timing and 
"probable variations" are given for this example problem. Calculations 
are performed using a general purpose computer program, Economic Systems 
Analysis (ESA), to compute the prospective rate of return, before or after 
income taxes. Uncertainty is addressed by evaluating the effects of using 
the optimistic or pessimistic cost estimates. The procedure used to determine 
rate of return is not detailed in the paper: no equations, no flow diagram, 
no computer program. It is not possible, therefore, to critique the method- 
ology further. 

Stout, 1973 . Assumed costs and benefits for a "typical" project are plotted 
as a function of time and the ROI calculated. ROI is defined (properly). 
Tanner. 1978. The author asserts that "simple 'rules of thumb' and a stream- 
lined cost analysis method can be applied to determine the potential economic 
return of a contemplated robot installation." A cost analysis form is included 
with numerical examples illustrating the calculation of: 


(1) 

and (2) 


Pawhart - Total Expenditures 

" Total Annual Net Future Savings 

ROI is that value of r such that: 


Current Value 
of Savings 




Total Annual Net Future S>v i ngs 
(1 + r)" 
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The calculation In (2) Is described as "ROI for the Depreciated (sic) Cash 
Flow Method". An alternative calculation for ROI Is given by 


(3) 


ROI = 


Profit After Taxes 
Investment Base 


The author's procedure assuines that cash flows remain constant over the 
planning horizon. Thus there Is no cash flow table for other than the 
“typical" year. Uncertainties are Ignored. 


IV. LIFE CYCLE (TOTAL SYSTEM) COSTS 

As a general principle, the economic consequences of proposed Investment 
In a robotics system should Include all significant costs that are likely to 
result from the Investment. ("Benefits" are reductions In costs that may be 
obtained when comparing any pair of alternatives, and hence are Implicit In 
this principle.) Economic consequences, or Impacts, should be measured over 
the total life cycle of the proposed syst«n.* Impacts should be estimated 
for each of the candidate systems as well as the existing process. 

For our purposes, economic consequences may be grouped Into three 
broad categories, as follows: 

!• Plant and Equipment 

1.1 The robot(s). Including sensors and Interlocks 

1.1.1 Initial cost 

1.1.2 Service life (not a cost) 

1.1.3 Residual value (net salvage value) at the end of the 
service life. 


♦In the literature of economic analysis, life cycle Is frequently referred 
to as the planning horizon. Strictly speaking, however, the planning 
horizon may be longer than the system life cycle, especially If signifi- 
cant economic consequences persist beyond the end of the life of the 
system 
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1.2 Associated tooling 

1.3 Spare pa**ts 

1.4 Property taxes 

1 .5 Insurance (property only) 

1.6 Energy reoul rements 

1.7 Tax consequences 

1.7.1 Investment credit 

1.7.2 Tax savings due to depreciation 

1.7.3 Gain (loss) on disposal 

1.8 Space requirements 

1.9 Installation (including rearrangement of existing facilities) 

1.10 Safety equipment (protective clothing, etc.) 

1.11 Programming 

1.12 Modification of existing equipment to ensure compatability with 
robot(s). 


2. Operation and Maintenance 

2.1 Operating labor 

2.1.1 Salaries/wages 

2.1.2 Fringe benefits (costs to employer) 

2.2 Maintenance labor (for periodic maintendnce) 

2.2.1 Salaries/wages 

2.2.2 Fringe benefits (costs to employer) 

2.3 Direct cost of injuries and illness (hospitalization, medical 
care, etc.) 

2.4 Absenteeism (cost of lost productivity) 
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2.4.1 Illness 

2.4.2 feigned Illness 

2.4.3 injury 

2.5 Training 

2.6 Supervision 

2.7 Insurance (personnel only) 

2.8 Overtime (not included in 2.1 and 2.2 above) 

2.8.1 Operating labor 

2.8.2 Maintenance labor 

2.9 Labor turnover 

2.9.1 Termination 

2.9.2 Recruitment 

2.9.3 Training 

2.10 Retooling and set-up costs for batch processing 

2.11 Maintenance tools and supplies 

2.13 Documentation (operation and maintenance) 

2.13 Costs of interrupted production not included in 2.10, especially 
down time. 

3. Product 

3.1 Required changes in product design 

3.2 Raw material requirements 

3.3 In-process inventory 

3.4 Effects of production rate on: 

3.4.1 other plant activities 

3.4.2 shipping schedules 


- 24 - 


3.5 Defective (sub- standard) product 

3.5.1 Scrap rate (not a cost) 

3.5.2 Net cost of handling and reworking defective product 

3.5.3 Costs due to undetected defective product released to 
customer (e.g., loss of good will, responding to customers' 
complaints, replacing returned products) 

In addition to the above, certain assumptions are required to complete 
discounted cash flow economic analyses. These Include: 

4.1 Income tax rates 

4.1.1 Federal 

4.1.2 State 

4.1.3 Local 

4.2 Engineering (and consulting) costs not Included above 

4.3 Cost of capital (to be used as the discount rate) 

All cost estimates should be expressed In terms of probability 
distributions when available and where appropriate. In the absence of the 
full distributions, however, only the principal statistics (mean, median, 
range and/or standard deviation) may be estimated. Point estimates or 
certainty equivalents should be used only when probabilistic estimates are 
unavailable. 
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